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ABSTRACT 



The NPS computer simulation model was modified to study 
the angular dispersion and ring of sputtering sites 
formation of sputtered copper isotopes. The sputtering ratio 
as a function of the ion mass and the microcrystallite size 
needed to contain low-energy events were also investigated. 

A shift in the ejection angles of sputtered Cu (63) 
isotopes toward the surface normal is observed. When the 
"primary knock-on atom" was a light isotope, all of the 
sputtered atoms were ejected more nearly normal to the 
surface. The predominance of sputtered atoms in the <10O> 
and < 1 0 1 > directions from the impact atom indicated an ease 
of momentum transfer in these directions. Results obtained 
from studies based upon a grid of impact points compares 
very favorably with the results obtained whan randomly 
selected impact points were used. The rings of sputtered 
atoms present a "liquid splash" effect. 
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I. 



INTRODUCTION 



The phenomena of physical sputtering have been 
extensively studied since it was first mentioned in 
Philosophical Transactions by Grove in 1352. Sputtering is 
the ejection of atomic material from a target resulting from 
bombardment by an energetic particle. Most scientific 
investigation has been restricted to ionic bombardment 
because the ion can easily be accelerated to any desired 
velocity. Sputtering has defied complete physical 
description because of complicated mechanisms, difficulty in 
treating the behavior of a many-body system, and an 
incomplete understanding of the interacting forces and 
potentials in crystal structures. 

Early experiments typically determined the "sputtering 
ratio" of single and polycrystalline materials. The 
"sputtering ratio" (interchangeable with sputtering yield 
Y) is defined as the number of target atoms ejected from the 
surface per incident ion. 

In 1923, Kingdon and Langmuir [1] bombarded thoriated 
tungsten with various ions in a glow discharge tube. This 
was a special case of sputtering since the thin surface film 
of thorium on a tungsten substrate was sputtered rather than 
the tungsten itself. The results of this experiment 
qualitatively suggested an ejection mechanism. 

A few years later, Blechschmidt and von Hippie [2,3] 
proposed a theory which described sputtering as an 
evaporation of surface atoms (an emission mechanism). 
Earlier, von Hippie [4] had found by spect roscopic analysis 
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that at least some sputtered atoms were in an excited state. 
The sputtering theory showed that atoms in the region of 
impact could rapidly acquire thermal energy if the kinetic 
energy of the incident ion was converted to thermal energy 
at the target surface. If an atom acquired a sufficient 
amount of thermal energy, it would then evaporate from the 
surface, some atoms evaporating while in an excited state. 

Direct application of statistical methods to sputtering 
was made by Harrison [5] who envisioned the interaction of 
two distribution functions (the crystal lattice and the ion 
beam) . These models, based on statistical methods, 
implicitly accept e jection-t y pe mechanisms. 

One of the most, important contributions to the study of 
sputtering was made by Wehner [6] in 1953. He discounted 
the evaporation theory and presented strong evidence for a 
momentum transfer process. Shortly after Wehner's findings 
were reported, Henschke [7] proposed a theory of sputtering 
based solely on classical theory, treating normal and 
oblique incidence sputtering separately. His concept of 
normal incidence sputtering required many-body collisions in 
which the ion was eventually reflected outward to sputter 
surface atoms. This theory was plausible for oblique 
incidence sputtering; however, the case for normal incidence 
sputtering required that the ion be reflected inside the 
crystal . 

Wehner and Rosenberg [8] reported that in the sputtering 
of polycrystalline metals at low ion energy (normal 
incidence) more material is ejected obliquely than in a 
direction normal to the target surface. This is due to the 
fact that less than a full (180°) overall direction change 
of momentum is more favorable for sputtering because this 
requires fewer number of successive collisions for atom 
ejection. In 1960, Wehner and Anderson [9] reported 
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spattering patterns from single crystal targets which showed 
an anisotropic spread of energy from a collision center and 
atoms which are preferentially ejected in close-packed 
directions. 

When an atom collides with another which has a higher 
mass, the scattering angle can vary between 0° and 180°. In 
a collision with a lower mass atom the scattering is 
confined to between 0° and 90° (thus backward reflection 
cannot occur) . In low-energy sputtering of multi-element 
targets. Winters and Sigmund [10] showed theoretically that 
a significant contribution to sputtering, called "single 
reflective collisions", arises from the reflection of lower 
mass atoms from underlying heavier atoms. For a detailed 
discussion of momentum reversal, see Harrison and flagnuson 
[11 ]. 

In 1962, Wehner et al [12] published a compilation of 
sputtering yields for metals and semiconductors in the 
100-600 eV range. In 1966, Wehner, Anderson, and KenKnight 
[13] published the results of a study that covered almost a 
decade of research in low-energy sputtering. A few of the 
numerous topics discussed included more accurate and 
improved techniques for analyzing the sputtered mass, the 
temperature dependence of sputtering yields of varioious 
metals bombarded by atomic and polyatomic ions in the 2.8 to 
10 keV range, the effect of oxygen on ion ejection patterns, 
and the spatial distribution of ejected atoms in metals and 
semiconductors. 

A few low-energy sputtering experiments [14, 15, 16, 17, 
13], which lend themselves to direct comparison with these 
simulations, investigate the behavior of the sputtering 
ratio as a function of ion energy, lattice orientation, 
target material, and ion angle of incidence. 



18 



Paralleling the experimental investigation of the 

sputtering yield and the development of mass spectrographic 

\ 

analysis was the attempt to model the sputtering mechanism 
through the use of computer simulation. In I960, Gibson, 
Goland, Milgram, and Vineyard [19], the first investigators 
to use computer simulation as a means of studying radiation 
damage, set forth the criteria that must be satisfied before 
the simulation method could be applied to- crystals. 
Reference [19] gives a very detailed description of the 
simulation method. Basically, they built a computer model 
to represent metallic copper and studied radiation damage 
events at low and moderate ranges of energy up to 400 aV. 
In their model, one atom in the stationary lattice was given 
an arbitrary kinetic energy and direction of motion, as 
though it had been struck by a bombarding particle, and was 
allowed to interact with the remaining atoms in the lattice. 
This resulted in cascades of independent binary collisions. 
It is well established that sputtering of solids by ion 
bombardment is the result of successive independent 
approximately binary collisions somewhat as if these 
collisions take place with or between separate gas atoms. 
Low-energy ion scattering studies [20] have shown that the 
approximately binary scattering behavior still holds at 
bombarding energies as low as those approaching the 

sputtering thresholds. 

In 1967, Harrison, Levy, Johnson, and Effron [21] used a 
computer to simulate the bombardment of f ace-centered-cubic 
(fee) copper single crystals by an argon ion. From their 
studies, the sputtering mechanisms in the fee crystal were 
isolated and identified. They reported, for ion energies 

less than 10 keV, the sputtering process occurred 

predominately within three atomic layers of the surface. 
The computer simulation included only a repulsive force 
between the atoms. The equations of motion for the 
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many-body system were computed using an "average- force' 1 
algorithm rather than a "central-difference" algorithm, see 
Harrison, Gay, and Effron [22]. In 1972, Harrison, Moore, 
and Holcombe [23] in continuing efforts to develop a more 
realistic computer .model incorporated a Cu-Cu potential 
function which included an attractive potential, surface 
layer relaxation, and surface binding energy for the (100) , 
(110) , and (111) copper planes. They reported improved 
agreement between experimental and simulation results, but 
the original interpretation of the low-energy sputtering 
mechanisms in fee copper crystals remained unchanged. 



Computer simulation of physical processes is still in a 
very infant state because the available processing speed of 
existing computers severely limits the sophistication of the 
models. However, past investigations were quite successful 
in identifying the mechanisms which occur when energetic 
particles interact with crystal atoms. The ability to 
observe these mechanisms is an advantage peculiar to the 
computer simulation because each crystal atom is identified 
by a number for the mathematical calculation, and its 
complete track can be labeled and recorded. The tracks of 
selected atoms can be examined to determine the dynamics of 
the atoms which exit the front surface. 
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II. RESEARCH OBJECTIVES 



A. ANGULAR DISPERSION OF ISOTOPES 

Recently, Wehner [24] reported that the isotope ratio in 
metal deposits obtained by sputtering a flat target at low 
ion energy (lOOeV) under normal incidence is a ^unction of 
the ejection angle. He investigated the bombardment of a 
Cu (63-65) target with 100 eV Hg ions . He reported that the 
lighter isotopes are enriched in a direction normal to the 
target surface. This research was a computer simulation 
investigation of Wehner* s studies of angular dispersion of 
isotopes in a monocrystalline copper lattice. 



B. LIQUID SPLASH EFFECT 

In preliminary NPS simulation studies, after analysis of 
data, a predominance of sputtered atoms occurred at certain 
distances from the center of the lattice representative area 
in which the ion impacts occur. Figures 1 and 2 [ 25 ] depict 
the top layer of a copper fee crystal respectively bombarded 
by neon and gold ions at 2 keV. The numbers represent the 
number of times that the atom was sputtered in 36 events 
(shots) . The dark. triangle in the center represents the 
target impact area. All 36 events are directed into this 
representative area. The darkened elipses highlight the 
atoms most often sputtered and suggest a ring outward from 
the center (especially for the Cu-Au interaction) . 
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Macklin and Metaxas [ 26 ] experimentally investigated the 
splashing of liquid drops on deep and shallow liquid layers. 
They published high speed photographs which showed liquid 
splashes that might be analogous to surface atoms being 
ejected from the surface and a few eventually sputtering. 

This research was an investigation of "ring formation" 
or "liquid splash" effect of sputtered atoms. Determination 
was made of the probability that a congregate of atoms at a 
given distance from the impact area reference will sputter. 



C. ENERGY CONTAINMENT 



As with any computer simulation, a trade-off must be 
made between the accuracy and detail desired and the size 
and time requirements imposed by computer availability. 
Because of this constraint a limit exists on the size of the 
monocrystalline lattice used in the simulation. This 
research was an investigation to determine the size lattice 
required to contain a low energy event. 



D. ION/ATOM MASS RATIO 



In July of 1976, Eer Nisse [27], reported experimental 
confirmation of a monotonic dependence of Y with the atomic 
number. He investigated the yield of a gold target impacted 
by 40 keV ions up to the atomic number of gold. This 
monotonic dependence is predicted by a sputtering theory of 
Sigmund [28]. However, computer simulation data of previous 
NPS studies [29] indicate that the monotonic dependence of 
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III. SIMULATION MODEL MODIFICATIONS 



The basic simulation model used in this research has 
been developed by the NPS group over the last decade [21, 
23, 31, 32]. A brief description of this model and the 
alterations made to it for this research follow. For 
detailed discussions of the simulation method, the reader is 
directed to references 19, 20, 21, and 23. A detailed 
description of the mathematical model has been published 
[ 22 ]. 

In the development of the computer model, the lattice 
unit (LU) proves to be a convenient unit of length, and is 
used extensively. One LU=a 0 /2, where a 0 is the cubic 
lattice constant, the atomic spacing as determined by x-ray 
crystallographic studies. For copper, a 0 =3.615 angstroms. 

Certain properties of the program and results are best 
considered in the framework of a specific run. Each run is 
begun by constructing a target microcr ystallite whose sites 
represent the equilibrium positions of copper atoms near the 
surface of a fee crystal. The target atom - target atom 
interaction is represented by the Born-Mayer Gibson type II 
(Gibson-II) repulsive potential function. 

Each ion of the beam is represented by a single neutral 
atom of the desired element [30]. The ion is placed at the 
point on its trajectory such that its velocity vector 
determines the point of impact. The term bullet and ion are 
used interchangeably throughout this report for the incident 
particle to avoid confusing it with target atoms. 
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The bullet-target potential function was always either 
the Gibson-II or that determined by the Kinetic Secondary 
Electron (KSE-B) Potential [30]. These two potential 
functions seem to represent a reasonable range of "hard" and 
"soft" potential functions. 



As reported by Harrison, Moore, and Holcombe [23], the 
attractive forces between the atoms in the crystal are so 
small, when compared to the forces which occur in even 
moderately energetic collisions, that they may be safely 
neglected in sputtering simulations. Also, in order to use 
the smaller reduced target areas, the ions must be normally 
incident. Even though the computer program still retained 
the capability to utilize attractive potentials and 
obliquely incident ions, they were not used in this 
investigation . 



A. MICROCRYSTALLITE SIZE 



Initial discussion to establish research objectives 
revealed the advantages gained from using a symmetrical 
lattice that could be folded and rotated so that atoms 
sputtered by ions incident at all possible impact points 
would become available for analysis. A lattice with each 
layer a square and having an odd number of atoms along each 
boundary was needed. The lattice generation program as 
previously used could handle this change with no 
alterations. 

Figures 3, 4, and 5 show a roughly isometric projection 
of the labeled atoms in a 9x4x9 (100) orientation, a 9x6x9 

(110) orientation, and a 9x4x9 (111) orientation of a fee 

microcrystallite respectively. This type of projection, 
first used by Harrison and Delaplain [31], with the distance 
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between layers in the Y direction exaggerated, and the atoms 
represented by 20° ellipses has proved to be a valuable tool 
for analyzing the computer programs and for displaying data. 
The first number in each lattice is atom number two, with 
atom number one being the bullet. 

From previous studies, it was clear that at the low 
energy levels considered the sputtering mechanisms were 
confined to the surface layers [21], For this reason, the 
(100) and (111) orientation microcrystallites were only four 
layers deep. The (110) orientation lattice generation 
program constructs a lattice which has significantly fewer 
atoms and its layers are closer together. In an attempt to 
increase the number of atoms in the lattice while keeping 
the required computer time below that of the (100) and (111) 
orientation programs, the (110) orientation raicrocr yst allite 
was set at six layers deep. 



3. CRYSTAL BOUNDARIES 



The previous simulation program contained a 
microcrystallite with established boundaries that were not 
consistent with the development of the potential functions. 
Since the potential function is truncated at the 
nearest- neighbor separation (ROE) , all boundaries (except 
the top surface) were set at ROE. This is justified by the 
fact that an atom outside these boundaries is farther from 
the crystal than the potential truncation radius. The upper 
surface boundary was set by starting a small distance from 
the lattice surface, and increasing the distance in 
succeeding runs until the total yield was independent of the 
distance. For computational efficiency this distance should 
be set as small as possible. 
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Figure 6 shows a drawing of the microcrystallite in 
solid lines and the plane boundaries in dashed lines. The 
boundaries which are defined in the Glossary (Appendix C) 
are shown as planes at the distance from the 
microcrystallite surface cited above. 



C. ISOTOPE GENERATION 



In order to study the angular dispersion of sputtered 
isotopes, a method for generating a microcrystallite with 
the correct isotope ratio and masses was needed. The method 
used assigned an additional variable to each atom that 
identified the atom "type". The atoms were classified 
according tc two types; type 1 - least abundant, heavier 
isotope and type 2 - most abundant, lighter isotope. 

Appendix A is a detailed discussion of how the correct 
isotope ratio was obtained and each atom properly identified 
with the correct mass. 



D. TARGET IMPACT AREAS 

A representative area must contain all possible impact 
points upon a particular surface, that is any point in the 
surface must be projectable into a point of the 
representative area. Figure 7 shows the intrinsic geometry 
of the three fee orientations and the complete sets of 
impact points for each representative area. Figure 8 shows 
the reduced sets of impact points from which the complete 
sets are generated by reflection and rotation. Figures 9, 
10, and 11 show an expanded view of the target impact areas 
and the actual impact points used. The individual impact 
points are those used in previous simulations for normal 
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incidence. This simplification, which greatly reduces the 
computer running time, is only possible when the ion is 
normally incident. 

E. LATTICE FOLDING 

In order to investigate the ring formation or ''liquid 
splash 11 effect, it was necessary to develop a means of 
folding the results obtained from the reduced target area. 
This requirement was satisfied by writing three subroutines; 
KFOLDA, KFOLDB, and KFOLDC. Appendix B gives a detailed 
description of the "folding” logic and the geometry of each 
orientation. Figure 12 shows how the (100) and (110) 
orientation reduced representative areas were folded to 
obtain the complete representative areas. Figure 13 shows 
how the (111) orientation reduced representative area was 
folded and summed to obtain the complete representative 
area. 



F. ATOM SPUTTERING PROBABILITY 

A means of taking the folded data and calculating a 
conditional probabilty that an atom would sputter from a 
particular ring was needed. This was accomplished by 
writing the subroutine KPROB. In this subroutine, the 
distance from each atom to the reference (0,0) point in the 
reduced representative impact area was calculated. then all 
atoms at each distance were checked and the atom sputtering 
count was summed. This sum was divided by the total number 
of sputs from all the shots. This quotient was the 
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conditional probabilty that a given ring at a known distance 
would sputter an atom given the atoms that had already 
sputtered . 

KPR03 performs the above calculations. However, a very 
interesting result is obtained by dividing the conditional 
probabilty obtained above by the number of atoms in a 
particular ring. This result was a conditional probability 
that an atom at a given distance from the impact point would 
sputter. 

G. EJECTION ANGLES 



The ejection angle is defined as the angle between the 
atom's velocity vector and the surface normal. A tally was 
kept of each atom's ejection angle, grouped in five degree 
intervals. This tally was recorded separately for each 
isotope. To provide a quick analysis of the data, a program 
(HISTO) was written to print out a histogram of the number 
of sputtered atoms ejected in each interval for each 
isotope . 



H. VARIABLE-SHOT SELECTION 



The previous program was only flex 
either a single-shot run or all 36 
the lattice and the ion's energy 
computer run time for a set of 36 sho 
(12 - 18 hours) . At this point, a 
"package" was developed. The feature 
to specify the number of shots desir 
impact point of the "package". The 
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particular package run was printed out. The packages had to 
be combined by hand calculations in order to get the total 
information for a particular run. Even so, this procedure 
allowed runs to be made with larger lattices and higher 
energies (i.e. the containment studies) that previously 
would have been prohibitively long. 



I. RANDOMLY SELECTED IMPACT POINTS 



To obtain more angular dispersion data for the 
low-energy, (100) orientation events, a modification was 
made to the simulation model. This modification replaced 
the systematic method of incrementing the impact point over 
a grid of points in the representative target area with a 
random-shot capability. Again, the IMSL GGUB routine was 
used to generate uniformly distributed, psuedo-random 
numbers. The first random number was labeled the x 
coordinate in the impact area. The second random number was 
labeled the z coordinate and tested to determine if it was 
inside the impact area z boundaries determined by the value 
of the x coordinate. If this test was passed sucessfully, 
these values for the x and z coordinates were used. 
However, if the test did not pass (the x and z coordinates 
were outside the impact area) , the two original numbers were 
discarded, two more were selected, and the process was 
repeated again. 

This random impact point did not affect the variable - 
shot selection feature. But, enough random numbers had to 
be generated to allow for those points that did not fall 
into the impact area. 
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The results of randomly bombarding the impact area will 
be discussed later. This feature was incorporated in 
response to other investigators suggestion that the grid of 
impact points be compared with a randomly bombarded impact 
area. 
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IV. RESULTS 



The three lattice orientation models for the copper 
target were utilized in computer simulation runs in which 
100-eV, 200-eV, 500-eV, and 1-keV normally incident ions 
bombarded the target. The ion masses varied from 4 amu 
(He + ) to over 200 amu (Hg + ) . The lattice was constructed 
with copper atoms of either the natural mass (63.54 amu) or 
the proper ratio of isotopes (63.93 amu and 64.93 amu) . The 
lattice size varied from a 9x4x9, (100) orientation of 162 
atoms to a 17x4x17, (100) orientation of 578 atoms. 

Computer runs were made for various permutations of ion 
energy, ion mass, lattice orientation, lattice size, and 
potential function. The summary of each run contained the 
number of atoms sputtered from each impact point, the number 
of times each atom sputtered, a histogram showing the 
angular dispersion of the sputtered atoms for each isotope, 
the number of sputs for each atom after the folding process, 
the distance from each ring to the impact area reference 
point, the number, of atoms and the sum of their sputs in 
each ring, and the probability that each ring contained a 
sputtered atom. 

In the sputtering summaries the numbers in the ellipses 
indicate the number of times that particular atom sputtered 
in the simulation. The dark areas, a triangle for the (100) 
and (111) orientations and a rectangle for the (110) 
orientation, indicate the location of the target impact 
area. The atom that is struck by the ion is called the 
"primary knock-on atom" (hereafter PKA) . Since different 
seeds were used to generate the random isotope sites, a 
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letter L or H denotes whether the sputtered atom was a 
light or heavy isotope. Wherever possible, the results for 
the Gibson-II and KSE-B potential functions were combined. 
In the figures displaying the sputtering summary and the 
mechanism tracing, some of the lower layers were omitted. 
However, all atoms were included in all calculations. In 
the 100-eV results all the sputtered atoms are shown. 
Wherever a direction in the lattice is specified, it is 
always with respect to that particular lattice axes. The 
top surface is always an XZ plane. 



A. 9X9 HICROCRYSTALLITE 



Initially, all simulations used a 9x9 
determine the sputtering ratio as a function of 
at 100 eV and 1 keV. The model was modified to 
isotope generator, the symmetrical folding of 
the angular dispersion of sputtered isotopes, an 
sputtering" probabilities. 
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The (100) orientation 
constructed in a lattice consisting 
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(1) H el ium Io ns . The results when 100-eV 
helium icns impact the target area are shown in figures 15, 
16, and 17. Figure 15 shows the sputtering summaries using 

(a) Gibson-II, resulting in a sputtering ratio of 0.444, and 

(b) KSE-B, resulting in a sputtering ratio of 0.111. The 
predominance of the sputtered atoms in the <100> and < 1 0 1 > 
directions is quite evident. Figure 16 shows the 
conditional probability that a ring of atoms, as a function 
of the ring distance from the impact area reference point., 
will contain an atom that sputtered for (a) Gibson-II and 
(b) KSE-B. The term "ring sputtering" probability, 
henceforth ring probability, was used as a label to mean the 
above conditional probability. The relative amplitudes of 
the ring probabilties could be misleading because the 
calulaticn does not take into effect the number of atoms in 
each ring. Again, the predominance of rings that contain 
atoms in the <100> and < 1 0 1 > directions is evident. The 
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peak of the KSE-3 curve was closer to the center than the 
Gibson-II peaks, which were spread out all the way to the 
lattice edges. Figure 17 shows the angular dispersion 
results for (a) Gibson-II and (b) KSE-3. The particular 
seed used in the psuedo-r andom number generator to locate 
the isotope random sites was such that only light isotopes 
sputtered in this case. 

(2) Neon Ions . The results when 100-eV neon 
ions impact the target area are shown in figures 18, 19, and 
20. Figure 19 shows that for the light mass of the neon 
ion, atoms on the lattice boundary are sputtered using (a) 
Gibson-II and (b) KSE-3. Gibson-II resulted in a sputtering 
ratio of 0.889, while the sputtering ratio for KSE-B was 
0.806. Figure 19 shows that the Ne + sputtered atoms were 
further away from the lattice center than the He + , for both 
(a) Gibson-II and (b) KSE-B. This was confirmed by the 
increased number of next-nearest neighbor atoms (in the 
< 1 0 1 > direction) sputtering. Figure 20 indicates that some 
of the heavier isotopes sputtered for (a) Gibson-II and (b) 
KSE-B. There were not enough sputtered atoms to notice any 
trend in the isotope angular dispersion. 

(3) A rgon Ions^ The argon simulations w a re 
used to investigate the ’’reflective collision” theory. 
Argon ions of 100-eV energy bombarded the target area using 
both of the potential functions. In the random isotope site 
locations, the PKA (atom 22) was a heavy isotope, while 
three (52, 58, and 67) of the four atoms directly underneath 
were light isotopes. Four simulations were run using each 
of the potential functions; (1) the random isotope sites, 
(2) forcing the PKA to a light isotope, (3) forcing the PKA 
to a heavy isotope and all four atoms underneath to a light 
isotope, and (4) forcing the PKA to a light isotope and all 
four atoms underneath to a heavy isotope. The remainder of 
the lattice isotopes were unchanged. 
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Figure 21 shows the sputtering summaries 
for all four combinations using (a) Gibson- II and (b) KSE-B. 
In all four cases, identical atoms sputtered independent of 
isotope locations. Gibson-II resulted in sputtering ratios 
of 1.056, while KSE-B resulted in sputtering ratios of 
0.944. Figure 22 illustrates the ring probabilty in all 
four cases for (a) Gibson-II and (b) KSE-B. Figure 23 
illustrates the angular dispersion of ejected atoms, using 
Gibson-II, for (a) case (1) , (b) case (2) , (c) case (3) , and 
(d) case (4) . No light isotopes were sputtered. Figure 24 
illustrates the same results for KSE-B. 

Since at least four atoms were ejected into 
a different five degree window for Gibson-II, the sputtering 
mechanisms were checked. The most active event occurred 
when the impact area was struck by the ion at impact point 
4060 (impact displaced 0.40 LU in the X direction and 0.60 
LO in the Z direction from the center of the PKA, hence 
4060) . Four atoms were sputtered at this point. A 
single- impact program was run at this point with the output 
listing the positions, velocities, kinetic energy, and 
potential energy of each atom with a total energy greater 
than 0.1 eV. The listing was printed once every five 

timesteps. A timestep is the time interval between dynamic 
calculations. It is variable and determined as the time 
required for the fastest moving atom to move 0.1 LO from its 
previous position. In real time the timestep is on the 
order of 10 -15 seconds. Figure 25 illustrates atom 
trajectories in the sputtering of atoms 17, 32, 57, and 77. 
In cases (3) and (4) , the sputtering mechanisms were the 
same displaying only a slight difference in the ejection 
angles and the velocities of the sputtered atoms. 
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The argon ion collides with atom 22 (PKA) 

and atom 27 by timestep 5 , pushing them down and under atoms 

17 and 32 respectively. The ion continued downward and 
collided with atom 67 by timestep 15. The ion reflected 
back up through the hole vacated by atoms 22 and 27. By 
timestep 15, atoms 17 and 32 were moving upward. Atom 67 
moved down colliding with atoms 107 and 112 by timestep 20. 
Atoms 107 and 112 squeezed [21] under atoms 102 and 117 
respectively. Atoms 102 and 117 were moving up by timestep 
30, and struck atoms 57 and 77 from below by timestep 40. 

By this time, atoms 17 and 32 had moved far enough above the 

lattice with sufficient upward velocity to be classified as 
M sputs' T . Atom 17 was ejected at an angle of 33.56° to the 
surface normal with an upward velocity of 2924 m/sec for 
case (3) , and at an angle of 33.36° with a velocity of 2974 
m/sec for case (4) . Atom 32 was ejected at an angle of 

33.65° with an upward velocity of 2926 m/sec for case (3), 

and at an angle of 33.56° with a velocity of 2918 m/sec for 
case (4) . Atoms 57 and 77 continued upward through the 
holes left by the sputtered atoms 17 and 32. As atom 57 
passed through this hole its trajectory was bent closer to 
the surface normal by atoms 12, 16, and 21. Likewise, atom 
77 was reflected upward by atoms 28, 33, and 37. By 

timestep 75, atoms 57 and 77 were far enough above the 
lattice surface with sufficient upward velocity to be 
classified as "sputs". Atom 57 was ejected at an angle of 

18.12° with an upward velocity of 2967 m/sec for case (3), 

and at an angle of 17.75° with a velocity of 2882 m/sec for 
case (4) . Atom 77 sputtered at an angle of 25.04° with an 
upward velocity of 2799 m/sec for case (3) , and at an angle 
of 24.75° with a velocity of 2859 m/sec for case (4). 
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It is noted that when the PKA was a light 
isotope, all of the sputtered atoms were ejected closer to 
the surface normal. Also noted is the remarkable symmetry 
of ejection angles and velocities of the sputtered atoms. 

(4) Cop p er Ions^ The results when 100-eV 
copper ions impact the target area are shown in figures 26, 
27, and 28. Figure 26 shows the sputtering summaries for 
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potential functions. These results further indicate a 
maximum sputtering ratio occurs near an ion/atom mass ratio 
of unity. Figure 27 illustrates a continued spread outward 
of the ring probabilty for (a) Gibson-II and (b) KSE-3. 
Figure 28 illustrates the angular dispersion results for (a) 
Gibson-II and (b) KSE-B. 

(5) Kr yp ton Ions_. The results when 100-eV 
krypton ions impact the target area are shown in figures 29, 
30, and 31. Figure 29 shows the sputtering summaries for 

(a) Gibson-II, resulting in a sputtering ratio of 0.694, and 

(b) KSE-B, resulting in a sputtering ratio of 1.139. Here, 
a maximum in the sputtering ratio was reached for KSE-B, 
while the sputtering ratio for Gibson-II dropped 
drastically. Figure 30 illustrates a continued outward 
spread of the ring probabilities for (a) Gibson-II and (b) 
KSE-B. It is noticed that the outward spread still occured, 
even though the Gibson-II yield decreased. Figure 31 
illustrates the angular dispersion results for (a) Gibson-II 
and (b) KSE-B. 
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(6) Xenon Ions^ The results when 100-eV xenon 
ions impact the target area are shown in figures 32, 33, and 
34. Figure 32 shows the sputtering summaries for (a) 
Gibson-II, resulting in a sputtering ratio of 0.556, and (b) 
KSE-B, resulting in a sputtering ratio of 1.028. Figure 33 
illustrates the ring probabilities for (a) Gibson-II and (b) 
KSE-B. The outward spread of the ring probability peaks is 
still seen for KSE-B, while Gibson-II remained constant. 
Figure 34 illustrates the angular dispersion results for (a) 
Gibson-II and (b) KSE-B. 

(7) Gold Ions^ T he results when 100-eV gold 
ions impact the target area are shown in figures 35, 36, and 
37. Figure 35 shows the sputtering summaries for (a) 
Gibson-II, resulting in a sputtering ratio of 0.528, and (b) 
KSE-B, resulting in a sputtering ratio of 0.778. Figure 36 
illustrates the ring probabilities for (a) Gibson-II and (b) 
KSE-B. Gibson-II has settled out with approximately egual 
probability of sputtering from atoms in the first two rings 
in the < 1 0 1 > direction. In the case of KSE-B, a shift from 
the first and second rings in the <100> direction toward the 
second ring in the <101> direction is noticed. Figure 37 
illustrates the angular dispersion results for (a) Gibson-II 
and (b) KSE-B. 

(8) Mercury Ions . These Hg + bombardment 
results are compiled from the sets of simulation runs used 
to investigate the ''reflective collision" theory. Figure 38 
shows the sputtering summaries for all four cases, in which 
the masses of the PKA and the four atoms underneath are 
varied. The results for (a) Gibson-II were identical in all 
four cases, resulting in a sputtering ratio of 0.528. The 
sputtering ratio was 0.743 using KSE-B in a lattice 
consisting of atoms with only the natural mass. The results 
for (b) KSE-3 depended on the order of the heavy over light 
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interaction. When the PK A was a heavy isotope, shown in 
figure 38(b), the sputtering ratio was 0.7229 However, when 
the PKA was a light isotope, atom 30 sputtered an additional 
time (denoted by ▲ ) and the sputtering ratio became 0.750. 
Figure 39 illustrates the ring probabilties for (a) 
Gibson-II, which remained relatively unchanged, and (b) 
KSE-B, which still displayed a shift to the second ring in 
the 101 direction. Figure 40 illustrates the angular 
dispersion results for Gibson-II. Figure 41 shows the 
angular dispersion results for KSE-B. Again, atom 30 
sputtered an additional time. The additional ejection angle 
is denoted by A . 



b. 1-keV Results 

The 1-keV events were obviously not contained in 
the 9x4x9 microcrystallite. However, the simulations were 
run to compare the sputtering ratios with experimental 
results. Also, the dependence of the sputtering on ion mass 
was investigated. Useful information can be obtained about 
the ring probabilties and isotope angular dipersion at this 
higher energy. The 162 atom, 9x4x9 microcrystallite, with 
the correct ratio of copper isotopes at random sites, was 
used. Even though all atoms were used in the calculations 
and the printout provided information on all atoms, only the 
top layer results are shown in the figures. Even in the 
worse case (Cu/Ar + , Gibson-II), more than 86 per cent of 
the sputtered atoms were from the top layer. 

The sputtering simulations for normally 
incident, 1-keV ions were run using He + , Ne + , Ar + , Cu + , 
Kr + , Xe + , Au + , and Hg + ions. Figure 42 shows the sputtering 
ratio as a function of the ion mass for both potential 
functions. As expected, the sputtering ratios have 
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increased. However, they are no longer smooth and 
continuous, but exhibit a randomness, as if each was 
independent of the others. 

All the 1-keV simulations used the same random 
isotope sites. The letters L and H, denoting light or heavy 
isotopes, appear only in the Cu/Ne + , KSE-B summary (more 
different atoms sputtered - but not the highest sputtering 
ratio) . 



(1) Ions • The results when 1-keV helium 
ions impact the target area are shown in figures 43, 44, and 
45. Figure 43 shows the sputtering summaries for (a) 
Gibson-II, resulting in a sputtering ratio of 1.914, and (b) 
KSE-B, resulting in a sputtering ratio of 2.833. In both 
cases, boundary atoms sputtered. Figure 44 illustrates the 
ring probabilities for (a) Gibson-II and (b) KSE-B. Both 
display a large action-reaction effect, as if a bomb were 
dropped at the center and everything out from the center was 
blasted upward. Figure 45 illustrates the angular 
dispersion results for (a) Gibson-II and (b) KSE-B. 

(2) Neon Ions . The results when 1-keV neon ions 
impact the target area are shown in figures 46, 47, and 48. 
Figure 46 shows the sputtering summaries for (a) Gibson-II, 
resulting in a sputtering ratio of 3.686, and (b) KSE-B, 
resulting in a sputtering ratio of 4.500. Host of the 
surface atoms were sputtered, indicating that the energy ( 
or momentum ) spread out over the entire lattice surface. 
Figure 47 illustrates the ring probabilities for (a) 
Gibson-II and (b) KSE-B. The KSE-3 simulation indicates a 
spreading outward of the sputtered atoms. The Gibson-II 
simulation not only exhibits increased sputtering of atoms 
further from the center, but it also clearly shows a shift 
in the number of sputtered atoms from those at 2.83 LU in 
the < 1 0 1 > direction to those at 3.16 LU in the <103> 
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direction. As mentioned earlier, this could be artificial 
because the ring at 3.16 LU is the first ring to contain 
eight atoms. Figure 48 illustrates the angular dispersion 
results for (a) Gibson-II and (b) KSS-B. 

(3) Arg on Ions^ The results when 1-keV argon 

ions impact the target area are shown in figures 49, 50, and 
51 . Figure 49 shows the sputtering summaries for (a) 
Gibson-II, resulting in a sputtering ratio of 3.806, and (b) 
KSE-3, resulting in a sputtering ratio of 4.611. Figure 50 
illustrates the ring probabilities for (a) Gibson-II and (b) 
KSE-B. The Gibson-II simulation continued to increase in 
the outward direction at the expense of the atoms at 3.16 LU 
in the < 1 0 3 > direction. However, the KSE-3 simulation 
reversed its trend and the ring probability for the close-in 
rings (1.41 LU and 2.00 LU) increased at the expense of the 

ring at 2.83 LU and the rings further out from the center. 

Figure 51 illustrates the angular dispersion results for (a) 
Gibson-II and KSE-3. 

(4) Cop pe r Ions^ The results when 1-keV copper 

ions impact the target area are shown in figures 52, 53, and 
54. Figure 52 shows the sputtering summaries for (a) 
Gibson-II, resulting in a sputtering ratio of 3.028, and (b) 

KSE-3, resulting in a sputtering ratio of 3.972. These low 

sputtering ratios not only disagree with a monatonic 

increase in the yield with mass, but are drastically low for 
the maximum to occur at a ion/atom mass ratio of unity. 
Figure 53 illustrates the ring probabilities for (a) 
Gibson-II and (b) KSE-B. The Gibson-II simulation showed an 
ubexpected change. The ring probability increased for the 

rings at 2.00 LU (<100> direction) and 3.16 LU (< 1 03 > 

direction). This trend was also displayed in the KSE-B 
simulation, where it starts to occur in the Cu/Ar + 

simulation. Figure 54 illustrates the angular dispersion 
results for (a) Gibson-II and (b) KSE-B. 
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The results when 1-keV 



(5) Kr££ton Ions_. 
krypton ions impact the target area are shown in figures 55, 
56, and 57. Figure 55 shows the sputtering summaries for 

(a) Gibson-II, resulting in a sputtering ratio of 2.833, and 

(b) KSE-B, resulting in a sputtering ratio of 3.917. Note 
that the sputtering ratio, for both potential functions, 
continues to decrease with increased ion mass. Figure 56 
illustrates the ring probabilities for (a) Gibson-II and (b) 
KSE-B. Eoth potential function simulations continue to 
exhibit the trend to establish a peak at 2.00 LU in the 
<100> direction and at 3.16 LU in the <103> direction. 
However, the KSE-3 simulation shows significantly larger 
ring probability at 3.16 LU. Figure 57 illustrates the 
angular dispersion results for (a) Gibson-II and (b) KSE-B. 

(^) Xenon Ion s^ The results when 1-keV xenon 
ions impact the target area are shown in figures 58, 59, and 
60. Figure 58 shows the sputtering summaries for (a) 
Gibson-II, resulting in a sputtering ratio of 2.333, and (b) 
KSE-B, resulting in a sputtering ratio of 3.472. 3oth 
potential functions continued to show a decrease in the 
sputtering ratio with an increase in the mass of the ion. 
Figure 59 illustrates the ring probabilties for (a) 
Gibson-II and (b) KSE-B. All the previous trends, except an 
increased ring probability at 3.16 LU were shattered. 
Figure 60 illustrates the angular dispersion results for (a) 
Gibson-II and (b) KSE-B. 

(7) Gold Ionsj. The results when 1-keV gold 
ions impact the target area are shown in figures 61, 62, and 
63. Figure 61 shows the sputtering summaries for (a) 
Gibson-II, resulting in a sputtering ratio of 2.778, and (b) 
KSE-3, resulting in a sputtering ratio of 4.167. A definite 
increase in the sputtering ratio occured. Figure 62 
illustrates the ring probabilities for (a) Gibson-II and (b) 



KSE-B. No new trends were noted, except that the ring at 
3.16 LU continues to dominate and the other ring 
probabilities fluctuate. Figure 63 illustrates the angular 
dispersion results for (a) Gibson-II and (b) KSE-B. 

' (8) Mercury Io ns . The results when 1-keV 
mercury ions impact the target area are shown in figures 64, 
65, and 66. Figure 64 shows the sputtering summaries for 

(a) Gibson-II, resulting in a sputtering ratio of 2.778, and 

(b) KSE-B, resulting in a sputtering ratio of 4.278. These 
simulations resulted in no change in the sputtering ratio 
for Gibson-II, and only a slight increase for KSE-B. Figure 
65 illustrates the ring probabilities for (a) Gibson-II and 
(b) KSE-B. Only insignificant changes occured compared to 
the gold ion results. Figure 66 illustrates the angular 
dispersion results for (a) Gibson-II and (b) KSE-B. 

c. 2-keV Results 



The 2-keV ion simulations were the earliest runs 
made. No isotope information was obtained. As a 
consequence, the only results are the sputtering ratio 
versus ion mass curve, figure 67, for both potential 
functions. They exhibit an even more erratic behavior than 
did the 1-keV simulations. The sputtering summaries 
frequently show atoms sputtered from the third layer, and 
rarely an atom sputtered from the fourth layer. At certain 
impact points as many as 25 atoms sputtered. The small, 
9x4x9 microcrystallite was effectively blasted apart by the 
2-keV ions. 
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2 • Qr ientat io n 1.0 0 SLl 

The effort devoted to (110) orientation simulation 
runs was very limited due to the higher priority assigned to 
the (100) orientation runs. As a result, only the argon and 
mercury ion simulations at 100 eV were run. Still needing 
investigation are the sputtering ratio as a function of ion 
mass and sputtering mechanism traces. 

The (110) orientation m icrocrystall ite was 
constructed in a lattice consisting of 123 copper isotopes. 
In the random isotope site locations, the PKA (14) was a 
light isotope, as were three (33, 36, and 37) of the four 
atoms underneath the PKA. Figure 4 shows an isometric 
projection, of the lattice and its numbering scheme. Figure 
10 shows an expanded view of the grid of 36 impact points 
used. A 9x6x9 lattice was used with no increase in computer 
time required to complete a run. 

Two sets of simulations were run to investigate the 
"reflective collision" theory. Four runs were made in each 
set of simulations for each potential function. The impact 
area was bombarded by 100-eV, argon and mercury ions. The 
four runs consisted of; (1) the random isotope sites (PKA 
was a light isotope) , (2) forcing the PKA to a heavy 
isotope, (3) forcing the PKA to a light isotope and all four 
atoms underneath to a heavy isotope, and (4) forcing the PKA 
to a heavy isotope and all four atoms underneath to a light 
isotope. The remainder of the lattice isotopes were 
unchanged. 
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a . 



Argon Ions 



The results when 100-eV argon ions impact the 
target area are shown in figures 68-75. Figure 63 shows the 
sputtering summaries using Gibson-II for all four cases. 
Case (1) and case (3) , shown in figure 63 (a) , resulted in a 
sputtering ratio of 0.839. Case (2) and case (4) , shown in 
figure 68(b), resulted in a sputtering ratio of 0.861. A 
predominance of sputtered atoms in the < 1 0 0 > and <101> 
directions is evident. Figure 69 illustrates the ring 

probabilities for (a) case (1) and case (3) , and (b) case 

(2) and case (4) using Gibson-II. Notice that an atom at 
the boundary ( in the < 1 0 1 > direction ) sputtered. Figure 
70 illustrates the angular dispersion results for (a) case 
(1) , and (b) case (3) using Gibson-II. Figure 71 shows the 
same information for (a) case (2) , and (b) case (4) . 

Figure 72 shows the sputtering summaries using 
KSE-B for all four cases. Case (3) resulted in a sputtering 
ratio of 0.750, while case (1) sputtered atom 26 one 
additional time (denoted by ▲ ) and had a sputtering ratio 
of 0.773. Case (1) and case (3) are shown in figure 72(a) . 
Case (2) resulted in a sputtering ratio of 0.778, while case 
(4) sputtered atom 8 one additional time and had a 
sputtering ratio of 0.805. Case (2) and case (4) are shown 
in figure 72(b) . Again, it is evident that only atoms in 
the <100> and < 1 1 0 > directions sputtered. Figure 73 
illustrates the ring probabilities, using KSS-B, for (a) 
case (1) and case (3) , and (b) case (2) and case (4) . The 
close-in trend, similar to that of Gibson-II, is evident, as 
are the sputtered atoms at the boundary in the < 1 0 1 > 
direction. Figure 74 illustrates the angular dispersion 
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results, using KSE-B, for (a) case (1) , and (b) case (3). 

Figure 75 shows the same information for (a) case (2) , and 
(b) case (4) . 



b. Mercury Ions 



The results when 100-eV mercury ions impact the 
target area are shown in figure 76, 77, and 78. Figure 
76(a) shows the sputtering summaries for all four cases, 
using Gibson-II. Case (1), case (2), and case (4) resulted 
in a sputtering ratio of 0.278. Case (3) sputtered atom 13 
one additional time (denoted by ▲) and had a sputtering 
ratio of 3.05. Figure 76(b) shows the sputtering summaries 
for all four cases, using KSE-B, which resulted in 
sputtering ratios of 0.222. Figure 77 illustrates the ring 
probabilities for (a) Gibson-II and (b) KSE-B. Figure 78 
illustrates the angular dispersion results for (a) Gibson-II 
and (b) KSE-B. 



No significant trends were noted in the (110) 
orientation. However, the ease with which atoms sputtered 
in the < 1 0 1 > direction is evident. 

3 . (Ill) Or ienta t io n 100e V 

Like the (110) orientation, the (111) orientation 
did not receive extensive attention. Both the sputtering 
ratio as a function of the ion mass, and sputtering 
mechanism traces need further investigation. 

The (HI) orientation microcrystallite is 
constructed in a 9x4x9 lattice consisting of 163 atoms. In 
the random isotope site locations, the PKA (22) was a light 
isotope, as were three (58, 63, and 67) of the four atoms 
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underneath the PKA. Figure 5 shows an isometric projection 
of the lattice and its numbering scheme. Figure 11 shows an 
expanded view of the grid of 21 impact points used. 

Two sets of simulations were run to investigate the 
"reflective collision" theory. Four runs were made in each 
set of simulations for each potential function. The target 
was bombarded by 100-eV argon and mercury ions. The four 
runs consisted of; (1) the random isotope sites (PKA was a 
light isotope) , (2) forcing the PKA to a heavy isotope, (3) 
forcing the PKA to a light isotope and all four atoms 
underneath to a heavy isotope, and (4) forcing the PKA to a 
heavy isotope and all four atoms underneath to a light 
isotope. The remainder of the lattice isotopes were 
unchanged . 



a. Argon Ions 

The results when 100-eV argon ions impact the 
target area are shown in figures 79-83. Figure 79 shows the 
sputtering summaries for (a) Gibson-II, resulting in 
sputtering ratios of 0.857, and (b) KSE-B, resulting in 

sputtering ratios of 0.762. For KSE-B, case (2), atom 16 
sputtered one additional time (denoted by A ) , resulting in 
a sputtering ratio of 0.310. Figure 80 illustrates the ring 
probabilities for (a) Gibson-II and (b) KSE-B. It was 

unexpected that, for KSE-B, all of the sputtered atoms are 

in the < 1 0 1 > direction in the ring at 2.83 LU. Figure 81 
illustrates the Gibson-II angular dispersion results for (a) 
case (1) and case (3) , and (b) case (2) and case (4) . 

Figure 82 illustrates the KSE-B angular dispersion results 
for (a) case (1) and (b) case (3) . Figure 83 shows the same 
results for (a) case (2) , and (b) case (4) . 
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b. Mercury Ions 



The results when 100-eV Mercury ions impact the 
target area are shown in figures 84-89. Figure 84(a) shows 
the sputtering summaries, using Gibson-II, for case (1) , 
case (2) , and case (4) , which resulted in sputtering ratios 
of 0.476. Case (3) had a sputtering ratio of 0.381 and is 
shown in figure 84(b) . Figure 85 shows the sputtering 
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85(b), have sputtering ratios of 0.952. Figure 36 

illustrates the ring probabilities for (a) Gibson-II and (b) 

KSE-B . Again, it is evident that the only sputtered atoms 

are in the ring at 2.83 LU, in the < 1 0 1 > direction. Figure 

87 illustrates the Gibson-II angular dispersion results for 

(a) case (1) , case (2) , and case (4) , and (b) case (3) . 

Figure 88 illustrates the KSE-B angular dispersion results 

for (a) case (1) and (b) case (3) . Figure 89 shows the same 

/ 

information for (a) case (2) , and (b) case (4) . 
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mercury ions bombarding a 13x6x13, (110) orientation 
microcrystallite containing 255 atoms are shown in figure 
90. Figure 90(a) is for Gibson-II, while figure 90(b) is 
for KSE-B. Comparison of figure 90 with figure 76, 
confirmed that the 9x5x9 microcrystallite did not contain 
the 100-eV, mercury event. Based on the original definition 
of containment, the 13x6x13 microcrystallite did not contain 
the Gibson-II simulation run. 

In the (111) orientation, the sputtering summaries for 
mercury ions bombarding a 13x4x13 microcrystallite 
containing 339 atoms are shown in figure 91. Figure 91(a) 
is for Gibson-II, while figure 91 (b) is for KSE-B. 
Comparison of figure 91 with figures 84 and 85, indicate 
that even though the 9x4x9 lattice had no boundary atoms 
sputtered, the sputtering ratio increases when a 13x4x13 
lattice is used. 

1 • 1 00 - eV Conta inme nt 

Based on the above information, the definition for 
"containment" was changed to: "an event is contained when an 
increase in the microcrystallite size does not result in any 
further increase in the yield." Earlier simulations 
resulted in larger yields for KSE-B. For this reason, KSE-B 
was used to investigate the size microcrystallite needed to 
contain a 100-eV event. 

The sputtering ratio as a function of the ion mass 
was examined for the (100) orientation . Figure 9 2 sho ws the 
results for 9x4x9, 11x4x11, and 11x6x11 microcrystallites. 
Additionally, simulations were run for copper (largest 
sputtering ratio) and mercury ions for 13x4x13, 13x6x13, and 
17x4x17 microcrystallites. Based on the new definition of 
containment, 100-eV helium and xenon events are contained in 
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a 9x4x9 microcrystallite because an 11x4x11 microcrystallite 
produces the same yield. However, the argon, copper, and 
krypton events were not contained. Further increase of the 
mi'crocrystallite size demonstrated that the 100-eV, copper 
event is contained by a 13x4x13 microcrystallite. Figure 92 
also illustrates that the maximum sputtering ratio occurs 
when the ion/atom mass ratio is approximately unity. Figure 
93 shows the same information in a plot of sputtering ratio 
versus lattice size for all eight ions. 

2 • 1 ~x! xl 7 Mi cr o c rystallite 

Asuming that all 100-eV events are contained by a 
17x4x17 microcrystallite, the effect of increasing the ion 
energy was investigated. The reason for this investigation 
was to determine how energetic an event could be contained 
by a 17x4x17 microcrystallite. And to compare higher energy 
events with experimental results. The reason for not 
increasing the lattice size without limit was obviously to 
keep the computer run time within reason. Again, only KSE-B 
with mercury and copper ions were investigated. 



a. Copper Ions 



The results when 100-eV copper ions impact the 
target area are shown in figures 94, 95, and 96. Figure 94 
shows the sputtering summary with a sputtering ratio of 
1.389. Containment for the 100-eV copper event had already 
been established. Figure 95 illustrates the "atom 
sputtering probability”, henceforth called the 
"probabilit y ” . The ’’probability” is the "ring sputtering 
probability” divided by the number of aroms in each ring. 
The dimensions for "probability" is sputs per atom divided 
by the total sputs at a fixed distance from the impact area 
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reference point. Figure 95 indicates an enhancement, of 
sputtered atoms in the <100> and <101> directions. Figure 96 
illustrates the angular dispersion results. 

Results for the 200-eV copper ion simulations 
are shown in figures 97, 98, and 99. Figure 97, the 
sputtering summary with a sputtering ratio of 2.889, shows 
that boundary atoms are sputtered, and that the event 
definitely is not contained. Figure 98 illustrates the 
"probability'’ results and indicates that an atom farther out 
had increased chances of sputtering. Figure 99 illustrates 
the angular dispersion results. A bias toward light 
isotopes closer to the normal is evident. 

The results for the 500-eV copper ion 
simulations, with a sputtering ratio of 4.028, are shown in 
figures 100, 101, and 102. Again, boundary atoms are 
sputtered and a slight shift outward of the "probability" is 
noticed. Figure 102 illustrates the angular dispersion 
results which showed that more light isotopes are ejected 
more nearly normal to the lattice surface. 



the 1-keV, copper ion results with a sputtering 
ratio of 7.083, are shown in figures 103, 104, and 105. 

Figure 103, the sputtering summary, shows a very large 
number of the surface atoms sputtered. Boundary atoms in 
every direction are sputtered. This increased spreading 
outward was confirmed in figure 104, which shows the 
"probabilities" relatively level. Figure 105 illustrates 
the angular dispersion results. 



Obviously, the 200-eV, 500-eV, and 1-keV copper 
events are not contained, even in a 17x4x17 
micro cry stallite . 
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b. 



Mercury Ions 



The results for the 100-eV mercury ion 

simulations are shown in figures 106, 107, and 108. Again, 

the 11x4x11 lattice has already been determined to contain 
the 100-eV, mercury event. Figure 106 shows the sputtering 
summary with a sputtering ratio of 0.833. Figure 107 

illustrates the "probability", and like the copper 

\ 

simulation, was predominant in the < 1 0 1 > direction. Figure 
108 illustrates the angular dispersion results. 

The 200-eV mercury ion results, with a 

sputtering ratio of 1.722, are shown in figures 109, 110, 

and 111. Figure 109, the sputtering summary, shows an 
outward spreading of the sputtered atoms. However, no 

boundary atoms were sputtered and the 200-eV, mercury event 
is tentatively classified as "contained". Figure 110 

confirms the outward shift in "probability". Figure 111 
illustrates the angular dispersion results, which indicate a 
greater number of heavy isotopes sputtering than the isotope 
ratio would justify. 

The 500-eV mercury ion results, with a 

sputtering ratio of 3.667, are shown in figures 112, 113, 

and 114. Figure 112, the sputtering summary, indicates that 
boundary atoms are sputtered. Figure 113 illustrates the 
continued outward spreading of sputtered atoms. Figure 114 
illustrates the angular dispersion results, and the number 
of heavy to light isotopes closer to the correct ratio. 

The 1-keV mercury ion results, with a sputtering 
ratio of 6.139, are shown in figures 115, 115, and 117. 

Figure 115 shows that boundary atoms are sputtered. Figure 

116 shows the "probability" still spreading outward. Figure 

117 illustrates the angular dispersion results. 
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The "probability" figures for both the copper 
and mercury ions are remarkably similar and display the same 
trends. 

c. Sputtering Mechanism Traces 

The containment studies indicated that, as the 
lattice size was increased, the sputtering ratio increased. 
Investigation into the similarities of the 17x4x17 and 9x4x9 
microcrystallites, indicated that atom 108 (17x4x17) 
sputtered twice in 36 events. This corresponds to atom 40 
in the 9x4x9 microcrystallite, which did nor sputter in 36 
events. Sputtering mechanism trace runs were performed to 
determine the sputtering mechanism involved in the 
sputrering of atom 108. Atom 108 sputtered when ions 
impacted at points 50 ( impact point displaced 0.50 LU in 
the Z direction only) and 60 ( impact point displaced 0.6 LU 
in the Z direction only) . 

Figure 118 illustrates the atom trajectories in 
the sputtering of atom 108 from impact point 50 ( the 
mechanism is the same for impact point 60) . The copper ion 
struck atoms 74 (PKA) and 91 forcing them down by timestep 
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45. As atom 108 was going upward, it also had a velocity in 
the positive Z direction of greater than 6000 m/sec. Atoms 
116, 117, and 125 bent atom 108 upward quite radically. The 
velocity in the Z direction had been reduced to less than 
700 m/sec by timestep 50. By timestep 60, atom 108 was far 
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enough above the lattice surface, with sufficient velocity 
to be classified as a sput. Atom 108 was ejected at an 
angle of 19.23° to the surface normal with an upward 
velocity of 3206 m/sec. 

Quite evident is the role atoms 116, 117, and 
125 have in channeling atom 108 upward. In the 9x4x9 
lattice, atom 40 ( a boundary atom) corresponds to atom 108. 
And nothing equivalent exists for atoms 116, 117, and 125. 
Investigation was then made to determine whether an atom 
corresponding to atom 103 could be sputtered in the 9x4x9 
lattice by using the same impact point on a different PKA. 
In this simulation, atom 31, corresponding to atom 109 in 
the larger lattice, was expected to sputter. Figure 119 
illustrates the atom trajectories when the copper ion struck 
the new PKA (atom 13) at impact point 50. Atom 31 was 
ejected at an angle of 20.17° with an upward velocity of 
3127 m/sec. The mechanism is the same, and the importance of 
boundary atoms is very evident. 

Additional runs were made on the 9x4x9 lattice 
bombarding atom 13 at impact points 50 and 60. At both 
points, comparison is made between runs with (1) only atoms 
of the natural mass, , and natural mass, and (2) forcing 
atom 22 (91) to a heavy isotope and atom 67 (224) to a light 
isotope. In all four cases, the sputtering mechanism is rhe 
same. The differences in the ejection angles and upward 
velocities, for each impact point, are less than 0.4° and 75 
m/sec respectively . 

3 • Io n/At om Mass Ra tio 

Sigmund [28] proposed a theory that the maximum 
sputtering yield increases as a function of energy. 
However, results of this investigation strongly indicate a 
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dependence on the ratio of the ion mass to the target, atom 
mass. A rearrangement of Sigmund's theoretical curves 
(figure 16b, ref. 28) to show the sputtering ratio as a 
function of the ion/atom mass ratio results in figure 120. 
The curves shown by Sigmund are predicted for 200-eV, 
600-eV, 1-keV, and 5-keV ions bombarding an aluminum target. 
Figure 120 illustrates that Sigmund's theory would predict 
the sputtering ratio to increase as a function of the 
ion/atom mass ratio. Figure 121 shows experimental results 
obtained by Wehner and is presented by Sigmund ( figure 12, 
ref. 28 ) to support his theory. This figure has been 
replotted to show the sputtering ratio as a function of the 
ion/atom mass ratio. The results were obtained by 
bombarding a copper target with 100-eV, 200-eV, and 600-eV 
ions. Experimental results plotted in this manner show a 
decrease in the sputtering ratio with an ion/atom mass ratio 
greater than unity, in disagreement with the theory proposed 
by Sigmund. Results obtained in the present investigation 
completely agree with the shape of experimental data shown 
in figure 121. 



C. RANDOM IMPACT POINTS 



In response to other investigators' suggestion that the 
target impact area be bombarded randomly rather than the 
systematic grid of impact points, the computer program was 
slightly modified. Since additional data were desired for 
mercury events, Hg + was used as the incident ion. The 
KSE-E potential function was used. Earlier results revealed 
that the 100-eV, mercury ion event was contained in an 
11x4x11 microcrystallite. 
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Figure 122 illustrates the (100) orientation impact area 
that was struck at the 144 randomly selected impact points. 
A few gaps are evident, but the randomness is very apparent'. 
Figures 123-128 show a comparison of the results from a 
randomly bombarded impact area and the results from the grid 
of impact points. Since the number of events were not the 
same, the sputtering summaries indicate a sputtering ratio 
for each atom. Figure 123 shows the normalized sputtering 
summary for the randomly selected impact points. Figure 124 
shows the same information using the grid of impact points. 
The four groups of 36 randomly selected impact points had a 
mean sputtering ratio of 0.812 with a variance of 0.0048. 
Additional runs ( not shown ), combined with the four 
previous groups to give a total of 21 groups of 36 randomly 
selected impact points ( 756 shots ) , resulted in a mean 
sputtering ratio of 0.782 with a variance of 0.016. This 
compares very favorably with a sputtering ratio of 0.778 
using the grid of impact points. Figures 125 and 126 show 
the "probability" results for the randomly selected impact 
points and the grid of impact points respectively. It is 
quite evident that they are almost identical. Figure 127 
illustrates the angular dispersion results for randomly 
selected impact points. A definite shift of the lighter 
isotopes toward the surface normal is observed. Figure 128 
illustrates the same information for the grid of impact 
points. The absence of the shift of the light isotopes 
toward the normal is noted. 
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V. 



CONCLUSIONS AND RECOMMENDATIONS 



A. SIMULATION MODEL 



The more sophisticated a simulation model becomes, the 
more uncertain are its results. Some of these uncertainties 
are now beginning to appear in this model. 

The lattice generators can now construct randomly 
positioned isotopes. Symmetrical unfolding of the event 
results is the reverse of the folding required to reduce the 
representative impact area. The model can now display the 
ejection angle of the sputtered isotopes in a histogram for 
visual comparison. It now determines the distance from the 
surface layer center to each ring of sputtered atoms, the 
number of atoms and the number of times each atom sputters 
in each ring, and the probability than that each ring 
contained an atom that was sputtered. Other modifications 
were made in the model to approximate more nearly a 
meaningful physical decription. These include changes in 
the distances from the lattice surface at which calculations 
are truncated, the depth of lattice needed to prevent any 
change in the yield, and a systematic determination of the 
distance above the surface at which an atom is declared zo 
have sputtered. 
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3. CONTAINMENT 



The containment of an event is paramount in a complete 
description of the sputtering process. However, it is less 
significant when studying mechanisms and trends. 

Based on the definition of containment ( no increase in 
yield with an increase in the microcrystallite size ) , the 
100-eV results showed that a 9x4x9, (100) orientation 
microcrystallite contained only the helium and xenon events. 
The 100-eV, mercury event is contained in an 11x4x11 
lattice. To contain a 100-eV copper event, which had the 
largest yield, a 13x4x13 lattice was needed. Sputtering 
mechanism traces revealed the importance of boundary atoms 
in deflecting and channeling atoms in the sputtering 
process. By choosing the PKA to be an atom other than the 
one at the surface layer center, the effect of an increased 
lattice size could be duplicated. 

Investigation with a 17x4x17 microcrystallite revealed 
that Gibson-II potential function simulations are more 
difficult to contain than when using the KSE-B potential 
function. Only the 100-eV, copper event was contained using 
Gibson-II, while using KSS-B the 100-eV and 200-eV, copper 
events were contained. It is quite obvious that to contain 
any higher-energy events would require a microcrystallite so 
large that the computer run time would be prohibitive. 
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c. 



ANGULAR DISPERSION 



Because the low-energy events resulted in so few 
sputtered atoms, the angular dispersion results for any 
particular ion are inconclusive. A trend of light isotope 
enrichment is observed when the results from numerous runs 
are combined. For example, the histograms of the four sets 
of 36 randomly selected impact points show a definite shift 
of the lighter isotopes to smaller ejection angles. In his 
report on the enrichment of lighter isotopes in the normal 
direction, Wehner [24] noted only a 0.6% enrichment in Cu 
(63) . The number of sputtered isotopes was not sufficient 
to detect this small difference. 

The "reflective collision" theory, which Wehner suggests 
to be the cause, was observed to have very little effect on 
the sputtered atoms. Investigation revealed that in most of 
the simulation runs, the same atoms sputtered independent of 
its isotope type or the various positions of the lattice 
isotopes. When the PKA was a light isotope, the atoms 
that did sputter were ejected more nearly normal to the 
surface. Varying the isotope type of the PKA, atoms that 
were known to sputter, and some of the target atoms that 
they collided with resulted in differences in the ejection 
angle normal of less than 0.4°. 

D. RING FORMATION 



do 

the 



Results definitely show that the sputtering mechanisms 
produce a formation of rings of sputtered atoms around 
impact area. The results indicate that the transfer of 
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momentum along the <100> and < 1 0 1 > directions in the surface 
layer is relatively easy. The predominance of sputtering of 
next-nearest neighbors, in the < 1 0 1 > direction is evident. 
For example, at 100-eV, all simulations using an ion of the 
mass of argon or greater sputtered the same atoms, but with 
different frequencies. The particular atoms that sputter do 
not appear to be dependent on the mass of the incident ion. 
However, the number of times the atom sputters does appear 
to be dependent on the ion mass. Ring and atom sputtering 
probabilities clearly show the sputtered atoms from certain 
site directions and the distances of the sites from the 
center of the PKA. 

A trend was noted: as the energy and mass of the ion is 
increased, the sputtered atoms come from sites which appear 
to spread outward from the impact area. These higher-energy 
events give the appearance of a "liquid splash", with the 
size of the crater dependent on the ion energy and mass. 

3. ION/ATOM HASS RATIO 



The results of this investigation do not support a 
monatonic increase in the sputtering ratio with increased 
ion mass as implied in figures of Sigmund[28] and Eer Nisse 
[27]. All the plots showing the sputtering ratio as a 
function of the ion mass, clearly indicate a maximum near an 
ion/atom mass ratio of one. At ion/atom mass ratios greater 
than one, the sputtering ratio appears to become asymtoptic 
to a lower value. This may be an artifact of the computer 
model. These plots also highlight the containment problem. 
The behavior of the sputtering ratio, for larger ion/atora 
mass ratios, is erratic because the microcrystallite does 
not contain the event. 
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The general shape of the sputtering ratio curve, when 
plotted as a function of the ion/atom mass ratio, is in 
agreement with the experimental data of 
co-workers. Results conclusively show 
mass ratio is an important factor in the 
the sputtering ratio. 



Wehner and his 
that the ion/arom 
determination of 



F. RANDOMLY SELECTED IMPACT POINTS 



The sputtering ratio obtained using the grid of impact 
points is extremely close to that obtained when randomly 
selected impact points were used. The similarity between the 
sputtering summaries and the atom sputtering probabilities 
are quite evident. The quantity of data obtained with the 
randomly selected impacr points was required to confirm that 
light isotope enrichment was present. 

Based on the comparison of results obtained from the 
randomly selected impact points and from the grid of impact 
points, the systematic method of "walking" through the grid 
seems completely justified. 



G. FUTURE RESEARCH 



Even the more than 1000 hours of computer time used in 
running these simulations did not begin to encompass the 
permutations and combinations of parameters that need to be 
investigated . 
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The (110) and (111) orientations need to be studied in 
depth. The few sputtering mechanisms observed were the same 
as those for the (100) orientation. However, the sputtering 
ratio as a function of the ion mass, containment 

requirements, ring formation, and angular dispersion of 
isotopes need further investigation. It is highly 

recommended that the effect of momentum transfer along 
certain directions be investigated. 
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APPENDIX A 



ISOTOPE GENERATION 



All atoms in the lattice were one of two types; type 1 - 
the heavier, least abundant isotope, or type 2 - the 
lighter, most abundant isotope. The average atomic weight 
of copper is 63.54 amu. The isotope weights are 62.93 amu 
and 64.93 amu with natural abundances of 69. 1% and 30.9% 
respectively . 

The lattice was generated with either the isotope masses 
and their corresponding abundance or the weighted average of 
the isotope masses with an abundance of 100%. Initially, 
all atoms were forced to the type 2 atom in the lattice 
generator. If the isotope lattice was desired, the 
International Mathematical and Statistical Library (ISML) 
routine GGOB was called. This routine generates a uniformly 
distributed set of psuedo - random numbers between zero and 
one. Each of these random numbers was assigned to a 
particular atom and was tested to see if it was less than or 
equal to 0.309. If the number was less than this test value 
(the abundance of the least abundant isotope) , the atom was 
changed to a type 1 atom. A seed for the random number 
generator, for the particular size lattice to be generated, 
was chosen by trial and error which would produce the 
correct number of atoms above and below the test value to 
assure the correct isotope ratio. 
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A feature was designed to allow any number of the atom 
masses (types) to be selectively changed after the lattice 
had been generated. This ' fearure was incorporated to 
investigate the changes and effects of picking heavy/light, 
heavy/heavy, or light/light collision atoms. This would 
facilitate the dynamic investigation of the ''reflective 
collision" theory proposed by Winters and Sigmund [10]. 

The dynamics of the simulation remained unchanged except 
that all calculations that were dependent on target atom 
mass were modified to use the correct mass. 

The simulation was now capable of generating a lattice 
with the correct isotope ratio and masses with randomly 
selected sites. 
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APPENDIX B 



LATTICE FOLDING 



As mentioned in the simulation model development, a 
lattice that was symmetrical about the folding axes was 
required. Once this was achieved by the lattice generators, 
it was necessary to fold the lattice in a manner such that, 
after the folds, the reduced representative impact area was 
identical to the complete representative impact area. After 
these folds' were performed correctly and in the right 
sequence, the results were identical to that obtained if the 
complete representative impact area had been bombarded 
systematically. 

After each fold, the sputtering count for all atoms that 
were symmetric about the folding axis were summed. The 
summing process started with the first atom in a layer to 
the highest numbered atom in that layer After this process, 
the sputtering count for each atom was twice what the sum 
was since the sputs were summed each time a symmetric atom 
was encountered. This doubling effect of each fold was 
corrected by scaling the sputtering count of each atom after 
all the folds in a layer were completed. All necessary 
folds were then performed on each successive layer. 
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A. (100) ORIENTATION FOLDING 

Figure 12 shows the (100) orientation reduced 
representative impact area folding geometry. The reduced 36 
impact point area, when folded, represents the complete 298 
impact point area. 



S ince 


the 


(100) 


orien tation 


im 


pact 


area 


is a 45°-45° 


triangle , 


three 


folds 


were required 


to 


complete 


the folding 


process . 


The 


first 


fold was along 


the 


diangon The first 



fold was along the diagonal axis, the second along the 
horizontal axis, and the third along the vertical axis. The 
diagonal fold process for all layers was the same. However, 
the horizontal and vertical folds were modified depending on 
which layer was being folded. 

After all the folds in a layer for the (100) orientation 
were complete, the sputtering count for each atom was 
divided by eight to correctly scale the results. 

B. (110) ORIENTATION FOLDING 

Figure 12 also shows the (110) orientation reduced 
representative impact area folding geometry. The reduced 36 
impact point area, when folded, represents the complete 144 
impact point area. 

Since the (110) orientation reduced impact area is a 
rectangle, two folds were required to complete the folding 
process. The first fold was along the horizontal axis and 
the second fold was along the vertical axis. The diagonal 
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fold was not necessary. The horizontal and vertical folds 
(exactly the same as for the (100) orientation) were 
modified depending on which layer was being folded. 

After all the folds in a layer for the (110) orientation 
were complete, the sputtering count for each atom was 
divided by four to correctly scale the results. 



C. (Ill) ORIENTATION FOLDING 



The (111) orientation was by far the most difficult to 
determine. Figure 13 shows the (111) orientation reduced 
representative impact area folding geometry. The reduced 21 
impact point area, when folded, represents the complete 126 
impact point area. 

Since the (111) orientation impact area is a 30°-60° 
triangle three folds along the diagonals resulted in the 
overlap shown in Figure 13. However, if the set of three 
folding sequences were summed the results could be scaled to 
represent the complete impact area. The number one fold was 
along the positive slope diagonal axis, the number two fold 
was along the negative slope diagonal axis, and the number 
three fold was along the vertical axis. 

Before each sequence of folds, the sputtering count for 
each atom was reinitialized to the original value. The 
first sequence (1,2,3) was performed and upon completion the 
sputtering count for each atom was divided by eight to 
correctly scale the results. The second sequence (1,2,1) 
was performed and again the sputtering count for each atom 
was divided by eight for scaling. The third sequence 
(1,3,2) was performed and again the sputtering count for 
each atom was divided by eight for scaling. 
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was 



Finally, the results obtained from each sequence 
summed and the final sputtering count for each atom 
divided by four to correctly scale the results due to 
summing process. 

Since the geometry of the (111) orientation 
complicated and not fully developed, the folding process 
performed only for the first layer. 



was 

the 



was 

was 
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APPENDIX C 



GLOSSARY 



ABON1 , 


ABUN2 




The abundance of the target's isotopes. 


AC 


A distance measurement used in impact. point 




generati on . 



ALFA , ALFA2 





The cosine of the angle between vectors R1,R2, 




squared. 


ALPHA 


Input Morse potential parameter. 


ANGLE 


The angle between the direction of motion of an 




ejected atom and the surface normal in degrees. 


AYE 


The sputtering ratio. 



AXT, AYT, AZT 





Components of the random thermal displacement. 


BENRGY 


The lattice binding energy in ev. 
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BIND 



Negative of the total potential energy (TPOT) at 
time zero. Retained for energy conservation 
checks . 

BHAS 

The mass of the bullet in amu. 

BQZ 

A fee crystal constant - (1/ROOT (3) ) . 

BOL 

The bullet heading Hollarith. 

BX , BZ 

Unsealed x,z coordinates of the impact point. 

BXMAX, BXflIN 

X direction boundaries of the target impact area. 

BZMAX, BZHIN 

Z direction boundaries of the target impact area. 

CELS 

A frictional force multiplier when frictional 
forces are included in the program. 

CGB1 , CGB2 

Morse potential parameters. 

CGD1, CGD2 

Morse potential parameters. 

CGP1 , CGF2 

Morse force function parameters. 
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CFO, CF 1 , CF2 

Force parameters of the cubic fit between lorse 
and Born-flayer functions. 

CO, CORF 

"Phonon" displacement correlation parameter used 
as an option in WARM. 

CORN 

1 - CORF . 

COX, COY, COZ 

Directional cosines of the incoming ion's velocity 
vector, relative to the surface normal. 

COXI, CO II , COZI 

Negative values of COX, COY, COZ. 

CPO, CPI, CP2 , CP3 

Potential parameters of cubic fit between Morse 
and Born-Mayer functions. 

CVBD 

A ratio used to avoid repeated division. 

CVD 

Distance conversion factor - angstroms to meters. 

CVS 

Energy conversion factor - ev to joules. 

CVED 

A ratio used to avoid repeated division, CVE/CVD, 
used to convert exponential potential functions to 
forces . 
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CVM 



Mass conversion factor 



- amu to 



kg . 



CVR 

Conversion factor - lattice units to angstroms. 

DBXM AX , DBZMAX 

Upper boundary limits for the target impact area. 

DCON 

Input Morse potential parameter. 

DFF 

The distance closer than ROE that an atom is to 
the bullet. Used when forces are turned on or off. 

DIST 

The distance between any two atoms. 

DNN (I) 

Distance to the Ith ring in lattice units. 

DBX , DRY, DRZ 

Components of DIST. 

DT 

Length of each timestep in seconds. 

DTI 

Number of lattice units the most energetic atom 
can move in one timestep. 

DTOD 

DT over distance - a ratio used to avoid repeated 
division . 
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DT0M1, DTOH2 

DT over target isotope mass - a ratio used to 
avoid repeated division. 



DTOMB 

DT over bullet mass - a ratio used to avoid 
repeated clivision. 

DX(I) , DY (I) , DZ ( I) 

Change in the position of the Ith atom. 



EHAX 

The velocity squared of the fastest atom. Used to 
make ETEST. 



ETEST 



The kinetic energy of the most energetic atom. 



ETHR 



Energy threshold to escape the surface. 



EV 



Primary energy of the bullet in ev. 



EVR 

Primary energy of the bullet in kev. R indicates 
a read-in and saved value. 



EXA, EXB 

Input Born-Mayer potential function for the target 
atom - target atom interaction. 



FA 

The temporary component force increment on an 
atom . 
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FAC 





The decrement size used in SETUP in locating the 
ion impact point. 


FINISH 


The elapsed program running time in seconds. 


FLAT 


A lattice unit in angstroms. 


FM 


A small number used in checking potential energy 
zero point. 


FM AX 


The maximum frequency of occurrence in the 

histogram intervals. 


FOBCE 


A numerical value of the force function for a 
specific pair of atoms. 


FPTC 


The corrective force at ROE. 


FBC 


Numerical value of the target force function at 
ROE. 


FB AN 


A real variable for number of passes completed at 
termination . 


FBANB 


The reciprocal of FRAN. 
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FREQ1, 


FREQ2 

Number of events in each interval of the 

histogra m . 


FSCAL 


A histogram scale factor. 


FTDI 


A parameter used to determine DT by maximum energy 
method . 


FIA 


4 

A Born-Mayer force function parameter. 


?X(I) , 


FY(I) , FZ ( I) 

Components of the total force on the Ith atom. 


HBHAS 


Half Bullet Mass - a ratio used to avoid repeated 
division . 


HDTOD 


Half DT over distance - a ratio used to avoid 
repeated division. 


HDTOM 


Half DT over most abundant isotope mass - a ratio 
used to avoid repeated division. 


HDTOMB 


Half DT over bullet mass - a ratio used to avoid 
repeated division. 



HTHA$1, HTMAS2 

Half target isotope mass - a ratio use 1 to avoid 
repeated division. 
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IB0X1, 


IBOX2 

Ejected atom angle tally keeper. 


ICOUNT 


Flag used in sequencing. 


IDEEP 


The first frozen layer. in the target 

microcrystal . 


IET 


Elapsed time in 2.6 microsecond blocks. Used to 
determine actual elapsed time. 


IFFM 


The atom number of any target atom that has its 
mass constrained due its isotope ratio. 


IFLAG 


Flag used in sequencing. 


IFB 


The isotope type (heavy or light) of any target 

atom that has its mass constrained due to its 
isotope ratio. 


IFT 


A number used as a constant when generating 
pseudo-random numbers. 


IH 1 


Alphanumeric array of the program heading. 


IH2 


Alphanumeric array of the Morse parameters. 
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IHB(1) , 



IHB(1) , 


IHB (2) , IHB (3) , 136(4), IHB (5) , IHB (6) 

Alphanumeric array describing the bullet. 


IHBPOT 


The type of ion-atom potential. 


IHBOLL 


The bullet elemental symbol or name. 


IHS 


Alphanumeric array for the type and orientation of 




the crystal. 


IHT 


Alphanumeric array describing the target. 


IHT ( 1) , 


IHT (2) , IHT (3), IHT (4) , IHT (5), IHT (6) 




Alpha numeric array describing the target. 


IHTARG 


The target elemental symbol or name. 


IHTPOT 


The type of target atom - target atom potential. 


IL ATT 


A flag used in the latrice generators. 


ILAYSR 


A variable used to identify which layer is being 




folded . 


ILL 


The number of the last atom in the first layer. 
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IMF 


The mass of each target atom with the correct 




isotope ratio. 


IN 


The number of intervals in the histograms. 


INDEX 


A variable used in determining phase in a dynamics 




cycle step. 


INT 


The width of the intervals in the histograms. 


IOR 


The target plane of incidence. 


IOUT 


Flag used in sequencing. 


IPASS 


The number of passes (shots) completed. 


IPNDM 


The impact point number. 


IPNUMS 


Scaled impact point number. 


IPROG 


Input variable used to skip attractive potential. 


IRAN 


Random number seed used in NORMAL subroutine which 




generates pseudo-random numbers. 
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IRANGE 





The upper limit in the histograms. 


IREAD 


The number of atoms which have the isotope type 
constrained externally. 


ISEEDA , 


ISEEDB, ISEEDC 




Seeds used in the uniform random number generator 
which determines the isotope ratio and 
distribution . 


ISET 

ISHOT 


Flag used in sequencing. 



A parameter used to shut down the program. 



ISNUH 


The atom number of a sputtered atom. 


ISOGEN 


A flag used as an option to generate the target 




atom masses with isotope or single element values. 



ISPX, ISPZ 





Initial start point x,z coordinates. 


IT 


Unsealed x coordinate used in lattice generation. 


ITEHP 


The lattice temperature in degrees Kelvin. 


ITIME 


Flag used in sequencing. 
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ITT 



An odd-even integer used to determine atom site 
location . 



IX, II, IZ 

The number of x,y,z, planes in the crystal. 

IXP, IIP, IZP 





Crystal dimensions in x,y,z. 


JN 


The numerical value of the histogram intervals. 


JOUT 


A histogram sequencing variable. 


JSCAL 

JT 


A histogram scale factor. 



Unsealed y coordinate used in lattice generation. 



JTS 


Variable used to establish atom site. 


JTT 


Variable used to establish atom site. 


K2S0M 


The square of the sputtering count which is used 




in calculating the variance. 


KATOM 


The number of atoms in a ring. 
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KFIRST 


A sequencing variable. 


KRING 


The highest numbered ring in a layer. 


KSPUT 


The number of sputs from all atoms in a ring. 


KSUM 


The total sputtering count before folding. 


KT 


Unsealed z coordinate used in lattice generation. 


LCUT (I) 


Used to identify an I th atom which can be omitted 




from the force and potential calculations by 




changing LCUT(I) from zero to one. 


LD 


The highest numbered atom in a mobile layer. 


LL 


The highest numbered atom in the entire crystal. 


LL2, LL3 


Twice and triple the value of LL . 


LNUM (I) 


The number of times the Ith atom sputters after 




all shots in the target impact area. 


LNUMF (I) 


The number of times the Ith atom sputters after 




all shots and the summary is complete and folded. 
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LNUMI (I) 





The number of times the Ith atom sputters after 
all shots but before folding the summary. 


LS 


Sum of the Miller index integers. 


LSS 





Variable used to identify the type of surface. 



LSDR 


The lattice surface type. 


LTYPE 


The type of lattice (i.e. fee ) . 


MAX 


A height scale factor used in the histograms. 


MCBO 


The number of simultaneous equations solved by 




CROSYM. 


HI 


NLAY - 1, a sorting variable. 


NCT 

NCTT , ND 


Increments the timestep interval. 



Data output interval in number of timesteps. Also 
the timestep interval between potential energy 
calculat ions . 



NF 1 


A sequencing variable. 


NFIRST 


The number of the first atom in a given layer. 


NLAY 


The highest numbered atom in a given layer. 


NLAY1 


The highest numbered atom in the first layer. 


NLINE 


The line number in the output block. 


NOTH 


A Hollarith consisting of blanks. 


NPAGE 


A page numbering variable. 


NRAN 


Variable for number of passes (shots) completed at 




termination . 


NBS 


A restart flag. 


NS 


The next timestep number at which output will be 




printed. 
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NSHOTS 





The number of passes (shots) to be completed in a 
program. Used to break up an extremely large, 
many-shot program into a few small programs with 
manageable computer run time. 


NSPUT(I) 


The number of atoms that sputter for the Ith 
impact point . 


NSR 


The first timestep where output is printed. 


NSS 


Flag to print final output status of current shot. 


NT 


The timestep number. 


NTT 


Timestep number limit before shut down. 


PAC 


Parameter for bullet force function correction. 


PBMAS 


Bullet mass in kg. 



PBXA , PBXB 



PFIV 


Input Born-Mayer potential function parameters for 
the bullet-target interaction. 

Point five (0.5). 
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PFPTC 





Bullet force function evaluated at ROE. 


PFRC 


Bullet-atom repulsive force interaction evaluated 
at ROE. 


pfxa 


Bullet force function parameter. 


PKE (I) 


The kinetic energy of the Ith atom. 


PNEG 


Total negative potential energy in the lattice. 


PNOM 


i 

Alphanumeric variable for the impact point number. 


POT 


The potential energy between two atoms. 


PPE(I) 


The potential energy of the Ith atom. 


PPOS 


Total positive potential energy in the lattice. 


PPTC 


Bullet potential function evaluated at ROE. 


PfiOB 


The probability that an atom will sputter from a 
given ring. 
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PTC 


Target atom potential function evaluated at ROE. 


PTE (I) 


The total energy of the Ith atom. 


QH 


A small number used in checking the kinetic energy 
zero point. 


QUIT 


Coutoff variable checked against the total 
potential energy which can be used as a shut down 
cr iter ion . 


HI, R IS 


Distance from (0,0) impact point to the initial 
bullet position. 


R1, R1T, 


R 1 Z 

The x,y,z components of R1S. 


R2 


Magnitude of vector from impact point to the first 
atom hit by the bullet. 


R3 


Magnitude of vector from impact point to bullet 
start position. 


R 1 , R2 


Scalar product of vectors R1, R2. 


R2SQ 


R2 squared. 
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RAD 





The angle between the direction of an ejected atom 
and the surface normal in radians. 


RADDEG 





Conversion factor - radians to degrees. 



RAN (I) 


A list of normally distributed random numbers. 


RANGE 


A scaling variable used in the histogram. 


RBI, RBZ 


Unsealed x,z coordinates of the impact area 




reference point. 


RDIST 


The distance between an atom and the center point 




of its layer. 


RE 

RINCR 


Input Horse potential parameter. 



The RANGE increments. 



RMIDPT 


The midpoint of each interval in the histogram. 



ROE, ROE2 





Nearest neighbor distance - squared. 


ROEA 


The maximum cutoff for the Born-Mayer potential. 
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EOEB 

The minimum cutoff for the Morse potential. 

ROEC, ROEC2 

The maximum cutoff for the Morse potential 
squared. 



ROEM 

Defines the region where the repulsive force 
function must be modified. 



HI (I) 



RY(I), RZ(I) 

X,y,z coordinates 



of the Ith atom at any time. 



RXBL , R YBL, RZBL 

Maximum negative x,y,z coordinates an atom may 
move before it is LCDT. 



RXBH, RYBM, RZBM 

Maximum positive x,y,z coordinates an atom may 
move before it is LCUT. 



RXF, RYF, RZF 

Final x,y,z coordinates of the Ith atom when it is 
LCDT. 

RXI (I) , RYI (I) , RZI (I) 

The Ith atom's initial x,y,z coordinates. 

RXK (I) , RYK (I) , RZK (I) 

Temporary x,y,z coordinates of the Ith atom during 
the force cycle. 

RIS , RYS, RZS 

Scaled x,y,z coordinates of the impact point. 
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SAVE 



One-half POT. 



SBX, SB Z 

Scaled x ,z coordinates of the impact point. 



SCALE 

A scaling variable for the frequency of events in 
the histogram. 

SCX, SCI, scz 

The x,y,z lattice scale factors used to convert, 
the generated lattice to match a real crystal. 

SCXR, SCZR 

The inverse of SCX and SCZ - a ratio used to avoid 
repeated division. 

SDIST 

The square of the distance from the center of a 
layer to an atom. 

SETIME 

A library subroutine used to keep track of run 
time . 

SLOW 

A cutoff velocity variable checked against a long 
DT . 

SPDX , SPDZ 

The fraction of the target impact area that the 
impact points are incremented. 
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SPX, SPZ 





The x,z components of the distance from the impact 
area reference point to the impact point. 


sscz 


A z scale factor used in the (111) orientation 
lattice generation. 


STD 


The standard deviation of the sputtering ratio. 


STOPER 





This variable stops calculating and saves 60 





seconds for compilation and output. 


SOR 


Alphanumeric variable describing the surface. 


TA, TAHP 


The thermal amplitude of the lattice. 


TAR 


Alphanumeric variable describing the target 

element . 


TARGET 


Alphanumeric variable describing the target 

element . 


TAS 


A scale factor - T AbP/ROOT (3.0) . 



TASF, TASN 

Random thermal amplitude scale constants. 
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TE 

The total kinetic and potential energy of all the 
atoms. 

TEMP 

The temperature of the lattice in degrees Kelvin. 

TESTE 

A termination test of energy. 

TFAC 

A time factor ratio used to determine DT by 
maximum force methods. 

THERM 

Mean thermal energy of a lattice atom. Used in 
shut-down tests. 

TIME 

The elapsed calculating time in seconds. 

TITLE1 , TITLE2 

The titles for the histogram outputs. 

TMAS1, TMAS2 

The target isotope masses in amu. 

TPKE 

The total kinetic energy of all the atoms. 

TPOT 

The total potential energy of all the atoms. 

TRON 

The time the program is to run in minutes. 
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V EL 

The velocity calculation - SQRT (ev/HBMAS) . 

vss 

A temporary storage variable for the velocity 
componen ts. 

VTEST 

A test velocity used to determine if the atom has 
sufficient velocity (energy) to escape the 
surface. 

VTS 

A scale factor - 0.8*VTEST. 

VX(I) , VY(I) , VZ(I) 

The x,y,z components of the Ith atom's velocity. 

US, VYS , VZS 

The x,y,z components of the bullet's initial 
velocity . 



XDIST 

The x component of the distance from the center of 
a layer to the atom. 

XMAX, X MIN 

The x planes to determine whether two atoms are 
close enough to interact. 



YLAX (I) 

The crystal relaxation of the Ith layer in the -v 
direction. 
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YMAX, THIN 

The y planes to determine whether two atoms are 
close enough to interact. 

ZDIST 

The z component of the distance from the center of 
a layer to the atom. 

ZE 

Absolute zero (0.0). 

ZBAX, ZBIN 

The z planes to determine whether two atoms are 
close enough to interact. 

ZSTOP 

Used to print out a warning header that a pass was 
terminated before completion. 
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APPENDIX D 



COMPUTER PROGRAM 



GLD3 IS 0LD1 MCOIFIEC TG INCREASE SPEED 
ROLD IS GLD3 WITH 'WARMED' LATTICE 
RING IS RPF MODIFIED TO FOLD LATTICE AND 
STUDY RING FORMATI ON /SPUTTER ING 
RMAS IS RING MODIFIED TO STUDY TARGET ISOTOPES 
RANG IS RMAS MOD. TG STUCY ANGULAR CISP5RSICN 



DIMENSION DX(60 0) ,DY(600) ,DZ(60C ) 

DIMENSION RXK (600 ) ,RYK(600) ,RZK(600) 

CIMENSICN PKE( 600) ,PTE (600) 

DIMENSION NSPUT (50) , I PNUM S ( 50 ) 

CIMENSICN I SNUM (20) , I S EQ ( 60 0 ) 

CIMENSICN RXF(600),RYF(600),RZF(6C0) 

CIMENSICN 180X1(20), IBOX2(20), RANNUM(250) 

CIMENSICN FF6Q1 (80) , F REQ2 ( 80 ) , RANGE ( 20 ) , TIT LEI (24) , 

X T ITL E2 ( 24 ) 

CIMENSICN HCTOM ( 2 ) 

CIMENSICN IFFM(600) ,IFM(600 ) 

CCMMON/CCMl/RX(oOC) , RY ( 600 ) ,RZ ( 60C ) ,LCUT(60C) , LL , LD 
CCMMCN/CCM2/RCE, RCE2, ROEM, AC,PAC ,PPTC,PTC, PFPTC ,FPTC, 

X FM,PFIV,TPCT 

CCMMCN/C0M3/EXA , SXB , FXA , PEX A , P EX 8 , P FXA , T PER M 
CCMMCN/C0M4/IX,I Y ,IZ, IXP,IYP,SCX,SCY,SCZ,IDEEP, 

X REX, RE Z 

C CM MCN/ COM4 A/ LNUM(6 0O) ,LNUMI(600) ,LNUMF(60C) 
CCMM.CN/C0M5/FX( 600) ,FY(600) ,FZ ( 600 ) , P PE ( 600 ) 
CCMMCN/C0M6/C0XI , COY I , CCZ I , RXS ,RYS ,RZS, FAC 
CCMMCN/CCM7/R1, LSS, SPX , SPZ , COX , CC Y ,COZ , Y L AX (20) , I LL 
CCMMCN/C0M8/VX(60C) ,VY(600) ,VZ(60C) 

CCMMCN/C0M8A/VXS ( 600 ) , V YS ( 600 ) , VZS (600) 

CCMMCN/COML 1/RXI <60C),RYI(600) ,PZI (600) 

CCMMCN/ C0M13/ RAN ( 1800) , IRAN, TASK, TASF ,LL2 , LL3 
CCMMCN/C0M14/ AXT (6 00) ,AYT( 600) , AZT(6O0) 

COMMON/ COM 16/BMAS ,TMAS1 ,TMAS2, A8UN1,ABUN2, I FNUM 
CCMMCN/ COM 16A/3X, BY,BZ,EDI SP,CVR 

CCMMCN/C0M17/TIME,NT,DT,DTI , TEMP, CELS, EV,NP AG E,NLINE 
COM MON/ COM 18/ I TE M P ,T AMP , CO RF 

CCMMCN/COM20/IHK20) , IHS( 10 ) , I HT ( 6 ) , IHB ( 6 ) , IHTARG ( 2 ) , 

X IHBUL L ( 2 ) 

CCMMCN/ COM2 1/IHBPCT (2 ) , IHT P CT ( 2 ) , SUR ( 2 ) , I F2 ( 10 ) , I CR , 

X L S UR 

CCMMCN/C0M22/R0EA, RCEB,RCEC,R0£C2,CP0,CP1,CF2,CP3,CF0, 
X CF1,CF2,CGC1,CGD2,CGB1,CG3 2,CGF1 ,CGF2, IPRCC- 
C CMMCN/C0M23/ SPDX , SPDZ ,DBX MAX , DBZMAX , SSCZ 
CCMMCN/ COM2 4/ DNN (100),KATOM(100),KSPUT(100), FROB ( 100 ) 
CCMMGN/ C0M24A/ KSUM,KRING,RDIST(6CC) 

CCMMCN/ COM2 5/ IMF ( 6 00) , ISOGEN 
CCMMCN/C0M3C/ ISEEDA, ISEEDB, ISEEDC 
CCMMCN/ COMA / A(4,5),MCRC 
9005 FORMAT (20A4) 

90 1C FORMAT ( 20A4 ) 

902C FORM AT ( 10 A4 , 1 3 , 2A4 , 1 2 ,3 1 3 , 2F3 . 0 , 2 F5 . 2 , I 2 ) 

903C FORMAT ( 2 A4 , A4 , F 8 . 3 , 2 A4 , 2 F 1 0 . 5 , F 8 . 5 , 3F7 . 4 , I 2 ) 

9025 FORMAT (F5.2,I5) 
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9040 FCRMAT ( 10A4 , 3F5 . 2 , 3F 8 . 5 , 1 1 ) 

905C FCRMAT (2X , 13, 15 ( 13, IX, 1 1, 2XM 

90 SG FCRMAT(5X , 3F5 . 3 ,4 I 5 ,2 F8 .6 , 2X , A4 , 3 F6 .4 ) 

9 1 C 5 FCRMAT ( 1 IX , F5. 0 ) 

96CC FCRMAT ( // ) 

9 6 C 1 FCF MAT ( 1H1) 

96 2C FORMAT (IX, 3 (15, 4X, II, 4X,3F 6.2,8 X) ) 

9630 FCRMAT(105X,4HPAGE,I3,/,1H1 ) 

964C FORMAT ( ' BINDING ENERGY =',F6.2,', TEST 



SPEED = • 



F 10 .2 » ' M/SEC ’ , ', SAFETY FACTOR ' SPEED = • ,F1C.2, 



a M/SEC' , / , 1H1 ) 

96 80 FCRMAT < 2X, • ATGM • , 2X, • LCUT* , 9X , • V X * , 9X , • V Y » , 9X , • V Z • , 
X 8X,'RXF',8X,'RYF',8X,'RZF',8X, 'PPE',8X,'PKE',8X, 

X 'PTE'/) 

9681 FCRMATt IX, • ATOM DX • , 6X , • D Y • , 6X , • D Z • , 7X 



• VY 




X F 10 .3 , ' TOTAL NEGATIVE POTENTIAL 
9695 FORMAT (3X, • ATOM £ POTENTIAL SURF* 
X'VZ',10X,'PKE',6X,'PPE',6X,'PTE 



IRAN = ' ,110 



ENERGY IS ' , F10 . 3 ) 
,5X , 'VX * , 7 > , 'VY* ,7X, 
KEV OCX COY' , // ) 



r • , 



LL 



- 1 



15,//) 



9701 FCRMAT ( ' IFT =',110,', - ,xw, , __ 

9800 FORMAT ( IX, 215 ,2A4, I4,A4,1X,3F9.1, 3F9.4, F5.1 , 1 X , 2F 6 • 3 ) 
98C5 FCRMAT( 5X ,' SPUTTERED ATOM / ANGLE WITH NORMAL'/) 

9810 F CRM AT ( // , ' AT',F7.1,' E V, THE FOLLOWING ATOMS WERE', 

X • SPUTTEREC: • / , 2014, // ) 

9815 FCRMATt 1X,5(5X, I3,3X,F7.4) ) 

982C FCRMAT ( • THERE ARE NO SPUTS AT »,F8.1,'EV') 

98 AC FORMAT ( 25X , 'SUMMARY OF TARGET IMPACT AREA'/,25X, 

X 'FOLDED TO REPRESENT TOTAL COVERAGE'/) 

9845 FCRMAT (IX, 'ATOM NUMBER / TIMES SPUTTEREC •) 

9850 FCRMAT ( ///, ' SPUTTERING RATIO WAS',F6.3, 

X ', THE VARIANCE WAS',F8.3,', FOR ',14, 

X • IMPACT POINTS, THE STD DEV CF THE MEAN WAS', 

X F6.3,//) 

9860 FCPMAT(2X, • STATISTICAL SUMMARY FCR ' , 14 , ' I M P AC T POINTS' 
X ,' AT ' , F 5 . 1 , 'KEV' ,///) 

9865 FORMAT ( ' NUMBER OF TIMES EACH ATCM SPUTTEREC'/) 

987 G FCRMAT (10(16,14) ) 

9875 FCRMAT( / / , ' BINDING ENERGY WAS',F6.2,- cv, ... . ___ 

'■ ‘ ‘ SPUTTERED ATOM COUNT IS', 15) 



• WAS ' , F10.2 , ' M/SEC, 



• EV, TEST S PEEC ' 



9880 FCRMATt IX, 3 (• 

9885 FCRMATt IX, 3 ( • 

X F6 . 3 , 1 X ) ) 

989C FCRMAT( 25X , 'PROBABILITY 
X // ) 

9895 FCRMATt 25X 
99 CC FCRMAT ( 2X , 

X / ) 

9905 FCRMAT ( 2X ,' ANGLE 
9910 FCRMAT “ “ 

9915 FCRMAT 

9920 FCRMAT ( ' 00-05 05-10 



RING DISTANCE ATOMS SPUTS PRCB 
NN( ' ,12,' )' »2X,F6.3,3X,I2,5X, 13, 3X, 



) ) 



•SPUTTERED 

INCREMENTS 



THAT A RING OF ATCMS SPUTTERS' 

ATOMS ANGULAR DENSITY'/) 

ARE I N ' , F4 . 2 , • CEGREE INTERVALS' 



(DEGREES) / NUMBER OF ATOMS'/) 

( 2 X , 19 ( 13, 3X), 13) 

(3X , 12 , 2 ( IX , F4. 2 ) ) 



X • 35-40 
X ' 75-80 
9925 FCRMAT ( ’ 

X • 70-80 
993C FORMAT ( ' .. 

X • CHANGEC 
9935 FCRMAT (IX 



45-50 



10-15 

50-55 



15-20 

55-60 



994 u 
9945 
995C 
9955 



9960 



40-45 
80-85 85-90* /) 

00-10 10-20 20-30 30-40 
80-90'/) 

NUMBER OF ATOMS CHANGEC 
AND ISOTOPE TYPE ARE' 
13(13, 2X ± I1,3X)) 



20-25 

60-65 



25-30 

65-7C 



30-3 5 ' 
70-75 ' 



40-50 5C-6C 60-70' 
WAS', 14/,' ATOMS', 



t 1J ( 1 J? 4A J i I ? Jm i , 

FCRMAT (IX , 'SCX = ' , F7 .4 , • SCY=',F7.4,' 

FCRMAT ( 3 ( I10,10X) ,15) 

FORMAT ( ' ISEED = ' ,110) 

FCRMAT ( IX, •*** THIS SHOT WAS TERMINATED 
X • WAS COMPLETED'/, 4X, ' ZST0P=1 THE CURRENT* 
X • TIMESTEP EXCEEDED THE LIM IT • / , 14X , • ANY 
X • ARE REMOVED FROM THE SUMMARY') 

FCRMAT ( IX ,' *** THIS SHOT WAS TERMINATED 



SCZ=', F7.4) 



BEFORE IT 



SPUTS ' 



3EFCRE IT* 
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X 

9965 F 0 F MA T ( 4X , 1 IMPACT PCINT ',15) 

99 70 FORMAT < 4X , ' BX=« ,F6.4, ■ B Z= • , F6.4) 



x - ■ 




77ivj runi-Mi \ - da- - , r o . n » , r o • H J 

9975 FCFMAT ( 8X , ' NOT IN TARGET APEA',/) 

CALL CV F LOW 
CALL SETIME 
READ (5,9005) TITLE1 
READ ( 5 , 90 C 5 ) TITLE2 
REAC(5,91C5) TRUN 
S TO PE R=TRUN*6 0.0-60.0 
FINISH=0.0 
I LATT =0 
NCFC=4 
NCTT = 10 
TESTE=2 . 0 
El ND=0. 0 
CVB=0.774E-9 
E ENR C-Y = 2 .0 
I ShLT=0 
L C =0 

FM=1 .0E-10 
P F I V=0 . 5 
C M =C .0 1 
SLCW=1. OE-14 
WPITE(6 ,9601) 

CC 1 1=1,20 
I ECX1 ( I ) =0 
I ECX2 ( I ) = 0 

1 CONTINUE 

DC 2 1=1, 60C 
LNUM( I ) =0 
LCLT ( I ) =0 
RXK( I )=C.O 
R YK ( I ) =0 - 0 
R ZK ( I )=C.O 
VX (I )=0 .0 
VY (I) =0.0 
V Z ( I ) =0 . 0 
R X ( I ) =0 .0 
RY( I ) = J .0 
R Z ( I ) =0 .0 
FX ( I )=0 .0 
F Y ( I )=0.0 
F Z ( I ) = 0.0 
PTE ( I ) =0.0 
PKE(I)=C.O 
PPE ( I )=0.0 
R X I { I ) = 0 . 0 
RYKI ) = 0.0 

2 R Z I ( I ) = C.O 

READ ( 5 , 9 CIO ) I H 1 

READ( 5, 9020 ) IHS , IOR,SUR,LSS ,IX,IY,IZ,RBX,RBZ,CTI,RCE2, 
X LTYPE 

READ (5,9030) IHBULL, BUL, BMAS , IFBPQT ,PEX A , P EXE, XXXX 
REAC( 5, 9030 ) IHTARG, TAR , TMAS 1, IHTPOT ,EXA , EXB ,FLAT, TMAS2 
X ,ABUN1 ,ABUN2, ISOGEN 

READ ( 5,9040) IH2 »ROEA,ROEB,ROEC , DC CN, ALPHA , RE , IPRCG 
5 REAC(5,9080)EVR,SPX»SPZ»NTT,NSR,NC»IDEEP , COX, COY, I TEMP 
X ,7A ,CG ,5THR 

RE AC (5, 994 5 ) I SEE DA , I SE ED 9 , I S EE CC , NS HOT S 

READ (5,9915) IR ANC-E , RM IDPT , R INCR 

READ (5,9050 I RE AD , ( I FFM ( J ) , I FM ( J ) , J = 1 , I RE AD ) 

IF(EVR.LE.O.O) GO TO 9999 

IF ( ISOGEN. NE.O) GC TO 33 

TMAS2=63. 54 

A EUN1 = C . 000 

AeLN2=l .000 

33 CONTINUE 
T AMP = TA 
CCRF=CC 
I PASS = 1 
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c 

c 

c 



c 

c 

c 



c 

c 

c 



c 

c 

c 



NCTT=ND 

IF( ILATT.NE .0) GO TO 10 

MAKE HEADING HOLLAR ITH 

IHE (1 )= IHBULLI 1 ) 

IHB(2) = IHBULLI 2) 

I HE (3 ) = IHT A RG ( 1 ) 

IHe(4) = IHTARG(2 ) 

I HE (5 ) = IHB PCT( 1 ) 

I He (6 ) = IHBPCT ( 2 ) 

IHT ( 1) = IHB ( 3) 

IHT (2 ) = IHB ( 4 ) 

IHT(3) =IH8 (3) 

IHT(4)=IHB(4) 

IHT (5) = IHTPCT(l) 

I HT ( 6 ) = IHT PCT ( 2 ) 

MAKE CONSTANTS 

RACCEG=57. 2957795 
P C E = SQRT < RC E2 ) 

EV=EVR* 1C0C .0 
NS=NSR 

FAC=0. 1 
PCEM=RC E-CT I 
CVR=0 . 5 *FLAT 
C V E = 1 . 6 02 IE-19 
CVM=1 .6604E-27 
CVC=CVR*1.0E-10 
C VEC=CVE/CVD 
RXBL=-RCE 
RVBL=-J .8 
RZBL=-RCE 
IXP = ( IX + 1 ) / 2 
I V F = ( IY + LJ /2 
1 2 P = ( IZ + 1 ) /2 
PEMAS=BMAS*CVM 
HeMAS=0.5*PBMAS/CV5 
PTMAS 1 = TMA S 1*CVM 
PTMAS2=TMAS2*CVM 
HTMAS 1= 0. 5*PTMA$1/CVE 
HTMAS2= 0. 5*PTMAS2/CVE 
VTEST=-SQRT (2.0*ETHR#CVE/PTMAS2 ) 
VTS=0.8*VTEST 

MAKE REPULSIVE POTENTIALS 

FXA=ALQG(-EXB*CVED )+EXA 
PFXA=ALCG(-PEXB*CVED)+PEXA 
AC = ALOG (CVEO+EXA 
P4C=AC-EXA+FEXA 
P TC = EXP ( EXA + EXB^ROE ) 

FRC=EXP (FXA+EXB^RCE) 
FPTC=EXF(AC+EXB*RCE) 

PFPTC = EXP ( PAC + PEXB*ROE ) 
PFRC=EXP(PFXA+PEXE*ROE) 
PPTC=EXP(PEXA+PEXE*ROE ) 

GENERATE THE LATTICE 

IF (IOR.EQ. 100) GC TO 6 
IF (IOR.EQ. 110) GC TO 7 
IF (IOR.EQ. Ill) GO TO 8 

6 CALL KF 100 
GC TO 9 

7 C£LL KF 110 
GC TO 9 

E CALL KF 1 1 1 
GC TO 9 
5 CONTINUE 

IF ( IOR. EQ.100 ) THETA = 1 . 0000 
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I F ( ICR. EQ. 110) TH6TA = C.OOOO 

IF( IOR. EQ.lll ) THETA= ( SQRT ( 1.5)/3.0)/($QRT(2. 01/2.0) 
I F ( ISOGEN. EC. 0) GC TO 12 
DC 12 1 = 1*1 READ 
J = I FFM ( I ) 

I N F ( J) = IFM( I) 

12 CONTINUE 

CC 65 1=2, LL, 3 
K = I + 1 
J = I +2 

65 WRITE <6, 9620 1 1, IMF(I),RXI(I),RYI(I),RZI(I),K,IMF(K), 
X RXI(K) »RYI(K),RZI(K),J,IMF(J) »RXI(J)»RYI(J),RZI(J) 
WRITE ( 6, 96C0) 

WRITE ( 6,9933) IREAD 

WRITE ( 6,9925) ( I FFM ( J ) , I FM ( J ) , J = 1 , 1 READ ) 

WRITE (6,9600) 

WRITE (6,9940) SCX,SCY,SCZ 

L 12=2*1. L 

LL3=3*LL 

R 1S = R1 

I LATT=1 

NFAGE=1 

BX=R8X*SCX 

ez=Rez*scz 

sex=Bx 

sez=bz 

RXEM= (IX-1 )*SCX+RCE 
R Y BM = ( I Y— 1 ) *SCY + RCE 
RZEM=(IZ-l)*SCZ+RCE 
10 CCNTINUE 

CCZ=1 . 0-CCX*CCX-CCY*COY 
C0Z=A8S (COZ ) 

CC Z = SQR T ( CC Z ) 

CCXI=-CCX 
CCYI=-CCY 
CCZI=-CCZ 
I FT = IT I ME ( I I) 

IRAN=2* IFT+1 

I F ( IRAN .1 T . 0 ) I RAN=-I RAN 
EXNIN=8X 
BXMAX=EX+DBXMAX 
eZN IN = B Z 
EZMAX=EZ+D8ZMAX 
I F ( IOR. ME. 100) GO TO 25 
IF(SPX.EQ.J.O) GO TO 25 
EZMAX=8ZMAX-10.0*SPX*SPCZ 
e ZMIN=BZMIN+10.0*$PX*SPCZ 
25 I F ( ICR. NE • 1 1 1 ) GO TO 30 
IFISPX.EQ.O .0) GO TO 30 
EZMAX=BZMAX-5.0*SFX*SPCZ 
E ZM I N= BZM I N+5.0*SFX*SPCZ 
30 CCNTINUE 

BX=(RBX+SPX )*SCX 
BY=(RYI ( 4 i + C .00 1 ) *SCY 
e Z=(RBZ+SPZ )*SCZ 
S CXR= 1.0/SCX 
SCZR = 1 . C/SC Z 

I F ( IOR. EQ.lll) SC ZR= 1 • O/SSC Z 
C 

C PREPARE TC WARM LATTICE 

C 

I FT = I T I ME ( I I) 

I RAN= 2* IFT + 1 

I F( IRAN .LT. 0) I RAN=- IRAN 
WRITS (6,9600) 

WRITE(6 ,9701) I FT , IRAN ,LL 
WR ITE (6 ,9640) E TH R , VTE ST , VT S 
IF(TAMP.EQ.C.O) GC TO 40 
TAS=TAMP/SQRT (3.0 ) 

CCRN=1 . C-CCPF 
T ASN = TA S*CC RN 
TASF=TAS*CCRF 
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EACH REPLICA STARTS HERE 

40 CCNTINUE 

IF( IPASS.GT .N SHOTS) GO TO 1000 
I F (TAMP .GT . C.O ) CALL WARM 
50 CCNTINUE 
T I ME = 0 . C 
NS=NSR 
NT=C 
N SS = 0 
NCT = 0 
INCEX=0 
ZSTCP=0 .0 
S PX = ( BX-S8X ) *SCXR 
SFZ=( EZ-SBZ )*SC ZR 
ISPX=( SPX + 0 .001 ) *100. 0 
ISPZ=(SPZ+0.001)*100.0 
I FNUM1C0* ISPX+I SPZ 
IFNUMS ( IPASS) = IPNUM 
R 1=R1S 

RECONSTRUCT LATTICE 

I SEC ( 1 ) =0 
CC 60 1=2, LL 
LCLT ( I ) =0 
I S EC ( I ) =0 
RX ( I) =RXI ( I ) 

R Y ( I ) =R Y I ( I ) 

R Z ( I ) =R Z I ( I ) 

P T E ( I ) = 0 . 0 
PPE ( I ) = C .0 
VX(I) = 0.0 
V Y ( I ) = 0.0 

V 2 ( I ) = 0.0 
60 CCNTINUE 

IF (TAMP .ME .C.O) CALL WARM 
R XS = EX 
R Y S = BY 
R ZS = EZ 
R 1 = R 1 + F AC 
CALL SETUP 
LCL T ( 1 ) =0 
RX I ( 1 ) =RX ( 1 ) 

RYI ( 1 ) = RY ( 1 ) 

RZI (1 ) = RZ(1 ) 

P F E ( 1 ) = 0 . 0 

VEL = SQRT ( EV/HBMAS ) 

VX ( 1 ) =V EL*CCX 
VY ( 1 ) =V EL*C CY 

V Z ( 1 ) = V EL *C CZ 
F CT I =DT I *C V C 
CT=FDTI/VEL 

BEGIN TIMESTEP CYCLE 

ICC CTCD=DT/CVD 

F CTCD = 0 .5 *QTGD 

CTCM1=DT/PTMAS1 

CTCM2=CT/PTMAS2 

HCTCM( 1 ) =0 . 5*DT0M 1 

FCTCM(2 )=0.5*DT0M2 

CTCMB=DT/P6FAS 

FCTGMB=C.5*DT0MB 

2 CO CALL STEP 

I F ( INDEX. GT .0 ) GG TO 250 

FIRST HALF OF TIMESTEP 

2 1C INDEX = 1 \ 

1 = 1 
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I F (LCUT ( I ) . EQ.2 ) GC TO 215 
IF(LCUT(I) • N E . 0 ) GG TO 220 
R X K ( I ) =RX ( I ) 

RYK (I ) = RY ( I ) 

R ZK ( I ) = RZ ( I ) 

RX ( I ) =RX ( I ) +DTOD*(HDTOMB*FX ( I ) + VX (I ) ) 

RY ( I )=RY( I )+DTOD*(HDTCMB*FY ( I ) + VY< I ) ) 

RZ(I )=RZ(I )+DTOD*(hOTOME*FZ ( I ) + VZ ( I ) ) 

GC TO 220 

215 RX(I)=RX< I)+DTOD*VX(I ) 

R Y ( I ) =R Y( I ) +DTOD* VY ( I ) 

R Z ( I ) =R Z ( I )+DTOD*VZ(I ) 

22C DC 24 0 1=2, LL 

IF(LCUT(I ) .EQ.2) GO TO 235 
I F ( LCUT ( I ) . GT . 0 ) GO TO 240 
J= INF ( I ) 

RXK ( I ) = RX ( I ) 

R YK ( I ) = P Y ( I ) 

R ZK ( I ) = RZ { I ) 

PX( I ) =RX ( I ) +DTOC* (HDTOM (J ) *FX ( I ) + VX ( I ) ) 

R Y ( I ) =R Y < I ) +OTOO* (HDTCM ( J )*FY ( I ) + VY ( I ) ) 

R Z ( I ) =R Z ( I ) +DTOD~ (HDTCM(J)*FZ(I )+VZ(I )) 

GC TO 240 

225 RX(I)=RX<I)+DTOD*\/X(I ) 

P Y ( I ) =R Y ( I )+DTOD*VY(I ) 
RZ(I)=RZ(I)+DTOD*VZ(I ) 

24C CONTINUE 
GC TO 2C0 

SECOND HALF OF TIMESTEP 

250 CONTINUE 
I N C E X = 0 
Ef'AX = 0.0 
T I NE=TI NE+DT 
NT = NT + 1 
1 = 1 

I F ( LCUT ( I ) . Nc • 0 ) GO TO 260 
\SS=VX( I ) 

VX( I)=VSS + HCT0M8*FX(I ) 

R X ( I ) =R XK ( I ) + ( VX { I )+VSS )*HDTOD 
V, S S = V Y( I J 

V Y ( I )=VSS+FDT0M8*FY (I ) 

RY(I)=RYK(I )+(VY( I ) + VS S )=*HDTOD 
VSS=VZ( I ) 

VZ( I)=VSS+FOTOMB*FZ(I ) 

RZ( I )=RZK( I )+(VZ ( I )+VSS )*HOTOD 
FX ( 1 ) = 0 .0 
F Y ( 1 ) =0 . 0 
F Z ( 1 ) =0 . 0 

PKE ( I ) = VX( I )*VX( I ) +VY { I )*VY{ I )+VZ( I)*VZ < I ) 
EPAX=PKE( I ) 

K = I 

P KE ( I ) = PKE ( I)*HBMAS 
26C DC 290 1=2, LL 

I F ( LCUT ( I J . NE . 0 ) GC TO 290 
J = IMF < I ) 

VSS=VX ( I ) 

VX(I)=VSS + HCTOM(J)*FX( I) 

RX(I)=RXK(I ) + (VX( I )+VSS )*HDTOD 

V SS = V Y ( I ) 

V Y ( I )=VSS+HDTOM< J)*FY( I ) 

RY(I)=RYK(I ) + ( V Y ( IJ + VSS )*HDTOD 
VSS=VZ< I ) 

VZ(I)=VSS + hDTOM(J ) ^ F Z < I ) 

RZ(I)=RZK(I)+(VZ(I ) +VS S )*HDTOD 
FX ( I ) = 0 .0 
F Y ( I )=0.0 
FZ ( I )=0 .0 

PKE( I )=VX ( I )*VX ( I )+VY (I )*VY < I J+VZ( I )*VZ< I ) 
I F ( PKE ( I) .LT.EMAX ) GG TO 280 
ERAX=PKE( I ) 
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K = I 

280 I F ( IMF ( I ) . EC. 1 ) GC TO 282 
PKE(I ) = FTMA S2* PKE (I) 

GC TO 290 

282 PKE( I )=PTMAS1*PKE ( I ) 

2 9 C CONTINUE 

TEST FOR PRINT OR TERMINATION 

CALL GcTIME(IET) 

CALL SETIME 

F IN ISH= IET *2.65-05 +FINISH 
IF(FINISH.GT.STOPER) GO TO 370 
IF ( IMF( K ) .EC.l ) GO TO 293 
E1EST=HTMAS2*EMAX 
GC TO 254 

253 ET EST = HTMA S 1*EMAX 

254 CONTINUE 

IFIETEST.LE. TESTE ) GO TO 380 
N CT =NCT +1 

IF (NCT.NE.NCTT) GC TO 360 
NCT = 0 

CC 305 1=1, LL 

IF ( LCUT ( I ) . GT. 0 ) GO TO 305 
P F E ( I ) = C . 0 
3C5 CONTINUE 

CALL ENERGY 
3C7 CC 350 1=1, LL 

I F ( LCUT (I ) . GT . 0 ) GC TO 350 

I F ( RY ( I ).LT .RYSL.ANO.VY(I) .LT.O.O) GO TO 32C 
IF (RX (I ) .LT .RXBL .OR.RX ( I ) .GT.RXBM ) GO TO 31C 
I F ( R Z ( I ).LT.RZBL.CP.RZ(I) .GT.RZBM) C-OTO 313 
IF(RY(I ) .LT.RYBM) GO TO 350 
3 1C LCUT ( I ) =1 

RXF( I) = RX( I ) 

RY F ( I ) =RY ( I ) 

R 2 F ( I ) = R Z ( I ) 

PTE ( I ) = PKE ( I) +PPE ( I ) 

GC TO 350 

32C IF(I.EQ.l) GO TO 310 

IF ( PPE ( I ) .GT. 0.0001) GO TO 350 
PTEI I )=PKE( I) +PPE ( I ) 

IF (PTE ( I) .LT.O.O ) GO TO 350 
LCLT ( I) =2 
RXF ( I ) =RX( I ) 

RYF ( I ) = R Y ( I ) 

R ZF ( I ) =RZ ( I ) 

35C CONTINUE 

TEST FOR CUT FUT OF CURRENT STATUS 

36C CT=FDTI/SQRT( EM AX ) 

IF(NT.GE.NTT) GO TO 380 

GO TO NEXT TIMSSTEP 

IFINS.NE.NT ) GO TC 100 
GC TO 4 CO 
37C IS FUT = 1 
GC TO 4 C 0 
38C N S S = 1 

IF (NT ,GE .NTT) ZSTCP=1.0 

CALCULATE POTENTIAL ENERGY AND PRINT STATUS 
OF THE CURRENT TIM ESTEP 

4CC CONTINUE 
TF CT = 0 . 0 
P FCS = 0. C 
FNEG=0 .0 
CC 450 1=1, LL 

I F ( ISEQII ) .EQ.l) I SEQ ( I ) =2 
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I F ( LCUT ( I ) . NE .0 ) GC TO 450 
PPE(I)=C.O 
P TE ( I ) = 0.0 
45C CCMINUE 

CALL ENERGY 

555 PTE(1)=PKE(1)+PPE(1) 

T F K E = PK E ( 1 ) 

6CC CC 625 1=2, LL 

TFKE=TFKE+PKE( I ) 

IF ( PPE ( I ) .GT.0.0) GO TO 615 
FNEG=PNEG+PPE( I ) 

GC TC 620 

615 P FC-S = P PCS + P Pc ( I ) 

62C PTE ( I)=PKS( IJ+PPE (I) 

625 CONTINUE 

TE=TPKE+TPOT+BIND 

IF( (NSS.NE.O) .ANC. ( ZSTOP . EQ . 0 . 0 ) ) GO TO 5 5 0 
CALL KP PINT 
WRITE (6,9681) 

7 JO CC 750 1 = 1 , LL 

IF (PTE( I) .LT.0.1 ) GO TC 750 
IF ( ISEQ (I ) .EQ.O) I SEQ ( I ) = 1 
C X ( I )=RX( I)-RXI ( I ) 

C Y ( I ) =R Y ( I)-RYI ( I ) 

C Z ( I ) =R Z ( I)-RZI (I ) 

WRITE (6, 96 82) I,OX(I),DY(I),DZ(I),VX(I),VY( I),VZ(I), 
XPKE ( I ) , PPE ( I) ,PTE( I) , RX( I ) , RY ( I ) ,RZ( I ) , LCUT ( I ) , I SEC ( I ) 
75C CONTINUE 

WFITE<6 ,96 50 TPKE , TPCT ,TE , PPOS , PNEG 
IF ( ISHUT.NE .0) GO TO 972 
IFCZSTOP.NE. 0.0) GO TO 973 
WP ITE ( 6 ,963C ) NPAGE 
NFAGE=NF AGE+1 
N S = NS + N C 



GO TO NEXT TIMESTEP 

GC TO ICO 
95C CONTINUE 

TERMINATE REPLICA ANC THEN 

PRINT FINAL STATUS OF CURRENT PASS 



J = 1 
K =C 
I K = G 

CC 560 1=1, LL 

IF (LCUT ( I ) . LT . 2 ) GO TC 960 
IF (VY ( I ) .GT .VTEST ) GO TC 96C 
K =K + 1 

I S NLM ( J ) = I 
LNUM ( I ) =LNUM ( I )+l 
J=J + 1 

5 6 C CONTINUE 
K K = K 

NSFUT(IPASS)=KK 
CALL KP P I NT 
IF(KK.EC.O) GO TO 972 
WR ITE ( 6 ,98 10 ) EV, (ISNUM(I) ,I = 1,KK) 

W R I TE ( 6 ,9600 
WRITE(6 ,9655) 

DO 963 11=1, KK 
I = ISNUM ( II ) 

WRITE (6 ,9800 IPNUM, I , BUL ,TAR , I OR ,SUR ( 1 ) , VX ( I ) , VY ( I ) , V Z 
X( I) »PKE ( I) ,FPE( I) »PTE (I ) ,EVR ,COX ,CCZ 
968 CONTINUE 

CC 569 1=2 , LL 

I F ( LCUT ( I ) . EQ.2 .AND. V Y ( I ) . L T . VTS ) WR ITE (6 ,9300 ) 

X IPNUM , I , BUL , TAR , IOR , SUR (1),VX(I),VY(I),VZ(I),PKE(I), 
X PPE ( I ) , PT 6 ( I ) , EVR 
965 CONTINUE 
C 
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C DETERMINE SPUTTERED ATOM ANGLE 

IF(KK.EQ.O) GO TO 572 
WRITE (6,9600) 

WFITE(6 ,98G5) 

DC 570 1 = 1 , KK 
L = I SNUM ( I ) 

RAC=ATAN(SCRT(VX(L)**2 + VZ(L)**2 )/ (-VY(L) ) ) 

ANC-LE = PAC*RADDEG 
J = I F I X ( (ANGLE/RINCR)+1.0) 

IF(IMF( L) .EC. I ) I EOX 1 ( J ) = I BCX I ( J ) +1 
I F < IMF (LI .EC. 2) I 80X2( J ) = I 3GX2 ( J ) +1 
WRITE (6,9815) L, ANGLE 
97C CONTINUE 

I Ft ( ISHUT.EQ.l) .CR. (ZSTOP.EQ .1.0) ) GO TC 572 
WRITE( 6,9630 NPAGE 
NPAC-E=NFAGE + 1 
GC TO 950 

972 I F ( I SHUT .E Q .0 ) GO TO 973 
WRITE ( 6 ,9965) IPNUM 
WRITE(6 ,9600) 

W R I TE ( 6 ,9960 
NSPLTt IPASS ) =0 
I FASS = I FASS-1 
N SHOT S = NSH0TS-1 
GC TO 975 

973 IF(ZSTCF.EC.O) GO TO 974 
WRITE(6,9965) IPNUM 
WRITE (6 ,96 00 

WRITE( 6 ,9955) 

NSPLT(I PASS )=0 
I F A SS = I PASS-1 
NSHCTS=NSHQTS-1 
GC TC 975 

974 WRITE (6, 5820 EV 

975 WR ITE ( 6 ,9630 ) NPAGE 
NFAGE=NPAGE+1 

DEFINE NEW REPLICA OR TERMINATE 

95 C IPASS = I PASS + 1 

IF ( ISHUT.NE.O) GO TO 1000 
EZ=6Z+SPDZ 

IF (eZ.LE.BZMAX) GO TO 40 

P X = PX + 5 PDX 

IF (BX .GT.BXMAX) C-Q TO 1000 
I F ( I0R.EQ.1C0) BZMIN=BZMIN+SPDZ 
IF(IGR.EC.lll) 3ZMIN=EZMIN+SFDZ/2 . 0 
E Z = EZMI N 

IF ( I0R.EQ.1C0) BZMAX=EZMAX-SPDZ 
IF(IOR.EG.lll) BZMAX=EZMAX-SPDZ/2»0 
GC TO 40 

PRINT SUMMARY OF ALL SHCTS 

1000 CONTINUE 

CALL KPRINT 
WRITE (6,9600) 

NRAN=IPASS- 1 

WRITE( 6 ,986C) NRAN , EVR 

WR ITE (6 ,9870 (( IPNUMS(I) ,NSPUT(I) ) ,1 = 1, NRAN) 

K S L M=0 
K 2SUM = 0 

DC 1005 1 = 1 , NRAN 
KSLM = KSLM + NSPUT( I ) 

K2SLM = K2SUM+NSPUT (I )**2 
1005 CONTINUE 
FRAN=NR AN 
FRANR = 1 .O/FRAN 
AVE=FLOAT (K$UM) *FRANR 
VAR=FLCAT (K2SUM) *F R ANR-AVE**2 
STD=VAR*FRANR 
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STC = SQR T ( STC ) 

WRITE (6,9600) 

WPITE(6,9875) ETHR , VTEST ,KSUM 
WRITE(6 ,9850) AVE ,VAR , NRAN , STD 
WRITE ( 6 » 96C0 ) 

WRITE (6,9865) 

WPITE(6 ,9870 ( I , LNU M < I ) , I =1 , LL ) 

WRITE (6,9600) 

WRITE ( 6 , 9S 20 ) IREAO 

WRITE (6,9925) ( I FFM( J ) , I FM ( J ) , J = I , IREAO ) 

IF ( I0R.NE.100) GO TO 1007 
WR ITE ( 6 ,9950 ISEEDA 
10C7 IF(IOR.NE.llO) GO TG 1008 
WRITE (6 ,9950) ISEECB 
1008 IF ( IOR.NE.l 11 ) GO TO 1009 
WR ITE(6 ,9950 ISEECC 
LO 09 CONTINUE 

W R ITE ( 6 » 9630 ) NPAGE 
NFAGE=NPAGE+1 
CALL KPNT 

FOLD FINAL SUMMARY 

IF (IOR.EQ.lOO) GC TO 1011 
IF (I OR. EQ. 110) GC TO 1012 
IF (I OR. EQ. Ill) GC TO 1013 

1011 CALL KFCLDA 
C-C TO 1C 14 

1012 CALL KFCLDB 
GC TO 1014 

1012 CALL KFCLOC 

1014 WRITE (6,9600) 

WRITE (6,9840) 

WRITE (6,9845) 

WPITE(6 ,9870 ( I , LNUM ( I ) , I = 1 , LL ) 

WRITE (6,9600) 

DETERMINE RING SPUTTERING PROBABILITY 

1015 CALL KPPCB 

WR ITE ( 6 ,9890) 

WRITE(6 ,9860 
CC 1020 1=1 »KR I NG , 3 
J = I + 1 
K = I + 2 

WRITE(6 ,9885) I , DNN ( I ) , K ATOM ( I ) , K S PUT ( I ) , PRO e ( I ) , J , 

X DNN ( J) ,KATOM(J) ,KSPUT (J) , PROB ( J) , K , DNN ( K ) ,K ATOM ( K ) , 
X KSPUT (K) , PROB (K ) 

1020 CONTINUE 

WRITE ( 6 , 96 CO ) 

WRITE (6,9620) NPAGE 
NF AGE = N PAGE+1 
CALL KPNT 
WRITE ( 6 , 9 6 CO ) 

DRAW ISOTOPE HISTOGRAMS 

1025 CC 1030 1= 1 , 1 RANGE 
FPEC1 ( I )= IBCX1 ( I ) 

F F EG2 ( I ) = I BCX2 ( I ) 

1020 CONTINUE 

RANGE ( 1 )=RM IDPT 
CC 1040 1 = 2 , 1 RANG E 
J=I-1 

RANGEt I )= RANGE (J ) + RINCR 
104C CONTINUE 

IN =-I RA NGE 

CALL H I STO ( IN, FREG1, RANGE, TITLED 
WRITE l 6 , 9 6 CO ) 

WRITE (6,9895) 

WRITE (6,9900) RINCR 

IF ( IRANGE.EQ.9) WRITE (6,9925) 
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IF ( I RANGE .EQ.I8 ) WR I TE ( 6 , 9920 ) 

WRITE (6,9910) ( IB0X1(J) ,J=1 , IR ANGE ) 
WRITE (6,9600) 

WRITE ( 6 , 96C0 ) 

CALL HI STO ( IN, FREC2, RANGE, TITLE 2) 
WRITE (6,9600) 

WRITE (6,9895) 

WRITE (6,9900) RINCR 

IF ( IRANGE.EQ.9) WRITE (6,9925) 

IF ( IRANGE .EQ.18 ) WRITE (6 ,9920 ) 

WRITE (6,9910) ( I BGX2 ( J ) ,J= 1 , IRANGE ) 

WRITE (6,9600) 

WRITE (6,9620) N RAGE 
NFAGE=NFAGE+1 
IF(ISHUT.EC.O) GO TO 5 
9999 STOP 
E f' C 
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SUBROUTINE SETUP 



C 

C 

C 

C 

C 

C 



THIS S/R FINDS THE STARTING POINT FCR EACH 
IMPACT POINT IN THE TARGET AREA 



CCMMCN/CQM1/RX( 60 C ) , R Y ( 600 ) , RZ ( 6 0C ) ,LCUT (60C ) » LL , LD 
CCNMCN/CCM2/RCE, RCE2, RQEM, AC, PAC , P PTC , P TC , P FPTC ,FPTC , 
X FN,PFIV,TPCT 

CCMMCN/C0M6/C0XI , COY I ,CCZI , RX S , R Y S ,RZS , F AC 
CCMMCN/CCM7/R1, LSS ,SPX,SPZ ,COX ,CCY ,CCZ,YLAX (23 ) , I LL 
CGMMCN/COMii/RXI (60C),RYI(600) , RZI (600) 

1 = 1 

100 R 1 = R 1-F AC 

IF(Rl.LT.O.C) GO TO 300 
R 1X=R1-CGXI +RXS 
F 1Y = R1*CCYI+RYS 
F1Z=R1*CCZI+RZS 
X M IN = R IX- POE 

Y N IN = R1Y-P0E 
ZM IN = R1Z-R0E 
X M AX = F1X + R0E 

Y MAX = R1Y + RQE 
ZM AX = R1Z+P0E 
CC 195 J = 2 » LL 

IF (LCUT ( J ) . GT • 0 ) GC TO 195 
I F (RY ( J ) .GE .YMAX ) GC TO 195 
IF (RY( J ) .LE .YMIN) GO TO 195 
IF (RX(J ) .GT.XMAX) GO TO 195 
IF(RX(J).LE.XMIN) GC TO 195 
IF(RZ(J ) .GE.ZMAX) GO TO 195 
IF(RZiJ).LE.ZMIN) GO TO 195 
C F X =RX ( J ) - R IX 
C R Y=R Y ( J )-RlY 
CP Z=RZ ( J )-R 1Z 
D I ST =0 R X*^ 2 +DR Y^ 3 * 2 + DRZ *^2 
I F (CIST .GE.R0E2) GO TO 195 
GC TC 2 C 0 
195 CONTINUE 
GC TO ICO 
2 0 C CONTINUE 

R2SQ= (RXI(J )-RXS )*#2+(RYI (J) -RYS ) **2+ (RZ I < J )-RZS) **2 
R2=SCRT (R2SC ) 

R 1 P 2= ( P 1X-RXS)* ( RXI ( J ) -RXS J + ( R 1 Y-RYS ) * ( R Y I ( J ) -RYS ) + 

X (P1Z-RZS ) * (RZI ( J )-RZS ) 

A LFA = R1 F2/ ( R1*R2 ) 

A LFA2 = A LFA* AL FA 

R2=R2*ALFA+SQRT(R2SQ*ALFA2-R2SQ+RCE2) 

RX (1 )=R3*CCXI+RXS 
R Y ( 1)=R2*CCYI+RYS 
RZ(1)=R3*CCZI+RZS 
RETURN 

3CC RX ( 1 ) =R XS 
R Y ( 1 ) =R Y S 
RZ( 1) =RZS 

97CC FC R M AT ( 4 1 H THIS DIC NCT LOCATE A PROPER START POINT) 
ie P IT E ( 6,9700) 

RETURN 

ENC 



SIEROUTINE hARM 



C 

C 

c 

c 
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8Y ITS THERMAL AMPLITUDE 



CCMMGN/C0M1/RX(6G0) ,RY(60J) , RZ ( 6 00 ) , LCUT (60 C ) ,LL,LD 
CCMMCN/C0M13/ RAN ( 1800) , I RAN , T A SN , TASF , L L2 , LL3 
COMMON/ COM 14/ AXT (6 CO) , AYT(600) ,AZT(600) 

A XT ( 1 ) = 0 
AYT ( 1 ) =C 
A 2T ( 1 ) = C 

CALL NCFMAL (IRAN, RAN, LL3) 

CC 1100 1 = 2, LL 
AXT ( I ) =RAN ( I ) 

AVT(I)=FAN( I + LL) 

AZT ( I ) = R AN ( I+LL2 ) 

P X { I } =R X ( I ) +TASN*AXT ( I ) +T AS F*AX T ( I- 1 ) 

RV ( I )=RY( I ) +TASN* A YT { I )+TASF*AYT ( 1-1 ) 

R 2 < I ) =R Z ( I )+TASN*AZT( I )+TASF*AZT{ 1-1) 

1100 CONTINUE 
RETURN 
END 



SLERCUTINE KF100 



LATTICE GENERATOR FOR A FCC (100) ORIENTATION 
LATTICE IS DEVELOPED IN THE X DIRECTION, THEN 
THE Z DIRECTION, THEN FINALLY THE Y DIRECTION 



CCMMCN/CCM1/RX(60C) , R Y ( 600 ) , RZ ( 60C ) , LCUT (60 0 , LL , LD 
CCMMCN/C0M2/RCE, RCE2 , RGEM , AC , PAC , PPTC ,PTC, P FPTC ,FPTC, 

X FM,PFIV,TPCT 

CCMMCN/C0M4/IX, IY , IZ , IX P , I Y P , SC X , SCY , SC Z , I C E E P , 

, X F EX, RE Z 

COMMON /COM 7/Rl,LSS,SPX,SPZ, COX, COY, CCZ,YLAX (20) , ILL 
CCMMCN/C0M11/RXI (600),RYI(600),RZI(600) 

CCMMCN/C0M23/ SP D X , SPDZ ,U3XN AX , C E ZMAX , SSC Z 
COMMON/ COM 2 5/ IMF (600) , ISOGEN 
CCMMCN/CGM3 0/ ISEECA, I SEEDS , ISEECC 
DIMENSION R (600 ) 

9001 FC FMAT ( ' I S EED= ' , 1 10/ / ) 

9690 FCFMAT( 35X , 'SUMMARY OF ATOMS FCC (100) ORIENTATION' 

X / , 3 5X , 'MASS 1 IS THE HEAV I ER /L E AST ABUNDANT ISOTOPE' 
X / / , 3 ( ' ATOM MASS POSITION '),/) 

LC = C 
N =600 

CC 3 1=2, N 
I M F ( I ) = 2 
2 CONTINUE 

I $ EED= I SEED A 

WRITE (6,9001) ISEED 

WRITE (6,9690) 

IF (ISOGEN. EQ.O) GC TO 4 

GENERATE RANDOM MASSES 

CALL GGU3 ( ISEED, N, R ) 

CC 4 1=2, N 

IF(R(I) .LE.C.309) IMF ( I ) = 1 

4 CONTINUE 

CC 5 1 = 3, IY 

5 YLAX( I ) =0.0 

DETERMINE ORIENTATION PARAMETERS 

SCX = 1 .0 
S C Y = 1 « 0 
SCZ=1.0 
S F CX = 0 . 1 
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20 

27 

56 

59 

6C 

7 C 

8C 

82 

85 

90 

ICC 

1 1C 

120 



S F C Z = 0 . 1 
C E X N AX= 0 .5 0 1 
C E ZHA X= 1 . 00 1 
YLAX( 1) = -0. 12 
YLAX(2) =-0.03 

GENERATE THE LATTICE 

N = 2 
J T =0 
Y =-SC Y 

DC 60 J = l, IY 
Y=Y+SCY 
K7=C 
Z=-SCZ 

CC 59 K =1 , I Z 
Z=Z+SCZ 
I T =0 
X=-SCX 

CC 58 1 = 1, IX 

X=X+$CX 

I TT = IT+ JT + KT 

I F ( ITT — ( ITT/2) *2) 57,30,57 
PX I (M > = X 

RYI (M)=Y + YLAX( J ) 

R Z I ( M ) = Z 
H = H + 1 
CONTINUE 
IT = IT+ 1 
CCNT IMS 
KT=KT+1 
CCNT INUE 

I F ( JT .EC.OJ I L L = M - 1 
JT = JT + 1 

IF(JT.6C.I DEEP ) LC=M-1 
CONTINUE 
L L = M— 1 

IF(LD.EC.O) LD=LL-1 
IF(COY.EQ.O) GO TG 120 
C-C TC ( 70,80,90) ,LSS 
R 1=RCE/A6S (COY) 

GC TO 100 
I S F X = R8 X 
I SFZ=RE Z 

N=( IX/2 )^ISFZ+ISPX/2+3 
I F ( IX- ( IX/ 2 )*2 ) 83,85,83 
N=N+I SPX/2 
LCLT ( N ) = 1 
R 1 = 2 .0 
LCLT ( 2) =-L 
GC TC 100 
N = 2 

RX I (M)=R8X*$CX 
RYI (M)=-SCY+YLAX( 1) 

PZI <M)=(PBZ+1.0)*SCZ 

R 1= (ROE+SCY J/A3S ( COY ) 

CONTINUE 

RETURN 

R 1=1.5 

RETURN 

ENC 



SLEPOUTINE KF110 



LATTICE GENERATOR FOR A FCC (110) ORIENTATION 
LATTICE IS DEVELOPED IN THE X DIRECTION, THEN 
THE Z DIRECTION, THEN FINALLY THE Y DIRECTION 
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c 

CCMMCN/CGM1/RX(60C) ,RY(600) , RZ ( 6C C ) , LCUT ( 60 C ) ,LL,LD 
CCMMCN/CCM2/RCE, RCE2,R0EM, AC, PAC, PPTC,PTC, PFPTC,FPTC, 

X FM,PFIV,TPCT 

CCNMCN/C0M4/IX, I Y, I Z, IX P, I YP, SC X, SC Y, SC Z, I DEEP, 

X F EX , RE Z 

CCNMCN/ CCM7/R1 ,LSS,SPX,SPZ,C0X,CCY,CQZ,YLAX(20),ILL 
CCNMCN/COM11/RXI ( 60 0 ) , R Y I ( 6 GO ) , R Z I ( 600 ) 

CCMMCN/C0M23/ S? DX , SPDZ ,D8X MAX , DB ZMAX , S SCZ 
CCMMON/COM25/IMF( 600) , I SO GEN 
CCMMCN/CQM3G/ ISEEDA, ISEE08, ISEEDC 
CIMENSICN R (600 ) 

9001 FCFMATC ISEED=' , 1 10// ) 

965C FORMAT ( 35X » 'SUMMARY OF ATOMS FCC (110) CR I ENT AT I ON ' 

X / ,35X , ' MASS 1 IS THE H6AV I ER /L E AST ABUNDANT ISOTCPE’ 
X / / , 3 ( * ATOM MASS POSITION '),/) 

LC = C 
N = 6C0 

CC 3 1=2, N 
IN F ( I ) = 2 

3 CONTINUE 
ISEED=ISEECS 

WRITE ( 6 , 90C1 ) ISEEO 
WRITE (6,9690) 

IF ( IS0G6N.EQ.0) GC TO 4 
C 

C GENERATE RANDOM MASSES 

C 

CALL GGLB ( I SE ED , N , R ) 

CC 4 1=2, N 

I F ( R ( I ) .LE.C.309) IMF( I ) = 1 

4 CONTINUE 

CC 5 1 = 3, IY 

5 YLAX(I)=O.C 
C 

C DETERMINE ORIENTATION PARAMETERS 

C 

P C = 1 . 0/ SQRT (2.0) 

SCX=RG 

SO =R0 

SCZ=1.0 

SFCX=0. 2*SCX 

SFCZ=0.2*SCZ 

CEXNAX=$CX+0.01 

C EZMA X= SC Z + C . 0 1 

Y L AX ( 1 ) =-0 . 15 
YLAX(2) =-0.03 

C 

C GENERATE THE LATTICE 

C 

M = 2 
J T = 0 

Y =-SCY 

CC 60 J«l, IY 

Y = Y+SCY 
KT = C 
Z=-SCZ 

CC 59 K = 1 , I Z 
Z=Z+SCZ 
I T = 0 

x=-scx 

CC 58 1=1, IX 
x =x + 9r x 

IF ( I T— ( IT/2) *2) 21,11,21 

11 IF (JT-< jT/2)*2) 57,12,57 

12 IF (KT-(KT/2 )*2) 57,3C,57 

21 IF(JT-(JT/2)*2) 22,57,22 

22 IF(KT-(KT/2 )*2) 30,57,30 
30 PX I (M) =X 

RY I (M) = Y + YLAX(J ) 

R Z I (M ) = Z 
M = M + 1 
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57 CCNTINUE 
IT =1 T+ 1 

58 CCNTINUE 
KT=KT+1 

59 CCNTINUE 

JT = JT + 1 

IF ( JT.EC.IDEEP) LC=M-1 

60 CCNTINUE 
LL=M-1 

IF(LD.EC.O) LD=LL- 1 
I F ( COY . cQ . 0 ) GO TO 120 
GC TO ( 70*80,90) ,LSS 
7 C R 1=RCE/ABS (COY) 

GC TO 100 
80 ISFX=SPX 
I SPZ=SP Z 

N = ( IX+1 )/2*ISPZ/2+ISPX+2 
LCLT(N) =1 
R 1=RCE/A8S (COY) 

LCLT ( 2 ) =-l 
GC TO 100 
SC M = 2 

R X I (M)=(SPX+1.0)*SCX 
RY I (M) = -SCY+YLAX( 1) 

RZI <M) = (SPZ+1.0)*SCZ 
R1=(PGE+SCY )/A8S(CQY) 

ICC CCNTINUE 
11C RETURN 
12C R 1 = 1 . 5 

RETURN 

ENC 



SLERCUTINE KF111 



LATTICE GENERATOR FCR A FCC (111) ORIENTATION 
LATTICE IS DEVELOPED IN THE X DIRECTION, TI-EN 
THE Z DIRECTION, THEN FINALLY THE Y DIRECTION 



COMMON/ COM 1/RX( 60 C ) ,RY (600 ) ,RZ ( 6CC ) ,LCUT ( 60 C ) , LL , LD 
CCMMCN/CCM2/PCE, RCE2,RGEM, AC, PAC ,PPTC,PTC, PFPTC,FPTC , 

X FM ,PFIV,TPCT 

CCMMCN/CQM4/IX,IY , IZ , IXP, IYP,SCX,SCY,SCZ,IQESP, 

X R E X , RB Z 

CCMMCN/CCM7/R1 ,LSS, SPX, SPZ.CCX ,CCY ,COZ , Y LAX ( 20 ) , I LL 
CCMMCN/C0M11/RXI (600),RYI(600),RZI(600) 

CCNMCN/C0M23/ SP DX , SPD Z , DBXM AX , D 6 ZMAX , SS C Z 
COMMON/ COM 2 5/ IMF ( 6C0) , ISOGEN 
CCMMCN/CCM3 3/ ISEEDA, ISEED8, ISEECC 
DIMENSION R ( 600 ) 

9001 FORMAT ( ' IS EED= ' , 1 10/ / ) 

9693 FCFMAT( 35X, 'SUMMARY OF ATOMS FCC (111) ORIENTATION' 

X / ,35X , ' MA SS 1 IS THE HEAV I ER /L E AST ABUNDANT ISOTCPE' 
X / / , 3 ( ' ATOM MASS POSITION '),/) 

LC = C 
N = 600 

DC 3 1=2, N 
INF C I ) = 2 
3 CCNTINUE 
ISEEC=I SEECC 
WRITE (6,9001) ISEED 
WRITE (6,9690) 

IF (ISOGEN. EQ.O) GO TO 4 

GENERATE RANDOM MASSES 



CALL GGL 8 ( ISEED, N,R) 
CC 4 1=2, N 



non ooo 



A 

5 



30 

A 1 
A2 

A3 



A5 

50 

57 

5 £ 

5 9 

60 



7 C 
80 



90 



I F ( R < I ) .LE. 0.309 ) IMF(I) = 1 

CONTINUE 

DC 5 1 = 3, IV 

Y LA X ( I ) =0.0 

DETERMINE ORIENTATION 

SCX=1.0/SQRT(2.0) 

SCY=2.0/SQRT(3.0) 

SC2=SQRT (1.5) 

SSCZ=2.0*SCZ/3.0 
SFCX=0. 2 *SC X 
S F CZ=0 . 2 *S SCZ 
CEXMAX=SCX+C.01 
CEZMAX=SSC2+0.01 

Y LAX ( 1 ) =-0 . 05 
YLAX (2 ) =-0 . Cl 

GENERATE ThE LATTICE 



Y = -SCY 

DC 60 J=l, IY 

Y = Y + SC Y 

J T S = JT+ JT/ 3 
Z=-SCZ 

1/T-O 

CC~59 K = I , IZ 
Z=Z+SCZ 
I T= C 
X=-SCX 

DC 58 1=1, IX 

X=X+$CX 

I N = IT + J IS + KT 

IF( IN-( IN/2)*2) 57,30,57 
PX I (M ) = X 

RYI (M)=Y+YLAX(J ) 

I F ( JT-3*( JT/3) ) A 1 , A 5 , A I 

JTT= JT 

JIT =JTT-3 

IF(JTT) A3 , A5 , A2 

JTT = JTT +3 

ZF= JTT 

RZI (M ) = Z + ZP*SSC Z 

GC TO 50 

RZI (M) = Z 

N=M + 1 

CONTINUE 

I T = I T+ 1 

CCNT INUE 

KT =KT + 1 

CONTINUE 

JT= JT+ 1 

IF(JT .EG.IDEEP) LC=M-1 

CCNT INUE 

Ll=M-l 

IF(LD.EC.O) LD=LL- 1 
IF(COY.EQ.O) GO TO 120 
GC TC ( 70,80,9 0) ,LSS 
R1=RCE/A6S(C0Y) 

LCLT ( 2 ) =- 1 
GC TO 100 
ISFX=SPX 
I SF Z = SP Z 

N = IX*( I SPZ/2J + ISPX + 2 
LCL T ( N ) = 1 
R1=RCE/ABS (COY) 

GC TC 100 
M = 2 

RX I (M)=$PX*SCX 
RYI (M)=-SCY+YLAX( 1) 



PARAMETERS 
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RZI (M)= (SPZ+1.0)*SCZ 
Rl= (ROE + SCY ) / AB S ( COY ) 
ICC CCNTINUE 
11C RETURN 
1 20 R 1= 1 . 5 
RETURN 
END 



SUBROUTINE ENERGY 



THIS S/R CALCULATES THE MUTUAL POTENTIAL ENERGIES 
THIS S/R ALLOWS ATTRACTIVE POTENTIALS 
THIS IS THE MUNOKUR HIGH SPEED VERSION 



CCNMCN/CCM1/RX(60C) ,RY<600) ,RZ ( 600 ) , LCUT (60C 1 ,LL,LC 
CCMMCN/CCM2/RCE, RC E2 , ROSM , A C , PAC , P PTC , P T C , P FPTC ,FPTC, 

X FM,PFIV,TFGT 

CCMMCN/COM3/EXA ,EXB,FXA,PEXA,PEXE,PFXA,TFERN 
COMMON/ COM5/FX( 60 0 1 ,FY(600) , FZ ( 6 00 ) , PPE ( 6 00 ) 
CCMMGN/COM22/ROEA,RCEB,ROEC , ROEC 2 , CPO ,C P 1 ,CF2,CP3 ,CF0 , 
X CF1,CF2,CGD1,CGC2,CGB1,CGB2,CGF, I PROG 
1 = 1 

IF (LCUT ( I) .NE.O) GC TG 200 
I F = I + 1 

INITIAL POTENTIAL ENERGY INTERACTION 
MAXIMUM DISTANCE BOUNDARIES 

XMAX = RX ( I ) +ROE 
YMAX = RY ( I ) +ROE 
ZM AX = RZU1+R0E 
XMIN = RX ( I l-ROE 
YMN = RY (I l-ROE 
ZMIN = R Z ( I l-ROE 

TEST TO DETERMINE IF BULLET AND ATOMS 
ARE CLOSE ENOUGH TO INTERACT 

DC 195 J=IP,LL 
I F ( LCUT ( J 1 . GT . 0 1 GC TO 195 
IF (RY< J ) .GE .YMAX 1 GO tq 195 
IF(RY(J l.LE.YMIN) GO TO 195 
IF(RX(J).GT .XMAX ) GC TO 195 
IF(PX(J l.LE.XMIN) C-0 TO 195 
IF(RZ( Jl.GE.ZMAX) GC TO 195 
IF{RZ(J).LE.ZMIN) GO TO 195 

TEST TO DETERMINE IF THE ATOMS ARE INSIDE A 
CIRCLE OF RADIUS ROE FRCM THE 8ULLET 

C R X =RX ( J ) — R X < I) 

CRY=RY(J)-RY( I ) 

CPZ = RZ(J)-RZ( I) 

C I S T = C R X* * 2 +0 R Y *■ * 2 + D R Z * *2 
I F ( D I S T .GE. R0E2 1 GC TO 195 
C IST = SQPT (CIST) 

CALCULATICN OF THE POTENTIAL ENERGY 

PCT=EXP (PEXA+PEXB^OIST l-PPTC 
IF (POT .LE .C .00001 ) POT = 0.0 
T FCT = TPCT + PCT 
SA VE = PF IV*PCT 
PPE < I ) = PPE ( Il+SAVE 
PPE ( J) =PFE( J) +SAVE 
195 CCNTINUE 
2C0 CCNTINUE 
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CC 300 1=2, LD 

I F (LCUT ( I ) . GT.O ) C-C TO 300 
I F = I + 1 

XMAX = RX ( I )+R0E 
XMIN = RX (I J-ROE 
YMAX = R Y ( I ) + R0E 
YMIN = RY (I )-R0E 
ZMAX = RZ (I ) + R0E 
ZNIN = RZ < I )-R0E 

TEST TO DETERMINE IF THE ATOMS 
ARE CLOSE ENOUGH TO INTERACT 

CC 295 J=IP,IL 
IF(LCUT{J) .GT.O) GC TO 295 
I F ( RY ( J ) . GE .YMAX ) GO TO 295 
IF ( RY ( J ) .LE.YMIN) GG TO 295 
IF(RX(J ) .GE .XMAX ) GC TO 295 
IF(RX(J ) .LE.XMIN) GG TO 295 
IF (RZ(J ) .GE.ZMAX) GC TG 295 
IF(RZ( J ) .LE.ZMIN) GO TO 295 
CRX=RX( J)-RX( I) 

DRY=RY ( J )-RY( I) 

CRZ=RZ( J )-RZ( I ) 

0 IST=DRX**2+CRY**2+CRZ**2 
IFtCIST .GE.ROE2) GC TO 295 
CIST=SQRT (DIST) 

CALCULATICN OF THE FOTENTIAL ENERGY 

PCT = EXP ( EXA+EXB*DIST)-PTC 
IFtFOT.LE. 0.00001) P0T=0.0 
28 C TFCT = TPCT + PCT 
SAVE=PF IV*PCT 
PPE ( I ) = PPE ( D+SAVE 
PFE(J)=PPE( JJ+SAVE 
295 CONTINUE 
3CC CONTINUE 

DC 313 1=1, LL 

I F ( LCUT ( I ) . EQ . 0 ) GC TO 310 
T FCT = T P CT + P PE ( I ) 

310 CONTINUE 
RETURN 
ENC 



SLERCUTINE STEP 



THIS S/R CALCULATES THE FCRCES BETWEEN ATOMS 
ATTRACTIVE FORCES ARE ALLOWED BETWEEN ATGMS 
THIS IS THE MLNDKUR HIGH SPEED VERSION 



CCMMGN/COMl/RXI 60 0 ,RY(600) ,RZ(fc0C) ,LCUT(6CC) , LL , L D 
CCMMCN/CCM2/RCE, RCE2,RCEM,AC,PAC, FPTC,PTC, PFPTC ,FPTC, 

X FM , PF I V , T PCT 

CCMMCN/C0M3/EXA , EXB , FXA , P EX A , P E XB , PFXA , THERM 
CCMMCN/CQM5/FX(60C) ,FY(600) ,FZ(600),PPE(600) 

COMMON/ COM2 2/ ROE A ,ROcB,RQEC , ROEC 2 , CPO ,C P 1 , C F2 , C P3 ,CFO , 
X CF1,CF2,CGD1,CGC2,CGB1,CGE2,CGF,IPR0G 
1 = 1 

I F ( LCUT ( I ) .NE.O) GO TO 200 
I F = 1 + 1 

INTIAL FORCE DISTANCE BOUNDARIES 

XMAX = RX ( I J+ROE 
YMAX = RY ( I ) + ROE 
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ZMAX = RZ ( I ) + ROE 
XMN = RX(I)-ROE 
YMN = RY ( I )-ROE 
ZMN = RZ (I )-ROE 

TEST TO DETERMINE IF BULLET AND ATOMS ARE 
CLOSE ENOUGH TO INTERACT AND WHICH FCRCE 
FUNCTION TO USE 

CC 195 J= I Pt L L 
I F ( LCUT ( J ) . GT . 0 ) C-C TO 195 
IF ( RY ( J ) . GE .YMAX ) GO TO 195 
IF (RY(J).LE.YMIN) GO TO 195 
IF(RXtJ).GT.XMAX) GO TO 195 
IF(PX(J ).LE.XMIN) C-0 TO 195 
IF (RZ( J ) .GE.ZMAX ) GO TO 195 
IF(RZ(J ) .LE.ZMIN) GO TO 195 

TEST TO DETERMINE IF AT CM IS INSIDE A CIRCLE 
OF RADIUS ROE FROM THE BULLET 

CBX=RX( J)-RX( I) 

DP Y=RY ( J ) — R Y ( I ) 

CRZ=RZ( J )-RZ( I ) 

DIST=DRX**2+DRY**2+DRZ**2 
I F(DIST .GE.R0E2 ) C-0 TO 195 
C I S T = SQ RT ( C 1ST ) 

IF (DIST .GT.BCEM) GO TQ 165 

ACTUAL FORCE FUNCTION DETERMINED FROM THE 
BCRM-MEYER BULLET-ATOM POTENTIAL 

FCRCE=EXP(F FXA + PEXB^D 1ST ) 

GC TO 189 
165 CFF=ROE-DIST 

IF(DFF.LT.l.OE-lO) GO TO 195 , 

AN AVERAGE FORCE FUNCTION FOR THE BLLLET-ATCM 
INTERACTION AT TURNON 

FCRCE = ( EXP ( FAC+PEXB*DIST)-PFPTC) /OFF 
IFtFM.GT. FORCE) GC TO 195 
18 C FCC =FCRCE/CIST 
FA = FCO*DRX 
F X ( J ) = FX(J) + FA 
FMI) = FX( I ) — F A 
FA = FCC*DRY 
F Y ( J ) = FY ( JJ+FA 
FY ( I ) = FY ( I)-FA 
FA = FCD*DRZ 
F Z ( J ) = FZ ( J ) +F A 
F Z ( I ) = FZ ( I ) - FA 
195 CONTINUE 
2CC CONTINUE 

CC 200 I =2 t LD 

I F ( LCUT (I ) . GT. 0 ) GO TO 300 
I F = 1 + 1 

XNIN = RX ( I )-ROE 
XMAX = RX(I)+ROE 
YM IN = RY ( I )-ROE 
YN A X = RY ( I ) + ROE 
ZMN = R Z ( I )-ROE 
ZM A X = RZ (I J + R05 

TEST TO DETERMINE WHETHER ATOMS ARE CLOSE 
ENOUGH TO INTERACT AND WHICh FORCE FUNCTION 
TO USE. 

DC 295 J=IPtLL 
I F ( LCUT ( J ) . GT . 0 ) GC TO 295 
IF (RY (J ) .C-E .YMAX) GC TO 295 
IF(RYU).LE.YMIN) GO TQ 295 
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265 



28 C 



295 

3CC 



96 1C 
964C 

9645 

965C 



IF (RX(J ) .GE .XMAX) GC TO 295 
IF(RX(J ) .LE.XMIN) GO TO 295 
IF (RZ(J ) .GE.ZMAX ) GO TO 295 
IF (RZ( J ) .LE.ZMIN) GO TO 295 
CRX=RX ( J )-R X ( I ) 

C F Y = RY ( J ) — R V ( I) 

CRZ=RZ( J)-RZ(I) 

D IST=DRX**2+DRY**2 + DRZ**2 
I F ( C I ST .GE. R0E2 ) GC TO 295 
D I S T = SQRT ( D 1ST ) 

IF(CIST.GT.ROEM) GC TO 265 

ACTUAL FORCE FUNCTION DETERMINED FR CM THE 
CU-CU <Gie-II) POTENTIAL FUNCTION 

FCRCE=EXP(FXA+EX9*DIST ) 

GC TO 280 
G FF=RG£-D I ST 

IFIDFF.LT . I .0 E - 1 0 ) GO TO 295 

AN AVERAGE FORCE FUNCTION FOR THE A TCM-AT CM 
INTERACTION AT TURN ON 

FCRCE=(EXP(AC+EXB*C IST)-FPTC)/DFF 

IFtFM.GT. FORCE) GC TO 295 

FCC=FGRCE/D 1ST 

F A =FOD* CRX 

F X ( J ) = F X ( J ) +FA 

F X ( I ) = F X ( I)-FA 

F A = FCD*CR Y 

FY ( J ) =F Y ( J ) +FA 

FY(I)=FYU)-FA 

F A = FOD*DRZ 

F Z ( J ) = F Z ( J ) +FA 

F Z ( I ) = F Z ( I ) -FA 

CONTINUE 

CONTINUE 

RETURN 

ENC 



SIEROUT INE KPRINT 



THIS SUBROUTINE PRINTS THE PAGE FEACING 



CCMMCN/C0M1/RX(600) , RY(600),RZ(6CO),LCUT(60C),LL,LD 
COMMON/ COM3 /EXA, EXB , FXA , P EX A , P Exe , PFXA , T HER M 
COMMGN/COM4/IX, IY,IZ,IXF,IYP,SCX,SCY,SCZ,IOEEP, 

REX,REZ 

COMM ON/ C0M7/R1 , L SS , SPX , SPZ , CCX , CCY , COZ, YLA X ( 2 0 ) , ILL 
CCMMCN/COMI 1/RXI (6C0),RYI(600),RZI (600) 

COMMON/ C0M1 6/BMAS ,TMAS1 ,TMAS2, A BUN 1,ABUN2, I FNUM 
CCMMCN/C0M16A/BX, EY,BZ, EDISF,CVP 

CCMMCN/COMI 7/TIME ,NT , DT ,DT I ,TEMF , CELS , E V , NP AG E , NL I NE 
CCMMCN/C0M18/ IT EMP , T AM P , CORF 

CCMMCN/COM20/IHI (20) , I HS ( 10 ) , I HT ( 6 ) , IHB ( 6 ) , IHTARG (2 ) , 
IFBUL L ( 2 ) 

CCMMCN/C0M21/ IHBPCT ( 2 ) , IHT POT ( 2 ) t SUP ( 2 ) , IH2 ( 1 0 ) , I OR , 
LSLR 

CCMMCN/C0M2 2/R0EA ,RC£B,ROEC , ROEC 2 , C FO , CP I ,CF2,CP3 ,CFC , 
CFl,CF2»CGCl,CGD2,CGBltCGB2*CGFl,CGF2» IPRCC- 
FORMAT ( 2CX,20A4, / ,4CX, 10A4, ///) 

FORMAT ( 1 (',13,') PLANE ',2A4,' SURFACE, PRIMARY', 

• ENERGY = • , F7 • l , ' EV, CRYSTAL SIZE (',12,' X ',12, 

• X ',12,' ), IMPACT POINT ',19/) 

FORMAT ( ' TARGET -',2A4,20X,' BULLET -',2A4,7X, 
•LATTICE UNIT =',F7.4,' ANG') 

FORMAT ( ' MASS1 =', F7 . 2 , 2X ,' ABUNDANCE =',F5.2,6X, 
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X 'MASS =• ,F7. 2 ,9X, • LATTICE TEMPERATURE = • , A4 , 

X ' THERMAL AMPLITUDE =',F7.4) 

9652 FCPMAT ( ' MASS2 =' , F7 .2 ,2X , ' ABUNCANCE =',F5.2, 

X 29X, 'CORRELATION =',F5.3,//) 

9655 FORMAT ( 12X,2(2A4 , IX) , 2A4 , 3X , 5 HP EX A = , F9 . 5 , 2X , 5HPEXB = , 

X F9.5, 2X,5HPFXA= ,F9.5/ ) 

966C FORMAT ( 12H POTENTIAL , 2( 2A 4 , IX ) , 2A4 , 3X , 5HE XA =,F9.5, 

X 2X , 5H EXB =,F9.5 ,2X,5HFXA =,F9.5) 

9665 FORMAT ( BOH TARGET POINT ON CRYSTAL X =,F7.3,', Y = ', 
X F7 .3 » ' » Z =',F7.3,13H, COS TO X =,F8.5, 

X 12H , CCS TO Z = , F8 • 5 ) 

967C FG F MAT ( 3 OH PRIMARY START POINT (LU) X =,F7.3,5H, Y = , 
X F7.3t5H, Z = » F7 . 3/ ) 

9675 FCPMAT ( / ' TIMESTEP *,I4,40X,’ ELAPSED TIME (SEC) = ', 

X E10.4,', LAST TIMESTEP WAS =',E10.4/) 

9725 FORMAT ( ' CUT-OFF AT',F5.2,', WHEN R =',F6.3,' LU, MORSE 
X POTENTIAL PARAMETERS ARE ' , L0A4 , // , 10X , ' CGD1 = ', 

X F8 « 4 , ' , CGD2 =',F8.4,', CGBI =',F8.4,', CGB2 = ', 

X F8.4 , 1 , CGFI = ' , F8 . 4 , ' , CC-F2 =',F8.4,/) 

973C FORMAT ( ' WHEN',F8.4,' <R<',F8.4,' THE MATCHING POTENTI 
XA L PARAMETER ARE',//,3X,' CPC =',£12.4,', CPI =', 

X E 12 . 4 , '8 CP2 = ' , E 12 . 4 , ' , CP3 = • ,E12 .4 , /7X , • CFO =', 

X E 12 .4 , ' , CF1 = ' , E 12 • 4 , ' , CF2 =',E12.4,//) 

974C FCPMAT ( 11 H THE FIRST ,I3,22H ATOMS ARE BLOCKED, ,13, 
X 24H LAYERS ARE FREE TC MOVE,/) 

WRITE ( 6,9610) IH1,IHS 

WRITE ( 6,9640) I C R , SUR , E V , I X , I Y , I Z , I PNU M 

WRITE ( 6,9645) I HTARG , IHBULL , C VR 

WRITE (6,9650) TM AS 1 , ABUN 1 , BMAS , I TEMP , T AMP 

WRITE (6,9652) TM AS2 , ABLN2 , OCR F 

WRITE ( 6,9660) I HT , E X A , EXB , FXA 

WRITE ( 6,9655) I HB, PEXA, PEXB, P FXA 

WRITE ( 6,9665) B X , B Y , B Z , CO X , C C Y , CC Z 

WRITE ( 6,9670) R X I ( 1 ) , RY I ( 1 ) , R Z I ( 1 ) 

IF( IPRCG.NE.4) GO TO 130 

WRITE ( 6,9730) R CEA , ROEB , C PO , C P 1 , CP2 , C P3 , C FO , CF1 , CF2 
WRITE(6,9725)R0EC,R0EB, IH2 , CGD 1 , CGC2 , CG3 1 , C GB2 , C GF 1 , 

X CGF2 

IOC WRITE ( 6,9675) NT, TIME, DT 
NLINE=1 
RETURN 
ENC 



SLBROLTINE KPNT 



THIS SUBROUTINE PRINTS THE HEADING 
FOR OUTPUT THAT IS FILED 



COMMON/ COM4 /IX, I Y,IZ,IXP,IYP,3CX,$CY,SCZ,ICEEP, 

X P EX , RB Z 

COMMON /COM 16 /BMAS , TMAS1 ,TMAS2 , A BUM , ABUN2 , I FNUM 
CCMMCN/C0M16A/BX, BY,BZ, EDISP,CVF 

CCMMCN/C0M17/TIME ,NT,DT , DTI, TEMP, CEL S,EV,NPAGc,NLINE 
CCMMCN/C0M20/IHK20) , IHS( 10) , I HT ( 6 ) , IHB ( 6 ) , IHTARG ( 2 ) , 
X IHBUL L ( 2 ) 

CCMM0N/C0M21/ IHB POT (2 ) , IHT POT ( 2 ) , SUR ( 2 ) , I H2 ( 10 ) , I OR , 

X L SUR 

96 4 C FORMAT ( ' (',13,') P LANE ' , 10X , ' I ON ENERGY =',F7.1,10X, 

X 'CRYSTAL SIZE (',12,' X',I2,' X ', 12 ,')',// ) 

96 A 5 FORMAT ( ' TARGET -',2A4,20X,' BULLET -',2A4,6X, 

X ' POTENTIAL • , 2 ( 2 A4 , IX ) , 2 A4 ) 

96 5 C FORMAT ( ' MASS 1 =' ,F7 . 2 , 2X ,' ABUNCANCE =',F5.3,4X, 

X 'MASS =' ,F7.2,16X,2(2A4,1X) ,2A4) 

9652 FCPMAT (• MASS 2 =', F7 .2 , 2X ,' ABUNCANCE =',F5.3,//) 
WRITE ( 6,9640) 1 0 R , E V , I X , I Y , I Z 

WRITE (6,9645) I H TARG, I hBUL L , I HT 
WRITE (6,9650) TM A S 1 , ABUN 1 , EMAS , I HB 
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WRITE (6,9652) TMAS2,ABUN2 

RETURN 

END 



SUBROUTINE HISTO( IN , FRE G1 , R ANGE , T ITL E 1 ) 

C 

C 

C THIS SUBROUTINE CRAWS THE FISTORAM (FREQUENCY 

C VS. EJECTEC ANGLE) FOR THE TYPE I HEAVY ATOMS 

C 
C 

D IMENSICN FREQ1 (20 ) , RANGE ( 20 ) , JOLT ( 20 ) , TITLE 1(24) 
CIMENSICN P CUT ( 20 ) 

CATA NOTH/' •/ t K/'****'/ 

2 FORMAT ( ' 0 INTERVAL ' ,2X, 19 ( F5 . 1 , l X ) , F5 . 1 ) 

4 FORMAT ( //,8X,24A4/) 

5 FCRMATt 'OFREQUENC Y ’ ,2016) 

6 FORMAT ( '0 INTERVAL ' ,2X, 19 ( F5 .1 , IX) ,F5. 1 ) 

7 FORMAT ( ' , 12 8 ( ' -' ) ) 

WRITE (6,4) TITLE1 

J N = I N 
I N T = 0 

IF(JN.GT.O) GO TO 10 

I NT = 1 

JM-JN 

10 CC 12 1=1, JN 
12 JCLT( I) =FREC1 ( I ) 

WRITE (6,5) ( J0UT(I),I=1,JN) 

WRITE (6,7) 

C FIND LARGEST FREQUENCY 

C 

FMA X = 0 . 

CC 20 1=1, JN 

IF(FREQ1 (I ) .GT.FMAX) FMAX = F REQ 1 ( I ) 

2C CONTINUE 
C 

C SCALE 

C 

JSCAL=1 

I F (FMAX .GT .60. ) J SCAL= ( FMAX+59 . )/60. 

F SCAL=J SCAL 
CC 50 1=1, JN 
5C JCLT( I ) =NOTH 
MAX=FMAX/FSCAL 
CC 8 J 1=1, MAX 
>=MAX-( 1-1) 

CC 70 J = 1 » J N 

I F ( FREQ1( J J/FSCAL .GE .X ) JOUT ( J ) =K 
7 C CONTINUE 
IX=X*FSCAL 

£0 WRITE (6,2) I X , ( J CUT ( J ) , J = 1 , J N ) 

2 FCFMAT ( I6,4X,20(2X,A4) ) 

WRITE (6,7) 

IF( INT. EQ. 1 ) GO TO 16 

WRITE (6,3) (RANGE ( J ) , J = 1 , J N ) 

C*C TO 15 
16 DC 51 1=1, JN 
51 FCLT (I ) =RANC-E( I ) 

WRITE (6,6) (POUT (I ) ,1=1, JN) 

15 CONTINUE 

WRITE (6,100) 

ICC FORMAT! IX, 'MIDPOINT'/) 

RETURN 

ENC 
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SUBROUTINE FIS TO ( IN, FREC2, RANGE, TITLE2) 



C 

c 

c 

c 

c 

c 



c 

c 

c 



c 

c 

c 



2 

3 

4 

c 

6 

7 



1C 

12 



2C 



5 C- 



7 C 

8 0 



16 

51 

15 

ICO 



THIS SUBROUTINE CRAWS THE FISTORAM (FREQUENCY 
VS. EJECTEC ANGLE) FOR THE TYPE2 LIGHT ATOMS 



DIMENSION FREQ2 ( 2 J ) , R ANGE ( 2 0 ) , JOUT ( 20 ) , T ITL E2 ( 24 ) 
DIMENSION POUT (20 ) 

CATA NOTH/' ' / , K/ ' **** ' / 

F C F MAT ( I6,4X,2J(2X,A<,) ) 

FORMAT ( '0 INTERVAL ' ,2X , 19 ( F5 .1 , IX ) , F5. 1 ) 
FCFMAT(//,8X,24A4/) 

FORMAT ( ' OFREQUENC Y ' ,2016) 

FORMAT ( '0 INTERVAL • , 2X , 1 9 ( F5 . 1 , 1 X ) , F5 . 1 > 

FORMAT ( ' - ' , 128 ( ' - ' ) ) 

WRITE (6,4) TITLE2 
J N = I N 
I N T = 0 

IF(JN.GT.O) GO TO 10 

I N T = 1 

JN=-JN 

OC 12 1=1, JN 
JOLT ( I ) =FREC2 ( I ) 

WRITE (6,5) ( JOUT (I),I = 1,JN) 

WRITE (6,7) 

FINC LARGEST FREQUENCY 

F M A X = 0 . 

CC 20 1=1, JN 

I F ( F REQ 2 ( I ) .GT.FMAX) FMAX =F REQ2 ( I ) 

CONTINUE 

SCALE 

JSC AL=1 

IF(FMAX.GT.60.) J SCAL= ( FMAX+59 . )/60. 

F SCAL = J SC A L 
CC 50 1=1, JN 
J CUT ( I ) =NOTF 
MAX=FMAX/FSCAL 
CC 80 1=1, MAX 
X=MAX-< 1-1 ) 

CC 70 J = 1 , J N 

IF ( FREQ2 ( J )/FSCAL .GE .X ) JOUT ( J )=K 

CONTINUE 

IX = X 5 ! ! FSCAL 

WRITE (6,2) IX, (JCUT ( J) ,J=1 ,JN) 

WRITE (6,7) 

I F ( INT. EQ.l ) GO TO 16 

WRITE (6,3) ( RANGE ( J ) ,J = 1 » J M 

GC TO 15 

CC 51 1=1, JN 

PCLT (I ) =R ANGE ( I ) 

WRITE (6,6) ( POUT ( I ) , 1 = 1 , JN ) 

CONTINUE 
WRITE (6, ICC) 

FORMAT! IX, 'MIDPOINT'/) 

RETURN 

END 



SLEROUTINc KFOLDA 
C 
C 

C THIS SUBROUTINE STATISTICALLY FOLDS THE 

C TARGET IMPACT AREA FOR THE FCC (100) CRYSTAL 

C LATTICE TO REPRESENT FULL COVERAGE 

C 
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CCMMON/CCM1/RX(600) ,RY(60 3) ,RZ ( 60 0 ) ,LCUT (60 C ) , LL , LD 
COMMON/ COM4 /IX, IY,IZ,IXP,IYP,SCX,SCY,SCZ,IDEEP, 

X REX R B Z 

C CM MON/ COM4 A/ LNUM6J0) ,LNUMI(600) ,LNUMF(6GC) 

CC M MON/ COM2 1/ IH8PCT (2 ) » IHT PCT ( 2 ) , SUR ( 2 ) , I F2 ( 1 0 ) , I OR , 
X LSLR 
I LA YER= 1 
N F IRST = 2 



DETERMINE THE FIRST ATOM AND THE LAFGEST 

NUMBERED ATOM IN EACH LAYER 

10 GC TO (20, 30, 40, 50, 6C, 70) , ILAYER 

20 MAY =( (LL-l)/< IY/2) )/2+2 
K =NFIRST 
NLAY 1=N LAY 
GC TO 80 

30 NLAY=( ( LL-1 )/( IY/2 ) ) +1 
K=NF IRST +IX/2 
GC TO 80 

40 MAY = NIAYI + NLAY -I 
K = NF I RS T 
GC TO 80 

5C NLAY = NLAY1 + NLAY -2 
K=NFIRST +IX/2 
GC TO 80 

6C NLAY = NLAYI + NLAY -3 
MNFIRST 
GC TO 80 

7 C NLAY = NLAYI + NLAY -4 
K =N F I RS T +IX/2 
GC TO 8 C 

FIRST (DIAGCNAL) FOLD - ALL LAYERS 

80 DC SO J = NF I RST , NL A Y 

LNUM ( J ) =LNUN ( J) + LNUM(K) 

IF(J.EQ.K) LNUM( J) =2*LNUM( J ) 

LNUM ( K ) =LNUM ( J ) 

K =(< + I X 

IMK.GT .NLAY) K=K-NL AY 1 + 2 

SC CONTINUE 

DC 100 J=NF IRST , N LAY 
LNUM. ( J ) =LNUM(J ) / 2 
IOC CONTINUE 

SECOND ( HCR I ZONTAL ) FOLD - ODD LAYERS 



105 I COUNT = 1 
I S ET = 0 

GC TO ( 110,160, 110, IdO, 110,160) , ILAYER 
110 K=NLAY-IX/2 

CC 150 J=N F IRST , N LAY 
LNUM( J)=LNUM( J)+LNUM(K) 

IF(J.EQ.K) LNUM( J )=2*LNUM(J ) 

LNLM ( K ) =LNUM ( J ) 

ICCUNT= ICQLNT+1 
IF(ISET.EQ.l) GO TO 130 
I F ( ICCUNT.GT. IX/2 +1) GO TO 120 
K=K + 1 
GC TO 150 
1 2 C I SET = 1 
K =K- 1 X+ 1 
I C CUNT = 1 
GC TO 150 

13C IF (ICCLNT.C-T. IX/2) GO TO 140 
K=K + 1 
GG TO 150 
1 4 C I 5 ET = 0 
K = K- 1 X+ 1 
I CCLNT= 1 
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1 5 C CCMINUE 

C-C TC 210 

SECCNO (HORIZONTAL) FOL C 

16C K = N LAY- IX/2+1 

DO 200 J = NF IRST » NL AY 
LNUM J )=LNUI*( J) +LNUM (K ) 
IF(J.EQ.K) LNUM(J )=2 # LNUM (J ) 
LNLM( K) =LNUN( J ) 

I CCUNT = IC0UNT+ 1 
IF( ISET.EQ.l) GO TO 180 
IF ( ICGUNT.GT.IX/2 ) GO TO 170 
K = K + l 
GC TO 2C0 
170 I S ET= 1 
K =K- I X+ 1 
I CCLNT= 1 
C-C TC 2C0 

180 IF ( ICCUNT.GT. IX/2 + 1) GC TO 190 
K=K + 1 
C-C TO 2 CO 
190 I«ET=0 
K = K- 1 X+ 1 
I CCUNT= 1 

2 0 C CONTINUE 

2 1C DC 220 J = NF IRST ,NLAY 
LNLfM J )=LNUV(J) /2 
2 2 C CONTINUE 

THIRD (VERTICAL) FOLD - 

I CCUNT= 1 ' 

I SET = 0 

GC TO ( 23C, 280,230, 220,230, 280) , 
220 K =N F IR S T + I X/2 

DC 270 J=NF IRST,NLAY 
LNLM J)=LNUM( J)+LNUM(K) 
IF(J.EQ.K) LNUM(J )=2*LNUM (J ) 
LMH( K) =LNUN( J ) 

I COUNT = ICCUNT+ 1 
I F ( ISET.EQ.l) GO TC 250 
IF(ICGUNT.GT.IX/2 +1) GO TO 240 
K = K- 1 
GC TO 270 
240 I SET= 1 
K = K+ IX- 1 
I C CLNT= 1 
GC TO 270 

250 IF (ICGUNT.GT.IX/2) GO TO 260 
K = K- 1 
GC TO 270 
26C I S ET=0 
K = K + 1 X- 1 
I C CUNT= 1 
27C CONTINUE 
GC TO 330 

THIRD (VERTICAL) FOLC - 

280 K=NFIRST+IX/2-l 

DC 220 J= NF IRST T NLAY 
LNU/ I MJ)=LNUM(J) + LNUM(K) 
IF(J.EQ.K) LNUM( J)=2-LNUM( J ) 
LNLM( K)=LNUN( J ) 

I CCLNT = ICCUNT + 1 

I F ( ISET.EQ.l) GO TC 300 

IF ( ICOUNT .GT.IX/2 ) GO TO 290 



250 ISET=1 
K=K+IX-1 



- EVEN LAYERS 



GDC LAYERS 
I LAYER 



EVEN LAYERS 
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I COUNT = 1 
GC TO 320 

30C I F ( ICGUNT . GT. IX/2+1) GO TO 310 
K = K- 1 
GC TO 320 
3 1C I SET = 0 
K=K4IX-1 
I CCUNT= 1 
32C CONTINUE 

33C CC 340 J=NF IRST » N L AY 
LNLM( J) =LNUM( J ) /2 
34C CONTINUE 

I LAYER= ILAYER+1 
NFIRST=NLAY+1 

IF ( ILAYER.C-T. IY) GC TO 350 
GC TO 10 
35C RETURN 
5NC 



SLEROUTINE KF0L03 



THIS SUBROUTINE STATISTICALLY FOLCS THE 
TARGET IMPACT AREA FOR THE FCC (110) CRYSTAL 
LATTICE TC REPRESENT FULL COVERAGE 



CC^MCN/CGM1/RX(60C) ,RY(600),RZ(630),LCUT(60C),LL,LC 
CCMMCN/C0M4/IX, IY,IZ,IXP,IYP,SCX,SCY,SCZ,IDEEP, 

X REX , R8 Z 

CCNMCN/C0M4 A/ LNUPM60 0) ,LNU M I ( 6 00 ) , LNUM F (60 C ) 
CCMMCN/CGM21/IH6PGT (2 ) » IHTPOT(2 ) ,SUR(2 ) , IF2 (10 ) ,I0R , 
X LSLR 
I LA YER = 1 
N F I RST = 2 

DETERMINE THE FIRST ATOM AND THE LARGEST 
NUMBEREO ATOM IN EACH LAYER 

1C GC TO (20,3C,40,3C,40,3C) tILAYER 
2C NLAY=( ( IX+1 )/2 ) *( ( IZ + 1 )/2) + l 
N L A Y 1 =N L AY 
GC TO 50 

3C N LAY = NLAY + ( IX/2 )* ( IZ/2 ) 

GC TC 50 

AC NLAY=NLAY+NLAY1-1 
5 C I COUNT = 1 

GC TO (6C,50,60,90,50,5C),ILAY£R 

SECOND (HORIZONTAL) FOLD - OOD LAYERS 

60 K = N LAY- I X/2 

CC 80 J = NF I RST » NL AY 
LNLM( J)=LNUM(J)+LNUM(K) 

IF(J.EQ.K) LNUM ( J ) = 2*LNUM(J) 

LNLM ( K) =LNUM( J ) 

I COUNT = ICGUNT+1 

IF (ICCUNT.GT. IX/2 +1) GO TO 70 
K =K+ 1 
GC TO 3 C 
70 K=i< — IX 
I COUNT = 1 
80 CCNTINUE 
GC TO 120 

SECCND (HCRIZONTAL) FOLC - EVEN LAYERS 

SC K=NLAY-lX/2 +1 

CC 110 J=NFIRST, NLAY 
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LNUM( J)=LNUM( J) +LNUM(K ) 

IF(J.EQ.K) LNUM ( J ) = 2*LNUM(J) 

LNLM( K) =LNUM( J ) 

I CCUNT= ICOUM + 1 
IF(ICQUNT .GT.IX/2) GO TO 100 
K =K + 1 
C-C TO 110 
ICC K=K-IX+2 
I CCUNT = 1 
11C CONTINUE 

1 2 C DC 130 J=NFIRST,NLAY 
L Mi N ( J ) = LNUM ( J ) / 2 
120 CONTINUE 
I CCUNT = 1 

GC TO ( 140, 190, 140, 190, 140, 190) , ILAYER 

THIRD (VERTICAL) FOLD - CCD LAYERS 

1 AC K=NFIRST+IX/2 

CC 160 J = NF IRST , NLAY 
LNUM(J)=LNUM(J)+LNUM(K) 

IF(J.EQ.K) LNUM ( J ) = 2*LNUM(J) 
LNLM(K)=LNUM(J) 

ICCUNT=ICCUNT+1 

IF (ICQUNT.GT. IX/2+1) GO TO 150 
K=K-1 
GC TO 160 
1 5 C K = K+ I X 
I C CLNT = 1 
1 6 C CONTINUE 
GC TO 220 

THIRD (VERTICAL) FOLD - EVEN LAYERS 

150 K=NFIRST +IX/2 -1 

CC 210 J = NF IRST , N LAY 
LNLM(J)=LNUM(J)+LNUM(K) 

IF(J.EQ.K) LNUM ( J ) = 2*LNUM(J) 

LNLM( K) =LNUM J) 

I CCUNT= ICCUNT+1 
IF ( ICCUNT.GT .IX/2) GO TO 200 
K=K-1 
C-C TO 210 
2CC K =K + I X-2 
I CCUNT = 1 
210 CCNTINUE 

22C CC 230 J=NF IRST, NLAY 
LNLM(J)=LNUM(J)/2 
22C CCNTINUE 

I LA YE R= ILAYER + 1 
N F IRST = NLAY + 1 

I F < ILAYER. GT.IY) GO TO 240 
GC TO 10 
2 4 C RETURN 
ENC 



SLEROUTINE KFOLDC 



THIS SUBROUTINE STATISTICALLY FOLDS THE 
TARGET IMPACT AREA FOR THE FCC (111) CRYSTAL 
LATTICE TC REPRESENT FULL COVERAGE 
*** VALID ONLY FOR THE FIRST LAYER AND IX £ IY =5 



COMMON/ CCM1/RX( 600 ) ,RY(600) ,RZ(6CC),LCUT(6CC),LL,LD 
COMMON/ COM4 /IX, I Y ,IZ , IXP,IYP,SCX,SCY,SCZ,IDEEP, 

X R EX , RE Z 

COMMON/ COM4 A/ LNUM (600) , LNUM I (600 ) ,LNUMF (60 0 
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CCMNCN/C0M21/IH8PCT ( 2 ) , IHTPCK2 ) , SUR ( 2 ) , IH2 { 10 ) , I CR , 
X LSLR 
CC l J = 2 » L L 
LMJNF ( J ) = 0 
LNLM I ( J )=LNUM( J ) 

1 CONTINUE 
I CL T=0 

I T I R5 = 1 

NLAY=( < LL-2 ) / ( IY/2 ) ) / 2 + 2 

2 N F IFST = 2 + I X/2 
I CCUNT=C 
ISET=0 

K =N LA Y- I X/2 
J=NFIRST 
5 CCMINUE 



FIRST ( + SLOPE DIAGONAL) FOLD - FIRST LAYER 



L N U M ( J)=LNUN(J)+LNUN(K) 

IF (J.EG.K) LNUM ( J ) = 2*LNUM( J ) 
LNLM( K)=LNUP( J) 
ICCLNT=ICCLNT+1 
IF (ISET.EQ.O .AND. 



IF (ISET.EQ.O .ANC. 

IF (ISET.EQ.O .AND. 

IF (ISET.EQ.2 .ANC. 

IF (ISET.EQ.2 .ANC. 

IF (ISET.EQ.2 .AND. 

IF (ISET.EQ.2 .ANC. 

1 F ( IC OUNT . EG .4 .OR. 

X .CR. ICOUNT.EQ .17 ) GC TO 25 
IF ( ICOUNT.EQ. 20 ) GO TC 30 
J=J + 1 
K=K + 5 
GC TC 5 

1C J=J+( IX-3 ) / 2 
K = K- 1 X 
GC TO 5 
15 J=J+1 

K = K + ( IX + 1 )/2 
GC TO 5 

2C ISET=ISET+1 
I CCLNT = C 
22 J = J + 2 
K = K- 14 
GC TO 5 
25 K = K- 1 9 
J=J + 1 
C-C TO 5 

3C ISET=ISET+1 
I C CLNT = 0 

J=J + 1 

K =K + 5 
GC TO 5 

35 IF (I0UT.5G.1) GO TC 150 
GC TO (50,50,100) , I TI ME 



ICOUNT.EQ. 1) GO TO 10 
ICOUNT.EQ. 2) GO TO 15 
ICOUNT.EQ. 3) GO TO 20 
ICOUNT.EQ. 1) GO TO 22 
ICOUNT.EQ. 2) GO TO 15 
ICOUNT.EQ. 3) GO TO 10 
ICGUNT.EQ.4) GO TO 35 
ICGUNT.EQ.8 .OR. ICCUNT .EG . 13 



SECOND (- SLOPE DIAGONAL) - FIRST LAYER 



5C CONTINUE 
N F I RST= 2 
I F LAG = 0 
ICCLNT=C 
I S ET = 0 
I CL T = 0 
K =N LA Y 
J=NFIRST 
55 CONTINUE 

LNLM( J )=LNUM( J)+LNUM(K ) 

IF (J.EG.K) LNUM ( J ) =2*LNUM ( J ) 
L NLM ( K ) =L NUN ( J ) 

ICCUNT= ICOUNT+1 
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60 



65 



7C 

72 



75 



77 
8 C 

85 
$ c 



ICC 



lie 

1 2 C 

120 

140 

15C 



IF (ISET.EQ.O .ANC. 

IF (ISET.EQ.O .ANC. 

IF (ISET.EQ.O .AND. 

IF (ISET.EQ.2 .AND. 

IF (ISET.EQ.2 .ANC. 

IF (ISET.EQ.2 .ANC. 

IF (ISET.EQ.2 .ANC. 

IF (ICCUNT.EQ.4 .OR. 

IF (ICGINT.EQ.8 .OR, 

IF ( ICOUNT. EQ. 20 ) GO 
J = J + 1 
K = K-4 
C-C TC 5 5 
J = J + 5 
K=K-5 
GC TO 55 
j=J+l 
K = K-4 
GC TO 55 
I SET=ISET+1 
I C CUNT = 0 
J = J + 3 
K = K+3 
GC TO 55 
J = J + 2 
K = K +7 
GC TO 55 
J = J+1 
K =K + 1 1 
GC TO 55 
I 'ET=ISET+1 
I CCUNT = 0 
J=J + 1 
K = K-4 
GC TO 55 

GC TO (100,95,150) ,ITIME 

ICLT=1 

C-C TO 2 



ICOUNT. EC. 1) GO TC 6C 
ICOUNT. 6Q. 2) GO TO 65 
ICOUNT. EQ. 3) GO TC 7C 
ICOUNT. EQ.l) GO TC 72 
ICOUNT. EQ. 2) GO TO 65 
ICOUNT. EG. 3) GO TC 6C 
ICOUNT. EQ. 4) GO TO 85 
ICOUNT. EQ. 17) GO TC 75 
ICOUNT. EC. 13) GO TC 77 
TO 80 



THIRD (VERTICAL) FOLD - FIRST LAYER 



CONTINUE 

I CCUNT = 1 

I S ET = 0 

NFIRST=2 

K =NF IRST+IX/2 

DC 140 J=NFIRST,NLAY 

LNLV ( J ) =LNUM( J)+LNUM(K) 

IF (J.EQ.K) LNUM{ J)=2*LNUM( J ) 
LNUNt K ) =L NUN ( J ) 

I C CUNT= ICCUNT + 1 

IF (ISET.EG.l) GO TC 120 

IF (ICOUNT. GT. IX/2 +1) GO TO 110 

K =K - 1 

GC TG 140 

I SET = 1 

K = K + 1 X- 1 

I C CLNT = 1 

C-C TO 140 

IF (ICCLNT.GT. IX/2) GO TO 130 

K =K- 1 

GC TO 140 

I S ET = 0 

K =K + I X- 1 

I CCUNT = 1 

CONTINUE 

GC TO ( 150,150,50) , I T I ME 
IT INE= IT I ME+1 
I C LT = 0 

I CCUNT =0 

DC 152 J = 2 , N LAY 
LNLM( J)=LNUN(J)/8 
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LMHF ( J)=LN0MF( J ) +LNUM ( J ) 
152 CONTINUE 

IF ( I T IM E . GT.3 ) GC TO 160 
CC 155 J = 2 , MAY 
LNUM(J)=LNUMI(J) 

155 CONTINUE 
GC TO 2 

16C CC 170 J = 2 , NLAY 

LNUM( J)=LNUMF( J)/4 
1 7 C CONTINUE 
K =N LA Y+ 1 
CC 180 J = K , LL 
LNUM( J) =6*LNUM( J ) 

18C CONTINUE 
RETURN 
END 



SUBROUTINE K PROS 



KPRCB IS DESIGNEC TO STUDY THE PROBABILITY 
THAT A RING OF ATOMS (AT A GIVEN DISTANCE FROM 
THE LAYER CENTER) WILL SPUTTER 



COMMON/ C0M1/RX(60C) ,RY ( 600 ) f RZ ( 60 C ) , LOUT ( 60 C ) t LL , LD 
CCMMCN/C0M4/IX, IY , IZt IX F, IY P ,SCX , SCY ,SCZ , I DEEP, 

X F EX , RE 2 

CCMMCN/C0M4A/ LNUM(600) , LNU M I ( 6 CC ) ,LNUM F ( 60 C ) 

COMMON/ COM! 1/RXI (600),RYI(oC0) ,RZI(600) 

CCMM0N/C0M21/ Ih BP CT(2 ) , IHTPCT ( 2 ) , SUR ( 2 ) ,IH2(10),ICR, 

X L SUR 

CCMMC N/ COM2 4/ DNN(IOO) ,KAT0M(100) ,KSPUT ( 100 ) , PROB ( ICO ) 
COMMON/ COM2 4 A/ KSUM , K R I NG , R CIST ( 600 ) 

CC 2 1=1,100 
K ATOM ( I ) =0 
KSFUT ( I ) =0 
FRCE (I ) =0.0 
2 CONTINUE 

DETERMINE THE DISTANCE FROM EACH ATOM TO THE 
INITIAL (0,0) POINT OF TARGET IMPACT AREA 

CC 5 1=2, LL 
XC I ST =R X I ( I )-RBX*SCX 
ZD 1ST =R Z I { I )-RBZ*SCZ 
SCIST=X0IST**2 +ZCIST**2 
R C 1ST ( I ) =SD 1ST** . 5 
5 CONTINUE 
K F I RST= 1 
I L AYER= 1 
NFIPST=2 

IF (I0R.EQ.110) GO TO 170 
IF (IOR.EQ.lll) GO TO 200 

DETERMINE THE FIRST ATOM AND THE LARGEST 
NUMBERED ATOM IN EACH LAYER FOR 100 LATTICE 

10 GC TO ( 20,30,40, 5C ,60,70) , I LAYER 
20 NLAY =( (LL-1) /( IY/2) )/2+2 
N L AY1 =N LAY 
C-0 TO 8C 

3C NLAY=((LL-l)/(IY/2) ) + l 
GC TO 80 

4C NLAY=NLAY1 + NLAY -1 
C-C TO 8 0 

50 NLAY = NLAY1 + NLAY -2 
GO TO 8 C 

60 NLAY = NLAY 1 + NLAY -3 
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C-C TO 8 C 

7 C NLAY = NLAY1 + NL AY -4 
C-C TO 80 
C 

C SORT DISTANCES IN A LAYER INTO ASCENCING CRDER 

C 

8C N1=NLAY~1 

DC 100 I=NFIRST,N1 
N F 1= I + 1 

CC SO J = NF1 jNLAY 

IF (RCIST(J).GE.RDIST(I)) GO TO 50 
SA VE=RD I ST ( I) 

RD 1ST ( I J =R C 1ST < J ) 

RC 1ST ( J ) = S AVE 
SC CONTINUE 
ICC CONTINUE 
C 

C ESTABLISH A RING AT EACH CISTANCE (+ OR - DR) 

C 

K=KFIRST 

CC 120 I = NF IRST ,NLAY 
J = I + 1 

IF ( J.GT .NLAY) GO TO 110 

IF <RDI$T( J)-RDIST(I) .GT. 0.001) GC TO 110 
C-C TO 120 

11C CNN(K)=RDIST(I) 

K = K + 1 

120 CONTINUE 
C 

C REINITIALIZE THE DISTANCE FROM EACH ATOM 

C TO THE REFERENCE (0,0) POINT OF IMPACT AREA 

C 

CC 145 1=2, LL 
XC IST=RXI ( I )-RBX*$CX 
ZC IST = R Z I < I )-RBZ*SCZ 
SCIST=XDIST**2 + ZCIST**2 
R C I ST ( I ) = SC 1ST** . 5 
145 CONTINUE 
C 

C DETERMINE THE NUMBER OF ATOMS AND SFUTS FPCM 

C ATOMS IN EACH RING 

C 

K R ING=K -1 

CC 150 J = KF IRST , K R ING 
CC 140 I = NF IRST, NLAY 

IF<ABS<RDIST< I)-DNN(J) ) .LT. 0.001) C-0 TO 130 
GC TO 140 

130 KATCN ( J ) =K A TOM ( J ) +1 

KSPUT ( J ) =KSPUT ( J ) +LNUM ( I) 

14C CONTINUE 
C 

C DETERMINE THE PROBABILITY OF AN ATOM IN EACH 

C RING SPUTTERING 

C 

IF (KSUM .EQ. 0) KSUM = 1 

IF (I0R.EQ.100) PROB ( J ) =FLOAT < KSPL'T ( J ) )/ < FL CAT ( KSUM ) * 

X £ .0 ) 

IF <I0R.EQ.110) PROB (J)=FLQAT (KSPUTIJ) ) / (FLCAT(KSUM)* 

X 4.0) 

IF ( I OR . EQ . Ill) PROB (J)=FLOAT(KSPUT( J) ) /(FLCAT(KSLM)* 

X 6.0 
15C CCNTINUE 

ILAYER=I LAY ER + 1 
N F IRST = NLA Y + l 
K F IRST = K 

155 IF ( ILAYER.GT.IY) C-0 TO 160 
IF <I0R.EQ.100) GC TO 10 
IF < I OR . EG . 11 0 ) GC TO 170 
IF (I OR. EC. Ill) GC TO 200 
16C RETURN 
C 

C DETERMINE THE FIRST ATOM ANC THE LARGEST 
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C NUMBERED ATOM IN EACH LAVER FOR 110 LATTICE 

C 

17C GC TO ( 180, 190, 195, 190, 195, 190) , ILAVER 
ISC NLAY=( ( IX + 1 )/2)*( (IZ + l)/2) + l 
M A Y1=NLAY 
GC TO 8C 

190 N L A Y = NLA Y + ( I X/2 ) * < IZ /2 ) 

GC TG 80 

195 N LAY = NLA Y + < ( IX+1 ) /2 ) * ( ( IZ+ 1 ) / 2 ) 

GC TO ac 

DETERMINE THE FIRST ATOM AND THE LARGEST 
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2CG GC TO (210,220,230,240) ,ILAYER 
2 1C NLAY=( ( LL-2 )/( IY/2) )/2 + 2 
NLAY1=NLAY 
GC TO 8 G 

22C NLAY = ( ( LL-2 )/( IY/2) ) + l 
r r t n ft r 

22 C NLAY=NLAY1 + NLAY-1 
GC TO 8C 

24C NLAY=NLAY1+NLAY-1 
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Figure 1 - Sputtering Summary - (100) 2keV Cu/Ne+ KSE-B 




Figure 2 - Sputtering Summary - (100) 2 keV Cu/Au + KSE-3 
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Figure 3 - 9x4x9, (100) Orientation .iicrocrystall 
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igure 4 - 9x4x9, (110) Orientation Microcrystallite 
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Figure 5 - 9x4x9, (111) Orientation Microcrystal li 4 -? 
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Figure 6 - Microcrystallite Boundaries 




Figure 7 - intrinsic Geometry of Target Impact Areas 
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Seduced Representative Impact Areas 
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Figure 9 - Expanded Impact Area - (100) Orientation 
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Figure 10 - Expanded Impact Area - (110) Orientation 
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Figure 11 Expanded Impact Area — (111) Orientation 



(10 0 ) 




139 




Figure 13 
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Orientation Folding Geometry 
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Ficrure 14 - Spattering Ratio vs. Ion Atomic Number - 

(100) , 100-eV, 9x4x9 
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Figure 14 - Sputtering Ratio vs. Ion Atomic Number - (100), 100-eV, 9x4x9 
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Figure 15 - Sputtering Summary - (100), 100-97, Cu/H a+ 
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Figure 16 - Ring Probability - (100) , 100-eV, Cu/'i 
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Figure 17 - Angular Dispersion - (100) , 100-eV, Cu/He + 
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(a) Gibson-II 




Figure 18 - Sputtering Summary - (100), 100-eV, Cu/Ne + 
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Figure 19 - Ring Probability - ( 100) , 100-eV, Cu/Ne + 
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Figure 20 - Angular Dispersion - (100 ) , 100-eV, Cu/Ne + 
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(a) Gibson-II 




(b)KSE-B 



Figure 21 - Sputtering Summary - (100) , 100-eV, Cu/Ar + 
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Figure 22 - Ring Probability - ( 100) , 100-e'/, Cu/Ar + 
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Figure 23 - Angular Dispersion - (100) , 100-?V, Cu/Ar + 

Gibson- II 
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Figure 24 - Angular Dispersion - (100) , 100-e7, Cu/Ar + 

KSE-3 
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Figure 25 - Sputtering Mechanism Trace - 
(1 00), 100-eV, Cu/Ar + , Gibson-II 
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(a) Gibson-II 




(b)KSE-B 



Figure 26 - Sputtering Summary - (100), 100-eV, Cu/Cu + 
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Figure 27 - Ring Probability - (100) , 100-eV, Cu/Cu + 
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Figure 28 - Angular Dispersion - (100) , 100-e7, Cu/Cu + 
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(a) Gibson-II 
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Figure 29 - Sputtering Summary - (100) , 100-eV, Cu/K 
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Figure 30 - Ring Probability - ( 100) , 100-eV, Cu/K 
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Figure 31 - Angular Dispersion - (100), 100-eV, Cu/Kr + 
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igurs 32 - Sputtering Summary - (100) , 100-eV, CU/X5+- 
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Fiaare 33 - Ring Probability - ( 100) , 100-e7, Cu/Xe + 
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Figure 34 - Angular Dispersion - (100) , 100-eV, Cu/X 
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Figura 35 — Sputtering Summary (100) , 100 eV, Cu/Au + 
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Figure 36 - Sing Probability - ( 100) , 100-eV, Cu/Au + 
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Figure 37 - Angular Dispersion - (100) , 100-?V, Cu/Au + 
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Figure 38 - Sputtering Summary - (100), 100-eV, Cu/Hg + 
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Figure 39 - Ring Probability - (100) , 100-eV, Cu/Hg + 
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Figure 40 - Angular Dispersion - (100), 100-eV, Cu/Hg* 

G ibso n-II 
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Figure 41 - Angular Dispersion - (100), 100-aV, Cu/Kg + 

KSE-3 
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Figure 42 - Sputtering Ratio vs. Ion Atomic Number - (100), 1-keV, 9x4x9 
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Figure 43 - Sputtering Summary - (100) , 1-keV, Cu/He + 
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Figure 44 - Ring Probability - (100) , 1-keV, Cu/He + 
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Figure 45 - Angular Dispersion - (100), 1-keV, Cu/He' 1 ' 
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Figure 46 - Sputtering Summary - (100), 1-kaV, Cu/Ne + 
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Figure 47 - Ring Probability - (100), 1-keV, Cu/Ne + 
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Figure 48 - Angular Dispersion - (100) , 1-keV, Cu/Ne- + 
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Figure 49 - Sputtering Summary - (100), 1-keV, Cu/Ar + 
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figure 50 - Ring Probability - (luO) , 1-ksV, Cu/Ar + 
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Figure 51 - Angular Dispersion - (100), 1-keV, Cu/Ar + 
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Figure 52 - Sputtecing Summary - (100) , 1-keV, Cu/Cu + 
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Figure 53 - Ring Probability - (100) , 1-keV, Cu/Cu + 
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Figure 54 - Angular Dispersion - (100), 1-keV, Cu/Cu + 
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' Figure 55 - Sputtering Summary - (100), 1-xeV, Cu/Kr + 
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Figure 56 - Ring Probability - (100) , 1-keV, Cu/Xr + 
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Figure 57 - Angular Dispersion - (100) , 1-keV, Cu/Kr + 
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Figure 58 - Spattering Summary - (100) , 1-keV, Cu/Xe* 
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Figure 59 - Ring Probability - (100) , 1-keV, Cu/Xe + 
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Figure 60 - Angular Dispersion - (100) , 1-keV, Cu/Xe 1 " 
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Figure 61 - Sputtering Summary - (100), 1-ksV, Cu/Au + 
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Figure 62 - Ring Probaoility - (100 ) , 1-keV, Cu/Au + 
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Figure 63 - Angular Dispersion - (100) , 1-k?V, Cu/Au + 
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Figure 64 - Sputtering Summary - (100) , 1-k.eV, Cu/Ha + 
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Figure 65 - Ring Probability - (100), 1-keV, Cu/Hg + 
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Figure 66 - Angular Dispersion - (100), 1-keV, Cu/H: + 
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(a) Case (1) and Case (3) 




(b) Case (2) and Case (4) 



Figure 68 - Sputtering Summary - (110) , 100-eV, Cu/hr + 

G ibson-IT 
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Figure 69 - Ring Probabilty - (110) , 100-eV, Cu/Ar + 

Gibson-II 
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Figure 70 - Angular Dispersion - (110) , 100-eV, Cu/Ar + 

Gibso n-II 
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Figure 71 - Angular Dispersion - (110), 100-eV, Cu/Ar + 

Gibso n-II 
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(a) Case (1) and Case (3) 
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Figure 73 - Ring Probability - (110) , 100-eV, Cu/Ar + 

KSZ-3 
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Figure 74 - Angular Dispersion - (110) , 100-eV, Cu/Ar + 

KSE-B 



201 



Ejection Angle (Degrees) 

I 1 1 1 1 ! 1 1 1 i 1 1 1 1 1 1 1 1 1 

0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 



Heavy 



I 1 1 — I 1 1 1 1 1 1 1 h— H 1 1 1 1 1 1 



Light 



I i 1 1 1 1 1 1 1 1 ! 1 1 1 1 1 1 1 1 

(a) Case (2) 



Heavy 



I 1 1 1 1 1 1 1 1 1 1 1 i 1 1 1 1 1 1 



Light 



l 1 1— I 1 1 i 1 1 ! i 1 i 1 1 1 i 1 1 

(b) Case (4) 



Figure 75 - Angular Dispersion - (110), 100-eV, Cu/Ar + 

KSE-B 
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Figure 77 - Ring Probability - (110) , 100-eV, Cu/Hg+ 

KS2-S 
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Figure 78 - Angular Dispersion - (110), 100-eV, Cu/Ha+ 



205 




(a) Gibson-II 




(b) KSE-B 



'ire 7 9 - Sputtering Summary - (111), Cu/A 
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Figure 80 - Atom Probability - (111) , 100-eV, Cu/Ar* 
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Figure 8 1 - Angular Dispersion - (111) , 100-eV, Cu/Ar + 

Gibson-II 
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Figure 82 - Angular Dispersion - (111), 100-eV, Cu/Ar + 

KSE-B 
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Figure 83 - Angular Dispersion - (111), 100-eV, Cu/Ar* 

KSE-3 
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(a) Case (1), Case (2), and Case (4) 




Figure SU - Sputtering Summary - ( 111 ) , 100 - 27 , 
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Figure 85 - Sputtering Summary - (111) , 100-e’ 
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Figure 86 - Atom Probability - (111), 100-eV, Cu/Hg+ 
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Figure 87 - Angular Dispersion - (111), 100-eV, Cu/Hg + 

Gibso n-IT 
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Figure 88 - Angular Dispersion - (111) , 100-eV, Cu/Hg + 

KSE-B 
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Figure 89 - Angular Dispersion - (111), 100-eV, Cu/Hg + 

K5E-B 
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(b) KSE-B 



Figure 90 - Sputtering Summary - (110) , 100-eV, Cu/Hg+ 
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Figure 91 - Sputtering Summary - (111), 100-aV, Cu/Hg + 

13x4x13 
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Figure 92 - Sputtering Ratio vs. Ion Atomic Number 
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Figure 93 



Spattering Ratio vs. Lattice Size 
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Figure 95 - Atom Probability - (100) , 100-aV, Cu/Cu + , 
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Figure 96 - Angular Dispersion - (100), 100-eV, Cu/Cu + 
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Figure 97 - Sputtering Summary - (100) , 200-eV, Cu/Cu 
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Figure 98 - Atom Probability - (100) , 200-aV, Cu/Cu + 
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Figure 99 - Angular Dispersion - (100) , 200-eV, Cu/Cu + 
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Figure 100 - Sputtering Summary - (100) , 500-eV, Cu/Cu 
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Figure 101 - Atom Probability - (100), 500-eV, Cu/Cu 1 - 
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Figure 102 - Angular Dispersion - (100), 500-eV, Cu/Cu + 
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Figure 103 - Sputtering Summary - (100), 1-keV, Cu/Cu 



I SPUTS/ RING ] / l ATOMS IN RING x TOTAL SPUTS 




Figure 104 - Atoca Probability - (100), 1-keV, Cu/Cu + 
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Figure 105 - Angular Dispersion - (100), 1-keV, Cu/Cu + 
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Figure 106 - Sputtering Summary - (100), 100-eV, Cu/Hg 
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Figure 107 - Atom Probability - (100), 100-eV, Cu/Hg f 
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Figure 108 - Angular Dispersion - (100), 100-e'v, Cu/'^j + 
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Figure 109 - Sputtering Summary - (100) , 200-eV, Cu/Hg 
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Figure 110 - Atom Probability - (100), 200-eV, Cu/Hq + 
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Figure 111 - Angular Dispersion - (100), 200-eV, Cu/Hg+ 
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Figure 112 - Sputtering Summary - (100) , 500-eV, Cu/Hg 
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Figure 113 - atom Probability - (100) , 500-aV, Cu/Hg + 
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Figure 114 - Angular Dispersion - (100), 500-eV, Cu/Hg + 
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Figure 115 - Sputtering Summary - (100) , 1-keV, Cu/Hg 
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Figure 116 - Atom ProDability - (100) , 1-keV, Ca/Hg + 
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Figure 117 - Angular Dispersion - (100) , 1-ke7, Cu/Pg f 
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Figure 118 - 
( 100 ) , 



Sp utter ing Mechanism Trace 
100-eV, Cu/Cu + , KS 3-3 
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^igure 119 - Spattering Mechanism Trace 
(1 00), 100-eV, Cu/Cu + , XSE-3 
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Figure 120 - Sputtering Ratio vs. Ion/Target Atom Atomic Mass Ratio 
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Figure 121 - Sputtering Ratio vs. Ion/Target Atom Atomic Mass Ratio 
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Figure 122 - Placement of 144 Random Shots 
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Figure 123 - Sputtering Summary - Random Points 
(100), 100-eV, Cu/Hg + , KSE-3, 11x4x11 
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Figure 124 - Sputtering Summary 
(100) , 100-eV, Cu/Hg+, KSE-B 
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Figure 126 - Atom Sputtering Probaoility - Grid Points 
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Figure 127 



Angular Dispersion 
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Figure 128 - Angular Dispersion - Grid Points 
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